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1.0  INTRODUCTION 


There  has  been  considerable  interest  in  the  nonlinear  electrodynamical 
characteristics  of  microparticle  composites  (Refs.  1,2).  Of  particular  interest  are 
situations  in  which  one  or  more  of  the  microparticle's  Frohlich  modes  have  been 
resonantly  excited  by  incident  laser  light.  For  such  situations  excitation  of  this  surface 
mode  amplifies  the  incident  laser  intensity  which  gives  rise  to  enhanced  nonlinear  optical 
coefficients  (Ref.  3). 

Theory  for  surface-enhanced  phase  conjugation  and  coherent  beam  combina¬ 
tion  (or  coherent  two-wave  mixing)  is  relatively  primitive,  and  this  final  report  focuses 
on  calculating  the  third-order  susceptibility,  its  magnitude,  spectral  range,  response 
time,  and  possible  saturation  characteristics. 

The  amplitude  of  the  laser  electric  field  in  the  vicinity  of  a  metallic 
microparticle  is  dramatically  amplified  at  the  surface  plasma  frequency.  This  feature  of 
microparticle  electrodynamics  arises  from  a  resonant  coupling  of  the  laser  radiation  with 
the  electron  plasma  mode  and  has  been  used  as  a  means  to  greatly  enhance  a  variety  of 
nonlinear  processes. 

In  the  theoretical  studies,  the  electrodynamic  problem  of  metallic 
microparticle  composites  has  been  treated  globally,  calculating  the  multipole  moment 
induced  upon  the  microparticles  by  the  applied  field.  The  surface-plasma  resonance  as 
well  as  the  local  field  enhancement  of  the  microparticles  within  the  quasi-hydrodynamic 
formalism  was  obtained  directly.  In  the  theory  the  local  field  enhancement  factors  come 
in  naturally.  However,  in  other  quantum  mechanical  calculations  the  enhancement 
factors  have  to  be  added  at  the  end  for  the  local  field. 

The  linear  response,  which  includes  the  linear  fields,  linear  susceptibility, 
refractive  index,  and  extinction  coefficient  with  enhancement  at  the  surface  plasma 
dipole  mode  frequency  will  be  calculated.  The  linear  response  always  co-exists  with  the 
third-order  response  at  fundamental  frequencies.  In  particular,  the  linear  absorption  is 
an  important  limiting  factor  in  the  third-order  nonlinear  processes.  The  second-order 
re^HNise,  which  includes  the  second-order  fields  and  intensity  for  the  sum  and  difference 
frequency  generation  with  enhancement  at  the  surface  plasma  dipole  and  quadrapole 
mode  frequencies  will  also  be  calculated.  The  sum  and  difference  frequency  generations 
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are  intermediate  processes  to  nondegenerate  phase  conjugation  and  coherent  two-wave 
mixing  in  microparticle  composites.  The  dc  rectified  and  second  harmonic  generations 
are  intermediate  processes  to  degenerate  and  nondegenerate  phase  conjugation  and 
coherent  two-wave  mixing  in  microparticle  composites. 

Then,  the  third-order  response,  which  includes  the  third-order  fields  and 
susceptibility  with  enhancement  at  the  surface  plasma  dipole  frequency,  will  be 
calculated.  The  contribution  of  third-order  susceptibility  arises  in  microparticle 
composites  through  the  product  of  linear  fields  and  second-order  fields  as  well  as  in  the 
triple  product  of  the  linear  fields.  The  contribution  from  the  triple  product  of  the  linear 
fields,  which  is  the  commonly  known  direct  effect,  is  not  the  major  source  of  the  third- 
order  nonlinearity.  The  contribution  from  the  second-order  fields,  referred  to  as 
cascading  effect,  is  the  dominant  one.  The  time-dependent  third-order  susceptibility  will 
be  calculated  to  explain  the  transient  behaviors  in  the  wave-mixing  for  the  microparticle 
composites. 

To  be  more  specific,  the  nonlinear  electrodynamics  of  a  microparticle 
composite  are  examined  using  a  solid-state  nonlocal  continuum  theory  (Refs.  3,4,5).  The 
problem  is  formulated  in  terms  of  the  perturbative  response  of  a  particle,  which  consists 
of  a  degenerate  electron  gas  moving  in  a  positive  uniform  jellium  (a  uniform  positively 
charged  medium  against  which  the  electrons  move)  background,  to  laser  beams  within  the 
quasi-hydrodynamic  approximation  (Refs.  4,5).  The  first-order,  second-order  and  third- 
order  dynamic  responses  of  the  particle  to  laser  light  are  extracted.  Specifically,  the 
perturbed  electron  density  and  drift  velocity  are  determined  as  well  as  the  laser-induced 
dipole  and  quadrupole  moments.  The  electrodynamic  properties  of  the  local  fields  in  the 
vicinity  of  the  microparticles  are  also  resolved.  The  macroscopic  electrodynamic 
response  of  the  medium  as  a  whole  and  results  for  th^  intensity,  polarization,  and  angular 
dependence  are  presented. 

In  experimental  studies,  the  phase  conjugation  and  coherent  two-wave  mixing 
both  in  CS2  and  gold  colloids  were  performed  using  picosecond  pulses.  The  experiments 
done  in  CS2>  which  is  a  well-studied  simple  Kerr  liquid,  were  used  to  calibrate  the  more 
complex  gold  colloid  system.  From  the  degenerate  phase  conjugation  experiments  done 
in  CS2  and  gold  colloids,  the  degenerate  third-order  susceptibility  for  gold  colloids  was 
obtained.  The  transient  time-resolved  degenerate  two-wave  mixing  experiments  were 
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performed  to  investigate  the  coherent  coupling  coefficients  in  the  time  domain  for 
nonlinear  media.  The  polarization  resolved  experiments  were  also  p>erformed  to  study 
the  anisotropy  of  the  third-order  susceptibility. 

The  experimental  values  of  the  third-order  susceptibility  for  the  gold  colloids 
are  larger  than  the  calculated  values  from  the  free  earner  contribution  within  the  quasi¬ 
hydrodynamic  formalism.  It  indicates  that  the  gold  colloids  have  other  nonlinear 

mechanisms  involved.  However,  for  the  silver  colloids  the  third-order  susceptibility 

-9 

between  the  experimental  and  theoretical  values,  which  are  10  esu  per  sphere,  are 
comparable.  This  can  be  explained  by  the  fact  that  the  plasma  frequency  of  silver  falls 
just  below  the  threshold  for  d-electron  transitions  and  so  free  electron  behavior  is 
approximately  satisfied  for  silver.  For  gold,  the  threshold  of  core  electron  transition  lies 
at  a  frequency  below  the  plasma  frequency.  The  overlapping  of  valence  electron  and 
d-electron  transitions  prevents  the  existence  of  a  frequency  range  in  which  the  free 
electron  behavior  is  dominant.  Therefore,  the  third-order  susceptibility  of  gold  colloids 
will  result  from  the  d-electron  contribution  as  well. 

The  experiment  and  the  theory  for  transient  time-resolved  and  polarization- 
resolved  coherent  two-wave  mixing  in  CS2  are  in  excellent  agreement.  Qualitative 
explanation  for  the  transient  time-resolved  and  polarization-resolved  coherent  two-wave 
mixing  in  gold  colloids  will  be  given. 

In  search  of  alternative  nonlinear  optical  materials,  there  was  a  theoretical 
investigation  of  the  one-dimensional  conjugated  polymers.  A  one-dimensional  lattice 
model  Hamiltonian  was  used  to  calculate  their  third-order  polarization  and 
susceptibility.  The  nonresonant  nonlinearity  of  the  conjugated  polymers  below  their  band 
gap  is  on  the  order  of  10”^®  esu. 

The  final  part  of  the  report  gives  conclusions  and  recommendations  for 
microparticle  composites  as  well  as  conjugated  polymers  as  a  phase-conjugation  or  image 
amplification  mediums. 
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2.0  THEORY 


Consider  a  metallic  microparticle  composite  whose  particle  size  is  <  100  A. 

This  is  small  compared  to  any  of  the  radiation  wavelengths  so  that  the  Rayleigh  limit  is 
valid.  Also,  the  quasi-static  approximation  is  used  where  the  Maxwell  equations  are 
reduced  to  the  electrostatic  equation.  If  the  incident  and  nonlinear -generated 
frequencies  are  all  sufficiently  high,  the  ionic  (or  phonon)  contribution  inside  the 
microparticles  will  be  negligible.  It  is  assumed  that  each  of  the  microparticles  in  the 
composite  consists  of  a  rigid,  positive  jellium  background  with  charge  density  en^ 
accompanied  by  a  degenerate  electron  gas.  The  dynamics  of  the  electron  gas  is 
described  by  the  quasi-hydrodynamic  theory  with  damping.  Charge  neutrality  demands 
that  the  unperturbed  electron  density  be  n^  and  in  the  presence  of  laser  radiation  the 
perturbed  electron  density  is  denoted  by  n(r,t).  Note  that  Fuchs  and  Kliewer  have  shown 
that  the  hydrodynamic  and  self-consistent  field  approximations  yield  the  same  results  at 
all  experimentally  accessible  frequencies  (Ref.  6).  Thus,  Landau  damping  can  be 
accounted  for  with  a  simple  dissipation  model  characterized  by  a  frequency-independent 
damping  constant,  v.  It  is  also  assumed  that  microparticle  size  exceeds  20  A  so  that  the 
quantum-size  effect  is  not  important  (Ref.  7). 

The  hydrodynamic  equations,  which  are  the  Euler  equation  and  the  equation  of 
continuity,  govern  the  response  of  the  electron  gas  to  incident  laser  radiation.  Here,  the 
nonlinear  response  of  a  single-sphere,  which  consists  of  a  degenerate  electron  gas, 
confined  to  a  microsphere  of  radius  a,  embedded  in  a  passive  linear  dielectric  and 
irradiated  by  laser  beams  is  examined.  The  unperturbed  electron  profile  is  taken  to  be  a 
step  function.  Although  a  step  function  does  poorly  for  a  flat  metal  surface  (Ref.  8),  it  is 
used  here  for  several  reasons.  (1)  It  is  the  standard  model  used  for  modeling  the 
nonlinear  optics  due  to  the  resonant  plasmon  responses  of  metallic  microspheres 
(Ref.  9*)  and  using  more  complex  spatial  profiles  is  computationally  prohibitive.  (2)  The 
nonlinear  plasmon  response  appears  to  involve  the  electron  fluid  deep  within  the 
microsphere,  not  just  at  the  surface.  Thus,  on  physical  grounds,  the  frequency-mixing 
properties  of  these  resonances  should  not  be  overly  sensitive  to  the  details  of  the 


*In  their  calculation,  a  circularly  polarized  light  is  used  and  no  second-order  monopole 
response. 
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electron  profile  near  the  surface.  Mathematically,  this  is  reflected  in  the  fact  that  the 
dipole  and  quadrupole  moments  require  an  integration  of  various  moments  of  the  density 
over  the  entire  volume  of  the  microsphere.  (3)  Although  other  profiles  will  give  different 
results,  these  differences  are  masked  by  the  fact  that  numerous  suspension  parameters 
(geometric  shape,  particle  size  distribution  and  metallic  purity)  will  vary  significantly 
from  the  model  and  from  sample-to-sample. 

The  physical  situation  is  depicted  in  Figure  1.  The  beam  is  incident  normal 
to  the  surface  of  microparticle  composite,  is  elliptically  polarized,  has  a  polarization 

A  A  A  ^ 

vector  =  y  +  nz  and  a  propagation  vector  kj  =  k^x  within  the  medium.  For  a 
circularly  polarized  light,  n  =  i.  The  second  laser  beam  is  oriented  at  an  angle  with 
respect  to  the  first  beam,  has  a  unit  polarization  vector  62  -  (xcose2  ^ 

propagation  vector  k2  =  k2(xs1n02  -  2C0sa2)  within  the  medium,  where  62  =  ii/2  -  e^. 
Here,  kj  =  [ej^(wj)]^^^wj/c,  with  the  dielectric  constant  of  the  host  at  the  frequency 
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Within  each  microparticle,  laser  radiation  induces  an  electron  drift  velocity, 
v<r,t)  whose  dynamics  are  governed  by  the  Euler  equation 

n(r,t){3/3t  +  V  +  v(r,t)»v}v(r,t)  =  -en(r,t)E(r,t)/m  -vp(r,t)  (1) 

Here,  v  is  the  electron  collision  frequency  and  p(r,t)  is  the  quantum  pressure  of  the 
degenerate  electron  gas,  which  within  the  Thomas-Fermi  approximation  (Ref.  9)  is  given 
by 


P(r.t)  =  [^l^^^(h2/5m)n(r.t)^^^  (2) 

The  third  term  on  the  lefthand  side  of  the  Euler  equation  is  the  convective  current, 
which  arises  from  the  fact  that  the  conduction  electrons  are  a  continuous  medium  in  this 
model.  The  first  term  on  the  righthand  side  of  Equation  1  is  the  Lorentz  force,  the 
magnetic  term  having  been  dropped  within  the  electrostatic  approximation.  The  last 
term  is  due  to  quantum  pressure,  which  gives  rise  to  a  spatial  variation  (i.e.,  a  screening) 
in  the  electron  density  near  the  microparticle  surface  in  the  presence  of  radiation.  Thus, 
the  response  of  the  electron  gas  is  nonlocal.  The  perturbed  electron  charge  density, 
n(r,t)  obeys  the  equation  of  continuity, 

an(r,t)/3t  +  v[n(r,t)v(r,t)l  =  0  (3) 

The  incident  electric  component  of  the  laser  radiation  field  EQ(r,t)  is  given  by 

expI1(lc^T  -  ujt)  ^2^2  exp[1(IC2T  -  U2t)l  +  cc  (4) 

where  Cj,  Aj,  kj  and  Uj  (j  =  1,2)  are  the  incident  unit  polarization  vector,  electric  field 
amplitude,  wave  vector  and  frequency  of  the  laser  beam.  The  total  electric  field 
amplitude,  E(r,t),  in  Equation  1  is  driven  by  the  perturbed  charge  density  inside  the 
microparticle.  Within  the  quasi-static  approximation  E(r,t)  =  -  7*(r,t)  where  •(r,t)  is  the 
electrostatic  potential.  Outside  the  microparticle,  the  Laplace  equation  holds 

v^#(r,t)  *  0  (5a) 
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Inside  the  microsphere,  the  electrostatic  potential  obeys  Poisson's  equation 

v^*(r,t)  =  4nen(r,t)/e|^  (5b) 

where  the  electrostatic  potential  is  coupled  to  Equations  1  and  3  via  n(r,t) ,  and  is 
the  dielectric  constant  of  the  rigid  jellium  background.  Equation  5a  and  b  must  be 
supplemented  with  the  usual  electrostatic  boundary  conditions,  namely  that  ♦(r,t)  and 
ea#(r,t)/ar  are  continuous  at  the  surface.  Since  the  Euler  equation  and  the  equation  of 
continuity  are  coupled  together,  only  one  additional  boundary  condition  (ABC)  is 
required,  which  is  usually  chosen  to  be  the  acoustic  condition  to  ensure  that  no  electrons 
flow  across  the  boundary;  viz.,  the  radial  component  of  the  electron  drift  velocity  will 
vanish  at  the  surface.  These  boundary  conditions  can  be  followed  without  any  difficulties 
in  determining  the  linear  response  of  the  microparticle  to  laser  light.  However,  in  the 
next  order  of  perturbation  theory,  difficulties  can  arise  in  the  presence  of  the  monopole 
response  (Ref.  9).  Specifically,  if  one  demands  that  the  second-order  radial  velocity 
vanish  at  the  microparticle's  surface,  then  the  integral  of  the  second-order  perturbed 
electron  density  over  the  volume  of  the  microparticle  will  deviate  from  zero.  In  turn, 
this  implies  that  the  calculation  does  not  conserve  charge  and  reflects  the  fact  that 
there  are  basic  inconsistencies  in  describing  the  microparticle  response  by  means  of  a 
relaxation  model  (Refs.  10,11,12).  Here,  a  different  additional  boundary  condition  is  used 
for  monopole  responses,  namely,  that  charge  is  conserved  in  every  order  of  perturbation 
theory.  If  this  is  done,  then  the  second-order  radial  velocity  will  deviate  from  zero  at 
the  microparticle's  surface.  Calculations  in  Section  3.2.3  show  that  within  a  screening 
length  of  the  surface,  the  ratio  of  the  radial  velocity  at  the  surface  to  the  maximum 
velocity  is  on  the  order  of  v where  is  the  frequency  of  the  radiation  generated  by 
the  microparticle  via  nonlinear  mixing  processes.  Thus,  so  long  as  v/u^  «  1,  this 
inconsistency  is  not  significant. 

An  examination  of  Equations  1  and  2  reveals  four  nonlinear  terms  in  the  Euler 
equation  for  the  electron  drift  velocity:  (1)  the  electron  current,  (2)  the  convective 
current  (3)  the  Lorentz  force,  and  (4)  the  quantum  pressure. 

The  electron  density,  drift  velocity,  quantum  pressure,  and  the  electrostatic 
potential  can  be  decomposed  into  the  different  possible  frequency  components.  Thus,  if 
a(r,t)  is  any  of  these  quantitites,  then 
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(6) 


a(r,t)  =  Oq  +  a^(r,uj)  exp(-iujt)  +  a^(r,w2)  exp(-iu)2t) 

+  a2(r,u_)  exp(-iu_t)  +  a2(r,u_j^)  exp(-iu^t)  +  a2(r,o) 
+  a2(r,2u)j)  exp(-2iu)jt)  +  a2(r,2t02)  exp(-2iu>2t) 

+  a2(r,u)j)exp(-iu»jt)  +  a2(r,a>2)exp(-iu)2t)  +  c.c. 


where  is  the  equilibrium  value,  are  the  linear  response  at  fundamental  frequencies 
and  the  028  are  the  sum  (u»^  =  uj  +  U2),  difference  (u)_  =  frequency,  a  rectified  dc(o) 

and  the  second -harmonic  (2uij)  contributions  for  the  second-order  response,  and  03$  are 
the  fundamental  frequency  contributions  for  the  third-order  response. 

2.1  LINEAR  RESPONSE 

Inserting  the  perturbative  expansion,  i.e..  Equation  6  into  Equations  1  through 
3,  shows  that  the  induced  linear  density  satisfies  the  following  equation 

-  q^(u)l  n^  (r,u)  *  0  (7a) 

In  Equation  7a,  q<u)  is  the  self-consistent,  frequency-dependent  screening  wave  vector, 
which  is  given  by 

q(u)  ®  Iwp  -  u)(u  +  iv)]^^^/6  (7b) 

2  1/2  7  7 

where  Up  =  (4iin^e  /wejj)  ^  is  the  plasma  frequency  and  6^  =  0.6vp^  with  Vp  the  Fermi 

velocity.  In  the  low-frequency  limit,  the  screening  wave  vector  becomes  real  and 

reduces  to  the  Thomas-Fermi  value,  qjp  =  Wp/B.  In  this  limit,  the  field  barely 

penetrates  into  the  particle  and  the  electron  density  is  uniform  everywhere,  except  in 

the  immediate  vicinity  of  the  particle's  surface.  In  the  high-frequency  limit,  i.e., 

u  >  Up  ,  q(u)  has  an  imaginary  part,  the  electron  density  is  oscillatory. 

Explicit  expressions  for  the  first-order  electron  density,  electric  field  as  well 
as  the  drift  velocity  are  given  in  Appendix  A. 
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2.2 


SECOND-ORDER  RESPONSE 


To  determine  the  quadratic  response  of  the  electron  gas  to  two  nondegenerate 
laser  beams  oscillating  at  and  U2>  insert  the  expressions  for  the  linear  electron  density, 
drift  velocity  and  electric  field  into  Equations  1  and  3  and  extract  those  terms  that  are 
oscillating  at  for  sum  (u>^)  and  difference  (w.)  frequency  generation.  This  leads  to  the 
following  equation  for  the  second-order  electron  density  response  oscillating  at  u^: 

-  q^(u+)]n2(r.(i>^)  =  -6"^vS(r;wj,±w2)  (8) 

where  the  drive  terms  S(r;u|,  ±1^2)  are  defined  by 

S(r;(aj,±u2)  =  iwjnj^(r,wj)Vj(r,±ai2)  ±  ia>2nj(r,wj)Vj(r,±u2) 

+  no[vi(r.wj)*vl  Vjlr,±u»2)  +  n^lvj(r,±tt2)*vlvj(r,wj) 

+  e/inlnj(r,«j)E^(r,±«2)  +  nj(r,±u)2)Ej(r.«j)l 

+  (26^/3nQ)v(nj(r,«j)nj(r,±u.2)l  (9) 

The  linear  response  nj(r,u),Vj(r,u)  and  E^(r,u)  are  all  enhanced  at  the  surface  mode 
frequencies  u^.  The  drive  term  S(r;(iip±(d2)  may  be  doubly  enhanced  for  and  u>2  near  u^, 
since  S(r;u£,±tti2)  is  a  bilinear  product  of  linear  responses  at  and  ii>2«  For  compact  nota¬ 
tion,  denote  S(r;u^,  ±^2)  as  S^(r).  The  first  two  terms  in  S^(r)  arise  from  the  second- 
order  electron  current,  the  next  two  are  from  the  second-order  convective  current,  the 
fifth  and  sixth  terms  arise  from  the  second-order  Lorentz  force,  while  the  last  term 
arises  from  the  second-order  quantum  pressure.  The  explicit  expression  for  the  drive 
terms  are  given  by 

r 

SjCr)  •  -i-j  lrS*(r)  +  ,sj(r)  +  ;s!(r))  (10) 

V 

with 
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S*(r)  = 

+  (16i,/45)^/2s;^^(r)nsine2Y2o(e) 

-(2i./15)^/^S;^^  (r)ncos02lY2i(e.^)  -  Y2_i(e.^)l 
+  i(2it/15)^/^  S;2i+(r)sine2lY2i(0.^)+  Y2_i(e.*)l 
-i(2T,/15)^/^S;^22<''^‘^°^®2^^22^®**)  Y2_2(e.«)l  (Ha) 

S*(r)  =  (4Tr/45)^/2  ^  |_  Y2^(e) 

-  (2h/15)^/2  S;2i.(r)n  cose2  [Y2ie.^)  -  Y2,ie.^)] 

+  i  (2,r/15)^^2  S;2i+(r)  sin  62  ^  [Y2i(e.*)  +  Y2.1  (e.*)l 

-  i  (,./30)^/2  S*22(^)  cos  02  |j  lY22(e.*)  -  Y2.2(e.o)l  (Hb) 

Sj(r)  -  -  (2„/15)^/2  s±  (r.  ±  n  cos  82  ^  lY2i(e.^) 

-  Y2.i(e.^)l  +  1(211/15)^/2  S*(r.«j.  ±  0*2)  sin  02^1^ 

IY2i(e.*)  +  Y2_i(o.4)1  -1  (ir/30)^/2  ^^,5  62 

^^22^®*^^  “  Y2_2(0f^)l  (11c) 

where  the  different  components  (r)  are  given  in  Appendix  B  (Equations  B.l).  Also 
explicit  expressions  for  the  second-^der  electric  field  as  well  as  the  electron  density  and 
velocity  are  contained  in  Appendix  B. 
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Similar  expressions  can  be  obtained  for  the  dc  rectified  and  second-harmonic 
contributions. 

2.3  THIRD-ORDER  RESPONSE 

To  determine  the  third-order  response  of  the  electron  gas  to  two  nondegener¬ 
ate  laser  beams  oscillating  at  and  U2>  insert  the  expressions  for  the  linear  and 
quadratic  electron  density,  drift  velocity,  and  electric  field  into  Equations  1  and  3  and 
extract  those  terms  that  are  oscillating  at  (d|  and  ^2  fundamental  frequency  generation. 
This  leads  to  the  following  equation  for  the  third-order  electron  density  response 
oscillating  at  uj: 

-  q^(w^)ln3(r,u^)  =  -6‘^vT(r;u^)  (12) 

where  i  =  1  or  2  and  the  drive  terms  T(r;iit|)  are  defined  by 

T(r;«j^)  =  T^^^(r,u»2  +  «_)  +  +  0)  +  T^^^(r,-W2 

+  -  uj  2«j)  +  a)2>  “  «2) 

+  (i>^,  <iij,  -  wj^)  (13a) 

T(r;«2)  -  T^^)(r,«j  -  «_)  +  T^^^r,«2  +  0)  +  T^^^(r.-«j  +  u^) 

+  T^^^(r,  -  m2  +  2012)  u»2,  -  wj^) 

+  (r,  U2,  -  «2)  (13b) 

where 

T(2)(r,  -f  0)  ■  +  («l  -  «j))  +  T^^^(r,  +  (u2  -  **2^^ 

(13c) 
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In  Equation  13a,  the  first  cascade  term  is  given  by 

T^^^(r,u)2+w_)  =  iu2nj(r,ui2)v2(r,u_)  +  iw_n2(r,a)_)Vj(r,u)2) 

+  no(Vi(r,w2)-v]  V2(r.w_)  +  nQlv2(r,<.._)-v]Vj{r.a)2) 

+  e/m[n^(r,a»2)E2(r,a)_)  +  n2(r,a»_)E^(r.u)2)  1 
+  (2e^/3nQ)y(nj(r.a.2)n2(r,u)_)I 

(14a) 


and  the  first  direct  term  is  given  by 

T^^^(r,u»j,iu2.-«2^  '  nj(r.wj)  Ivj(r.(a2)  •  v)  v^(r,  -u»2)  +  c.c.) 

+  n^(r.«2)  lv^(r.  -^2)  * 

+  (vj(r,uj)-v)  Vj(r,  -<02)1 

+  n  )  l(j  (r.wj)  •  j)  v  (r.uj) 

+  v^(r,a)2)*v)vj)r,wj)]  (14b) 

and  similar  expressions  can  be  obtained  for  the  other  terms  in  Equations  13. 

The  linear  responses  n|(r,U|),  Vj(r,U|)  and  E|(r,up  are  all  enhanced  at  the  surface 
mode  frequencies  u^.  The  second-order  responses n2(r,w_),  V2(r,u>_),  and  E2(r,u»_) 
may  be  doubly  enhanced  for  and  U2  near  u^.  The  first  two  terms  in  Equation  14a  arise 
from  the  electron  current,  the  next  two  are  from  the  convective  current,  the  fifth  and 
sixth  terms  arise  from  the  Lorentz  force,  while  the  last  term  arises  from  the  second- 
order  quantum  pressure.  All  the  terms  appearing  in  Equation  14b  arise  from  the 
convective  current.  The  explicit  expression  for  the  drive  terms  can  be  worked  out. 
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In  Equations  13,  the  driving  terms  and  come  from  the  products  of  the 
second-order  fields  and  the  linear  fields  as  well  as  the  triple  products  of  the  linear 
fields.  The  driving  terms  and  can  be  further  broken  down  into  self-action  terms 
eind  cross-action  terms.  The  self-action  term  and  cross-action  term  in  T^^^  at  are  T^^^ 
(r,upu)|,-o)|)  and  T^^\r,ui  ,W2,-U2^»  respectively.  The  self-action  terms  in  T^^^  are 
T^^\r,u)j+(u|-un)  and  T^^V»-«|+2a>|),  which  come  from  the  product  of  the  self-rectified  and 
self-second-harmonic  fields  and  linear  fields.  The  cross-action  terms  in  T^^^  at  ui  are 
T'  Hr,w2+u_),  T'  ^(r,o»|+(u2-w2))»  T'  '(r,-u2+w^),  which  are,  respectively,  from  the 
product  of  the  difference  frequency,  cross-rectified,  and  sum  frequency  fields  and  linear 
fields.  The  same  classification  is  applied  for  T^^^  and  T^^^  at  u)2.  The  terms  in  T^^^ 
resulting  from  the  products  of  the  difference  and  sum  frequency  fields  and  linear  fields 
are  also  known  as  "coherent  coupling"  terms,  since  they  are  important  for  the  coherent 
waves. 


Accordingly,  all  the  third-order  responses,  which  include  the  third-order 
susceptibility,  polarization,  density,  electric  field,  and  velocity,  can  be  classified  from 
their  corresponding  driving  terms.  That  is,  there  will  be  self-action  and  cross-action 
effects  in  third-order  responses  resulting  from  all  possible  combinations  of  second-order 
fields  and  first-order  fields  as  well  as  the  triple  products  of  the  linear  fields. 
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3.0  ELECTRODYNAMICS  OF  MICROPARTICLE  COMPOSITES 


In  this  section,  the  macroscopic  linear  and  nonlinear  electrodynamic  responses 
of  a  very  dilute  microparticle  composite  are  investigated.  The  electrodynamic  response 
of  the  microparticle  composite  is  determined  by  the  total  field  generated  in  the  entire 
composite.  For  very  dilute  microparticle  composites,  each  microparticle  is  in  the  far 
field  of  all  the  others  and  the  single  scattering  approximation  is  valid.  Hence,  the  total 
field  of  the  composite  is  simply  the  sum  of  the  far  field  generated  by  each  particle. 

The  linear  and  nonlinear  electrodynamic  responses  of  the  microparticle 
composite  are  determined  by  the  induced  linear  and  nonlinear  polarization  of  the 
microparticles.  The  coherent  part  of  the  electrodynamic  response  is  determined  by  the 
phase-matched  amplitude  of  the  microparticles,  viz.  the  "optical  theorem"  (or  the 
"Ewald-Oseen  extinction  theorem"). 

3.1  LINEAR  ELECTRODYNAMICS 

3.1.1  Scattering  Amplitude  in  the  Far  Zone 

In  this  section,  the  first-order  electromagnetic  fields  in  the  radiation  zone, 
which  are  generated  by  a  single  microsphere  are  determined.  A  microsphere  at  point 
r*  s  (x,y',z')  scatters  the  two  incident  beams  to  generate  a  first-order  field  )  at 

the  incident  (up  frequency  at  point  r  =  (x,y,z),  where  r  »  a  in  the  radiation  zone.  The 
electromagnetic  field  associated  with  an  oscillating  dipole  of  the  microsphere  located  at 
r*,  is  given  by 

E^^hr.r'.u^)  »  expll4(r.r')]f{^^e'.*')/!r-r'|  (15) 

where  the  phase  4(r,r*)  is  given  by 

♦(r.r*)  «  k^|r-r'|  +  k^-r'  (16) 

and  kj  s  [ej^(uj)]^'^*wj/c  with  i  =  1  or  2,  the  integration  is  over  the  volume  V  =  4it/3  a^  of 
the  ^>here  at  r*.  In  the  quasi-static  limit  where  the  radiation  wavelength  x  »  a,  the 
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phase  factors  in  Equation  16  remain  nearly  constant  within  a  given  sphere.  Hence,  the 
scattering  amplitude  in  Equation  13  is  given  by 

=  k^(nxP^^hxn  (17) 

where  n  is  the  unit  vector  in  the  direction  (r-r*)»  and  the  vector  p^^^  is  the  first-order 
dipole  moment  tensor  induced  by  the  laser  light  in  a  microparticle  located  at  the  origin, 
cose'  =  (z-2')/r  and  tan*'  =  (y-y')/(x-x').  Note  that  the  dipole  moment  tensor  p^^^  can  be 
extracted  from  the  electrostatic  field  in  the  near  zone  for  a  small  scatter.  For  a 
microsphere  irradiated  by  laser  beam  at  (i>|,  the  radiation-induced  dipole  moment  is  given 
in  Appendix  A  (Equation  A.7).  Finally,  the  scattering  amplitude,  fj^\e',*'),  for  the 
electric  field  in  the  radiation  zone,  arising  from  a  single  microparticle  is 


fP^(0'»*')  =  kja^P^^^(wj^)  (6  Icose  sin*  -  nsine]  -  *  cos*} 

(18a) 

f^^^(e',*')  =  k|a^  P^^^(«2)  (0  [cos02cose  cos*-s1ne2s1nel  -  *  cose2Sln*} 

(18b) 

where  the  incident  frequency  is  at  and  the  polarization  P^^\u)|)  is  defined  in 
Appendix  A  (Equation  A.8). 

The  magnetic  field  components  are  then  given  by 


B^^^(r.r'.u^) 


n  X 


(19) 


so  that  the  total  intensity  Ig(r,r',W|),  for  a  single  microparticle  at  r'  is, 
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(20) 


3.1.2  Macroscopic  Electromagnetic  Fields  and  Intensities 

In  this  section,  the  total  linear  electromagnetic  fields  and  intensities  arising 
from  a  tenuous  collection  of  microp>articles  randomly  distributed  in  a  passive  linear 
dielectric  host  are  examined.  In  general,  a  collection  of  microparticles  will  scatter 
radiation  with  random  phase  with  respect  to  one  another.  The  resultant  total  linear 
electromagnetic  field  from  a  microparticle  composite  will  be  the  sum  of  the  coherent  (or 
averge)  field  and  the  incoherent  (fluctuating)  field  (Ref.  13).  The  phase  of  the  field  will 
fluctuate  randomly  for  arbitrary  directions,  so  that  the  volume-averaged  (or  ensemble- 
averaged)  incoherent  field  is  zero.  However,  along  a  particular  phase-matched  (or 
coherent)  direction,  the  different  particles  ail  scatter  in  phase  with  one  another  and  the 
resultant  average  field  is  coherent.  This  coherent  field  will  be  discussed  in 
Section  3. 1.2.1.  The  average  intensity  for  a  microparticle  composite  is  the  sum  of  the 
coherent  and  incoherent  intensities  and  will  be  examined  in  Section  3. 1.2.2.  A  tenuous 
microsphere  composite  with  the  particles  sufficiently  far  apart  so  that  they  do  affect 
each  other  will  be  considered.  Then,  the  far-field  and  single  scattering  approximations 
are  valid  and  the  total  field  at  a  point  r  is  the  sum  of  the  far  field  from  each  microsphere 
located  at  point  r',  as  depicted  in  Figure  1. 

3.1.2.1  Coherent  Field:  Linear  Dielectric  Function.  The  linear  coherent  field  at  is 
given  by  the  integral  (or  summation)  of  the  far  field  of  each  sphere,  which  is  given  in 
Equation  13,  over  the  entire  composite,  i.e.. 


E^^^(r,wi)  *  Jydr'E^^^Cr.r'.w^)  ,  (21a) 

1  L  1  L  ^  L 

-  N  dx'J^  ^dy'/f  ^  dz*  expIU(r,r')If{^)(e',^')/|r-r'| 

-jL 

(21b) 
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where  the  composite  is  assumed  to  reside  within  a  slab  of  dimensions  L  x  x  L_,  where 
L  «  and  LyL^  -►“>  but  L  remains  finite  and  N  is  the  microsphere  number  density. 

The  prefactors,  which  include  the  scattering  amplitude  f .  and  the  propagator  l/|r-r'| 
vary  slowly  with  r'  in  comparison  to  the  rapidly-varying  exponential  factor,  exp[i*(r,r')]. 
The  concern  is  with  the  limiting  case  in  which  r  »  L.  To  evaluate  the  integral  for  the 
y'-,  and  z'-coordinate  first  for  the  slab,  the  method  of  steepest  descents  (Refs.  15,16)  will 
be  used.  The  stationary  phase  point  for  the  microsphere  is  located  at 

The  coordinates 


85  =  0^  (22a) 

♦s  =  0  (22b) 


and  9j  are  the  incident  angles  at  wj. 

Equations  22a  and  b  indicate  that  the  phase-matched  direction  for  the  linear 
response  is  in  the  forward  direction. 

Then,  the  coherent  field  at  m2  in  Equation  21  is  given  by 


r(l) 


(r.u.)  =  -  efL 


sine. 


65.0) 


(23) 


where  e  is  unit  polarization  vector  of  the  coherent  electric  field  at  point  r  far  outside 
the  slab  and  f  =  N4iia  /3  is  the  volume  fraction  of  microspheres.  In  Equation  23,  the 
forward  scattering  amplitude  per  unit  volume  f  j  results  from  evaluation  of  the  integral  in 
Equation  21  at  the  stationary  phase  point  and  is  given  by 


T{^^(e^|0s.o)  -  2ir1k^P^^)(«^) 


(24) 
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This  general  relation  of  the  linear  coherent  field  given  in  Equation  23,  which  is  embodied 
in  the  forward  scattering  amplitude,  is  also  called  the  "Optical  Theorem"  or  the  "Ewald- 
Oseen  extinction  theorem." 

If  the  microparticle  composite  of  thickness  L  is  replaced  by  a  homogeneous 
dielectric  slab  of  same  thickness  with  effective  linear  susceptibility  xp  then  the  electric 
field  at  point  r  outside  the  slab  is  given  by 

=  2Ti1k^Lxj^((i)j)A^  /  e^(u^)sin0.  (25) 

where  the  uj,  Aj  amd  Oj  are  the  incident  frequency,  field  amplitude  and  angle, 
respectively. 

Equating  Equations  25  and  23,  the  effective  linear  dielectric  function  or 
susceptibility  of  the  microparticle  composite  with  volume  fraction  f  is  given  by 

e(ui)  =  (26a) 

(26b) 

where  the  microparticles  are  embedded  in  a  passive  linear  dielectric  material  with  the 
dielectric  constant  being  (uj). 

3.1.2.2  Total  Intensity;  Linear  Scattering.  The  total  intensity  from  a  slab  of 
microparticle  composite  of  thickness  L  is  given  by 

1  L 

=N  dx'  X;  dy'  X;  dx'  l|^)(r,r'.u,^)  (27) 

1 

where  N  is  the  microparticle  number  density  and  is  given  in  Equation  20.  The  total 
intensity  can  be  decomposed  into  two  parts,  viz,  the  coherent  and  incoherent  intensi¬ 
ties.  The  coherent  part  is  given  by 
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(28a) 


=  -T— kf  |P***(»,)i' 

sin^e^  ^  ^  ’ 


(28b) 


where 


=  r  X  and  is  given  in  Equation  23.  The  incoherent  part  is  given  by 

r  !•-  |fS^^(0'.4'') 

'f  ’  '  57 "  l'"''  n  “1''  c 


(29a) 


where  the  x',y',  and  z*  integrations  can  be  easily  evaluated,  and  the  incoherent  intensity 
is  given  by 


(29b) 


,  fl 

which  is  less  than  1  percent  of  the  incident  intensity  for  volume  fraction  f  <  10  . 

3.1.3  Numerical  Results  and  Discussion 

The  microparticles  to  be  considered  in  the  numerical  calculation  are  silver 
spheres  with  radius  a  >  20  A,  where  the  quantum  size  effect  is  not  important  (Ref.  14). 
The  background  dielectric  constant  e|j  and  the  plasma  frequency  Wp  of  the  silver  spheres 
are,  respectively,  3.378  and  3.8  eV.  The  damping  constant  v  inside  the  silver  sphere  is 
modified  by  the  free  path  effect  (Ref.  14)  and  is  given  by 

Vp 

V  -  Vb  +  /  (30) 
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where  vg  =  0.01  Up  is  the  bulk  damping  constant  and  Vp  is  the  Fermi  velocity  of  the 
electrons. 

For  comparison,  the  linear  resptonses  are  plotted  in  Figure  2  for  a  single  silver 
sphere  with  a  =  100  A  in  vacuum,  which  is  irradiated  by  0.3787  ym  laser  light  of  1  MW 
intensity.  The  0.3787  ynri  laser  light  can  excite  the  surface  dipole  mode  of  the 
aforementioned  silver  sphere.  Indeed,  as  shown  in  Figure  2(b),  the  electric  field  inside 
the  sphere  is  enhanced  by  the  surface  dipole  resonant  denominator  D(ug),  whose 
magnitude  is  about  10,  and  is  nearly  ir/2  out  of  the  phase  with  respect  to  the  incident 
electric  field.  Figure  2(c)  depicts  the  radial  velocity  throughout  the  microsphere. 

In  Figures  3  through  10,  the  linear  dielectric  properties  of  the  microsphere 
composites  are  shown.  In  particular,  the  calculations  are  done  for  the  parameters 
appropriated  for  silver  in  water  or  glass  substrate.  The  enhancement  near  the  surface 
plasma  mode  of  the  spheres  ise  clearly  shown  in  all  the  linear  dielectric  properties,  e.g., 
Figures  7  through  10.  The  higher  enhancement  factors  will  appear  in  the  second-order 
and  third-order  nonlinear  dielectric  properties  from  the  same  local  enhancement  field. 

All  the  nonlocal  nature  of  the  physical  quantities  will  approach  the  local  limit 
for  the  large  qja.  In  Figure  5  the  particle  size  dependence  of  the  linear  dielectric 
function  of  the  microsphere  composite  is  shown.  The  dielectric  function  has  an 
enhancement  from  the  nonlocality  of  the  small  sphere  limit  and  approaches  the  local 
limits,  which  are  indicated  by  the  arrows  in  Figure  5,  for  large  sphere  limit.  However, 
for  extremely  small  spheres  the  dielectric  function  decreases  and  approaches  the  pure 
host  limit.  Note  in  Figure  5,  the  nonlocal  nature  is  more  pronounced  for  the  smaller 
damping  parameter,  i.e.,  v  =  0.02  Up. 

3.2  SUM  AND  DIFFERENCE  FREQUENCY  GENERATIONS 
3.2.1  Scatto'ing  Amplitude  and  the  Radiation  in  the  Far  Zone 

In  this  section,  the  second-order  electromagnetic  fields  in  the  radiation  zone, 

which  are  generated  by  a  single  microsphere  are  determined.  A  microsphere  at  point 

^  =  (x,y',2')  scatters  the  two  incident  beams,  to  generate  a  second-order  field 

(2) 

'  (r,w^)  at  the  sum  (w^)  or  difference  {aj  frequency  at  point  r  =  (x,y,z),  where  r  »  a 
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(a)  Electron  density  profile 
n^  =  3.08  X  10‘®  cm“3. 


(b)  Radial  electric  field. 


Figure  2.  Linear  responses  of  a  single  silver  sphere  versus  r/a,  a  =  100  A.  The 

incident  laser  light  is  1  MW  at  0.3787  ym.  The  magnitude  of  the  incident 
electric  field  is  =  0.92  x  10^  stat  V/cm. 
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(a)  Surface  plasma  resonant  frequency  versus  sphere  radius.  The  host  material  is  either 
water  (e^  =  1.7689)  or  glass  =  2.23). 


(b)  Surface  plasma  resonant  frequency  versus  host  dielectric  constant  for  radius  being 
30  A  or  100  A.  The  damping  parameter  (v)  is  either  0.02  01^  or  0.1  Uq.  The  plasma 
frequency  is  3.8  eV.  The  background  dielectric  constant  ^|,)  of  the  sphere  is  3.378. 

Figure  3.  Surface  plasma  resonant  frequency  versus  sphere  radius. 
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(a)  The  local  field  enhancement  factor  at  the  surface  plasma  resonant  frequency  versus 
sphere  radius. 


(b)  The  local  field  enhancement  factor  at  the  surface  plasma  resonant  frequency  versus 
the  host  dielectric  constant.  The  parameters  for  the  ^>herical  particle  are  the  same  as 
used  in  Figure  3. 

Figure  4.  The  local  field  enhancement  factor. 
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u  /  (®co)Xiitr]  ay 


Figure  6. 


(a)  Real  part 


(b)  Imaginary  part 


The  host  dielectric  constant  dependence  of  the  real  and  imaginary  parts 
of  the  dielectric  function  of  the  microspheres  for  radius  30  a  or  100  A. 
The  parameters  for  the  spherical  particle  are  the  same  as  used  in 
Figure  3. 
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lm[4jtXi(“)/  Re[47i  Xi(w)/f] 


Figure  7. 


/  a>p 


(a)  Real  part 


0)  /  0>p 

(b)  Imaginary  part 


The  frequency  dependence  of  the  real  and  imaginary  parts  as  well  as 
absolute  value  of  dielectric  function  of  microspheres  in  water.  The 
values  of  the  dielectric  function  are  normalized  by  the  volume  fraction  f. 


26 


I  I  /  (o))Xirt^  I  U  /  (co)'X  Jit'jau 


Figure  8. 


©/©p 

(c)  Absolute  value 
Figure  7.  (Concluded) 


<*>  /  CDp 

(a)  Real  part 


The  fret^ncy  dependence  of  the  real  and  imaginary  parts  as  well  as 
absolute  value  of  the  dielectric  function  of  microspheres  in  glass. 
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(b)  Imaginary  part 


<c)  Absolute  value 


Figure  8.  (Concluded) 
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Figure  9. 


£0/  m 

p 


(a)  Refractive  index 


©/(Bp 

(b)  Extinction  coefficient 


The  frequency  dependence  of  the  refractive  index  and  the  extinction 
coefficient  in  water.  The  values  are  normalized  by  the  volume  fracticm 
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(D/tOp 


(a)  Refractive  index 


0.60  0.65  0.70  0.75  0.80  0.85  0.90  0.95 

©/(Dp 


(b)  Extinction  coefficient 


Figure  10.  The  frequency  dependence  of  the  refractive  index  and  the  extinction 
coefficient  in  glass. 


in  the  radiation  zone.  The  electromagnetic  field  aissociated  with  an  oscillating 
quadrupole  of  the  microsphere  located  at  r',  is  given  by 


(r.r- 


,u^)  =  exp[i<t.(r,r')lf^^^(e'.4.')/|r-r' 1 


(31) 


where  the  phase  ♦(r,r')  is  given  by 


♦  (r,r‘)  =  k^lr-r'l  +  (k^±k2)T' 


(32) 


1/2  ^ 
and  k^.  =  [ef,(w+)]  '  u^/c  the  integration  is  over  the  volume  V  =  4it/3  a^  of  the  sphere  at 

r*.  In  the  quasi-static  limit  where  the  radiation  wavelength  X  »  a,  the  phase  factors  in 

Equation  32  remain  nearly  constant  within  a  given  sphere.  Hence,  the  scattering 

amplitude  in  Equation  31  is  given  by 


)  »  -i(k^3/2)  (nxQ(n)/3)xn 


(33) 


where  the  vector  Q(n)  is  defined  by 


Q  (n  )  =  lQ  n  ^34) 

Qjjiy  is  the  second-order  quadrupole  moment  tensor  induced  by  the  laser  light  in  a 
microparticle  located  at  the  origin,  cose'  =  (2-z')/r  and  tan^'  =  (y-y')/(x-x').  Note  that  the 
quadrupole  moment  tensor  Q  can  be  extracted  from  the  electrostatic  field  in  the  near 
zone  for  a  small  scatter.  For  a  microsphere  irradiated  by  two  laser  beams  oriented  at  an 
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angle  of  ©2  with  respect  to  each  other,  the  radiation-induced  quadrupole  moment  tensor 
is  given  by 


(35) 


where  Explicit  formulas  for  are  given  in  Appendix  B  (Equation  B.17). 

Finally,  the  scattering  amplitude,  (©; «4>' )  *  for  the  electric  field  in  the  radiation 
zone,  arising  from  a  single  microparticle  is 


(e't*')  =  1kJa®I2iter^/(3e^(uj^)  +  2ejj)l{©I-sin02Sin2e'U2Q/2 


+  cos02Cos20'cos4iU2j1  -♦cos02cos0'sin4''u|j} 

(36) 


Note  that  in  the  collinear  and  copolarized  limit  (©2  =  k/2),  the  scattering  amplitude  has  a 
component  only  in  the  e-direction.  The  functions  are  dimensionless,  are  related  to 
weighed  averages  of  the  various  components  of  the  charge  density  over  the  sphere  and 
are  defined  in  Appendix  B.  Numerical  results  for  U^q  as  a  function  of  particle  sizes  will 
be  discussed  in  Section  3.2.3. 

The  magnetic  field  components  .w^)  are  then  given  by 


(r.r'.w^) 


n  X 


(r.r'.w^.) 


(37) 


so  that  the  total  intensity 


(r.r' 


,u^)  •  for  a  single  microparticle  at  r*  is, 
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(38) 


where  Q(n)  is  defined  in  Equation  33. 

3.2.2  Macroscopic  Electromagnetic  Fields  and  Intensities 

In  this  section,  the  total  second-order  electromagnetic  fields  and  intensities 
arising  from  a  tenuous  collection  of  microparticles  randomly  distributed  in  a  passive 
linear  dielectric  host  are  examined.  In  general,  a  collection  of  microparticles  will 
scatter  radiation  with  random  phase  with  respect  to  one  another.  The  resultant  total 
second-order  electromagnetic  field  from  a  microparticle  composite  will  be  the  sum  of 
the  coherent  (or  averge)  field  and  the  incoherent  (fluctuating)  field  (Ref.  13).  The  phase 
of  the  field  will  fluctuate  randomly  for  arbitrary  directions,  so  that  the  volume-averaged 
(or  ensemble-averaged)  incoherent  field  is  zero.  However,  along  a  particular  phase- 
matched  (or  coherent)  direction,  the  different  particles  all  scatter  in  phase  with  one 
another  and  the  resultant  average  field  is  coherent.  This  coherent  field  will  be  discussed 
in  Section  3.2.2. 1.  The  average  intensity  for  a  microparticle  composite  is  the  sum  of  the 
coherent  and  incoherent  intensities  and  will  be  examined  in  Section  3.2.2.2.  A  tenuous 
microsphere  composite  with  the  particles  sufficiently  far  apart  so  that  they  do  affect 
each  other  will  be  considered.  Then,  the  far-field  and  single  scattering  approximations 
are  valid  and  the  total  field  at  a  point  r  is  the  sum  of  the  far  field  from  each  microsphere 
located  at  point  r',  as  depicted  in  Figure  1. 

3.2.2.1  Coherent  Field;  Coherent  Generation.  The  second-order  coherent  field  at  u>^ 
is  given  by  the  integral  (or  summation)  of  the  faur  field  of  each  sphere,  which  is  given  in 
Equation  31,  over  the  entire  composite,  i.e.. 


“  J’ydr'E^^^(r,r',o»^)  (39a) 

1  L  1  L  ^  I 

-  N  / j  dx'/j  ^  dy'/j  *  dz'  expI1*(r,r‘)lf^^^(0',*')/jr-r’ I 
“7  “7  “-y  ■?  *-z  (39b) 
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where  the  composite  is  assumed  to  reside  within  a  slab  of  dimensions  L  x  Ly  x  where 
L  «  LytL^  and  LyL2  -►“»  but  L  remains  finite  and  N  is  the  microsphere  number  density. 
The  prefactors,  which  include  the  scattering  amplitude  and  the  propagator  l/|r-r'| 
vary  slowly  with  r*  in  comparison  to  the  rapidly-varying  exponential  factor,  exp[i*(r,r')]. 
The  concern  is  with  the  limiting  case  in  which  r  »  L.  To  evaluate  the  integral  for  the 
y*-,  and  z*-coordinate  first  for  the  slab,  the  method  of  steepest  descents  will  be  used 
(Ref.  16).  The  stationary  phase  point  for  the  microsphere  is  located  at 

The  coordinates 


=  Ix-x' l/s1ne^ 

(40a) 

o 

II 

+1 

-e- 

(40b) 

cose^.  =  ±k2Cose2/k+ 

(40c) 

and  are  the  coherent  angles  at  wj.  An  examination  of  Equation  39  reveals  that  there  is 
always  a  stationary  phase  point  for  the  sum-frequency  («^)  wave  or  the  difference- 
frequency  («_)  wave,  if  >  k2Cose2*  Thereby  it  restricts  the  range  of  angles  in  which 
the  second  laser  can  be  aligned  for  coherent  sum  or  difference  frequency  wave 
generation. 

The  last  integral  over  the  x'-coordinate  in  Equation  39b  is  over  a  finite  region 
and  will  give  rise  to  the  phase-matching  term.  Then,  the  coherent  field  in  Equation  38  is 
given  by 


^  *  '  (41) 

where  the  unit  polarization  vector  of  the  coherent  electric  field  at  point  r  far  outside  the 

AAA 

slab  is  ■  xcose^  -  zslne^  ,  jQ(z)  =  sinz/z  is  the  usual  phase-matching  term  and 
f  =  N4ira^/3  is  the  volume  fraction  of  microspheres.  In  Equation  41,  the  scattering 
amplitude  per  unit  volume  is  evaluated  at  the  stationary  phase  point  and  is  given  by 
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f 02 1  ®± “  i|er+lc|a^  I  cose2COs20^U2 j+2/3s i  ne2Si  n2©+U2Q  1  /  (3e (a)+)+2ejj) 

(42) 

while  the  phase  of  the  wave  ^^(r)  is 

♦+(»“)  =  k^Isine^x  +  cos0^.2l  (43) 


The  phase  factor  A^L/2  arises  from  phase  mismatch  between  the  polarization  driven  by 
incident  waves  and  nonlinear  wave,  which  is  given  by 


=  kj  ±  k2s1n02-  k^slne^ 


(44) 


An  examination  of  the  scattering  amplitude  evaluated  at  the  stationary  phase  point 
reveals  that  there  is  no  coherent  response  if  the  lasers  are  collinear  and  copolarized,  that 

IS,  ©2  “ 

3.2.2.2  Total  Intensity;  Second-Order  Scattering.  The  total  intensity  for  the  sum  and 
difference  frequency  generation  from  a  slab  of  microparticle  composite  of  thickness  L  is 
given  by 


*  ^ 

-  H  dx'  i;  dy'  n  dx'  (M) 

where  N  is  the  microparticle  number  density  and  is  given  in  Equation  38.  The  total 
intensity  can  be  decomposed  into  two  parts,  viz,  the  coherent  and  incoherent  intensi¬ 
ties.  The  coherent  part  is  given  by 
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(46a) 


I^^^(r.a)±)  =  I;-  X  •  r 


9  9  0  r  Ik.lVe^lr.r  +  p 

=  A2  - =7?  ^c(®2)  |jo(A+L/2)|2 

®  3ej^(o.^)+2ej2  2  °  - 


(46b) 


where 


±  i2 


1^(82)  =  |cos02Cos2e^U2j^  +  2/3sin62Sin20^U2Q| 


(46c) 


=  r  X  and  E^^^  is  given  in  Equation  41.  The  incoherent  part  is  given  by 


c  i"-  r-  r-  Ifi 

■?  '- 


|r-  r-| 


.|2 


(47a) 


where  the  x',y',  and  z*  integrations  can  be  easily  evaluated*  and  the  incoherent  intensity 
is  given  by 


-  fL 


3  A 


(47b) 


where 
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(47c) 


2  2 
1+(02)  =  cos^02U^j  +  8/9  s1n^02U|Q 


3.2.3  Numerical  Results  and  Discussion 

The  microparticles  to  be  considered  in  the  numerical  calculation  are  silver 
spheres  with  radius  a  >  20  A,  where  the  quantum  size  effect  is  not  important  (Ref.  14). 

All  the  parameters  for  silver  spheres  are  given  in  Section  3. 1 .3. 

Figures  11  through  13  depict  the  second-order  responses  of  the  silver  sphere 
for  difference  and  sum  frequency  generation.  To  illustrate  the  surface  enhanced  nature, 
the  frequency  of  the  first  incident  laser  light  is  always  tuned  to  the  surface  dipole 
resonant  frequency,  viz.,  01^=  in  order  to  maximize  the  enhancement  of  the  drive 
terms  for  the  second-order  fields.  The  frequency  of  the  second  laser  light  is  tuned  to  the 
vicinity  of  u^,  to  optimize  the  second  enhancement  factor  in  the  drive  terms.  Also,  the 
second  leiser  light  may  be  tuned  to  the  frequency  it>2  =  U5  -  <0q.  The  intensity  of  both 
incident  light  beams  are  1  MW/cm^. 

All  the  responses  in  Figures  11  through  13  for  the  small  difference  frequency 
and  the  sum  frequency  at  the  surface  quadrupole  resonance  are  well  screened,  since  the 
frequencies  are  below  the  plasma  frequency.  And  the  responses  in  Figures  11  through  13 
for  the  sum  frequency  generation  near  2u^  are  highly  oscillatory. 

For  quadrupole  (f,  =  2)  responses,  the  total  charge  is  zero  and  the  acoustic 
condition  can  be  imposed  as  an  additional  boundary  condition  (ABC).  However,  for 
monopole  (l  =  0)  responses,  one  must  ensure  the  charge  conservation  and  the  ABC  is  the 
charge  neutrality  instead  of  the  acoustic  condition.  Thus,  the  acoustic  condition  may  be 
violated,  which  implies  an  electron  flow  across  the  boundary  of  the  sphere. 

In  Figure  13,  the  acoustic  condition  is  indeed  violated  for  1  =  0.  For  the 
difference  frequency  generation  shown  in  Figure  13(a),  the  radial  velocity  at  the  surface 
is  nearly  equal  to  the  maximum  velocity  inside  the  sphere,  where  v/w_  -  1.  For  the  sum 
frequency  generation  near  2U5  shown  in  Figure  13(c),  the  radial  velocity  at  the  surface  is 
nearly  zero,  where  v/w^  ~  0.01.  Therefore,  the  ratio  of  the  radial  velocity  at  the  surface 
to  the  maximum  velocity  inside  the  sphere  is  on  the  order  of  v/m^  for  e  =  0.  Also,  the 
1  =  2  second-order  responses  in  Figures  11(e),  12(e),  and  13(e)  are  indeed  enhanced  with 
respect  to  the  1  =  0  counterparts  in  Figures  11(b),  12(b),  and  13(b). 
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(a)  (t>  1  =  i  =  0,  m  =  0; 

difference  frequency  genera¬ 
tion  with  a»_  =  0.01  Wp. 


(b)  u  1  =  Ugy  I  =  Oy  m  =  0; 
for  sum  frequency  generation 
with  =  lUq  =  0.3642  pm. 


0.80  0.85  0.90  0.95  1.00 


Figure  11.  Second-ordor  dimensionless  electron  density  profile  n2  versus 

r/a»  a  »  100  A.  * 
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0.02 


(d)  i  =  2y  m  =  0;  for 
difference  frequency 
generation  with 

oi_  =  0.01  Up. 


(e)  i  =  2,  m  =  0;  for 
sum  frequency  generation 
with  u^  =  Uq  =  0.3642  ym. 


n2,20*''“-*  0 

-0.02 
-0.04 

0.8 


(e) 


n2,20*'^'“q* 


0.6  - 
0.4  - 
0.2  - 
0  ■ 

_ i_ 

0  0.2 


0.4 


_i _ I _ i_ 

0.6  0.8  1.0 


(f)  i  =  2,  m  =  0;  for  sum 
frequency  generation 


=  AW.  ~  V.Ul 

0n904ym.  Th 
1.89  «*I0'«  cm 


0.80  0.85  0.90  0.95  1.00 


Figure  11.  (Concluded) 
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Finally,  the  coherent  and  incoherent  intensity  for  a  slab  of  tenuous  silver 
microsphere  composite  embedded  in  a  dispersive  host  wil  be  calculated.  The  configura¬ 
tion  is  shown  in  Figure  1.  The  dispersive  host  is  chosen  to  be  an  artificial  glass.  The 
first  incident  frequency  (i>|  is  chosen  to  be  the  surface  dipole  frequency  of  a  silver 
sphere  in  the  glass  host,  viz.,  uj  =  Ug  =  0.4304  ym  with  ef^(u)g)  =  2.15.  The  second  incident 
frequency  U2  is  either  tuned  to  the  vicinity  of  Wg,  viz.,  W2  =  Ug  -  0.05up  =  0.4608  ym  with 
e^(w2)  =  2.14.  Or  the  frequency  W2  may  be  the  difference  frequency  between  the  surface 
dipole  and  quadrupole  resonances,  viz.,  U2  =  -  Ug,  with  «  0.4  ym  and  e|^(u2)  =  1.9. 


(d)  A  =  2,  m  =  0;  u_  =  0.01  top. 


(e)  i  =  2,  m  =  0;  =  Uq. 


(f)  &  =  2y  m  s  0;  =  2u^  -  0.01  up. 


r/a 


r/a 


r/a 


Figure  12.  (Concluded) 


The  coherent  intensity  of  the  sum  (u^)  or  difference  (u_)  frequency 
generation  peaks  strongly  near  the  perfect  phase  matching,  viz.,  L  <  1,  where  A^  is 
defined  in  Equation  44  and  L  is  the  slab  thickness.  The  phase  mismatch  A^  defines  a 
coherent  length  =  l/A^.  Only  for  L  <  the  coherent  intensity  is  significant  and 
increases  more  or  less  quadraticly  with  L.  For  a  few  millimeters  thick  silver 
microsphere  composite,  the  sufficently  long  coherent  length  is  only  achievable 
under  the  perfect  phase  matching,  viz.,  A^  =  0. 

In  order  to  have  perfect  phase  matching,  the  host  dielectric  constant  at  the 
sum  (difference)  frequency  must  be  less  (greater)  than  C)^(uj)  and  e^(u2)«  The  incident 
angle  ^2  ^  determined  by  via  the  perfect  phase  matching  condition.  The  angle  62 
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will  then  set  the  coherent  angle  e+  through  the  stationary  condition  given  in 
Equation  40c.  It  is  no  loss  of  generality  to  restrict  the  values  of  02  to  the  range  between 
0®  and  90®,  since  the  range  90®  <  02  ^  180®  is  just  the  mirror  image  of  the  former  range 
with  respect  to  the  x-y  plane.  The  physical  values  of  02  are  going  to  limit  the  possible 
choices  of  the  host  materials,  viz.,  e^(u)^).  A  slab  of  silver  microsphere  composite  to  be 
considered  consists  of  radius  a  =  100  A,  thickness  L  =  1  mm  and  volume  fraction 
f  =  10'^.  Figure  14  shows  the  dependence  of  the  incident  and  coherent  angles  on  ej^((i)+). 
For  the  sum  frequency  generation  near  2a)g  or  at  the  possible  choices  of  e^(a>^)  are 
respectively  in  the  range  of  (1.08,  2.14)  or  (1.8,  2.13).  The  difference  frequency 
generation  at  <i)_  =  0.05  Up  requires  that  2.3  <  <  900. 

Figure  15  shows  the  dependence  of  the  coherent  and  incoherent  intensities 
upon  ej^(«+).  The  reduced  coherent  and  incoherent  intensities  are  1^(62/)  1^20 1 

*  *‘®sp€Ctively.  The  conversion  factor  to  the  corresponding  absolute 
intensity  is  given  in  the  figure  caption  for  each  case.  The  coherent  intensity  of  the  sum 
frequency  generation  near  2uig  is  about  1  W/cm  ,  which  is  10  dB  above  the  corresponding 
incoherent  intensity  for  1.2  <  e^(u^)  <2.1.  The  coherent  intensity  of  the  sum  frequency 
generation  at  Uq  is  0.05-0.2  W/cm^,  which  is  27  to  30  dB  above  the  incoherent  intensity 
for  1.75  5  ^  2.1.  The  coherent  intensity  at  the  difference  frequency  at  u)_  =  0.05  Up 

is  only  1-10  yW/cm^,  which  is  40  to  50  dB  above  the  incoherent  counterpart  for 
ej^(u)_)  >  2.4.  For  the  sum  frequency  generation,  the  incoherent  intensity  may  still  be 
detectable.  However,  for  the  difference  frequency  generation,  only  the  coherent 
intensity  can  be  observed  with  ease. 

In  Figures  16(a)  and  16(b),  the  size-dependence  of  the  dimensionless 
function  U^q  in  the  quadrupole  moments  for  dc  rectified  (U^q)  and  second-.  .;»rmonic 
(U20)  generation  for  the  case  of  a  single  intense  laser  irradiated  on  a  silver  microsphere 
composite  are  shown.  Both  figures  show  increasing  quadrupole  moments  for  decreasing 
particle  sizes  and  vanishing  quadrupole  moments  for  extremely  small  spheres.  These  are 
the  similar  size-dependent  characteristics  shown  in  Figure  5  for  the  linear  dielectric 
functions.  However,  the  quasi-hydrodynamic  theory  adopted  here  breaks  down  for 
extremely  small  particles,  where  the  quantum  effect  is  important  and  should  be  included 
with  care. 
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3.3  THIRD-ORDER  SUSCEPTIBILITY  IN  THE  FREQUENCY  DOMAIN 

For  the  third-order  susceptibility,  both  the  self-action  effect  of  the  laser 
fields  and  the  cross-action  effect  between  the  laser  fields  at  the  degenerate  and 
nondegenerate  frequencies  will  be  considered. 

The  self-action,  which  is  described  in  terms  of  the  intensity-dependent 
nonlinear  dielectric  function,  is  responsible  for  the  self-phase  modulation  and  self- 
focusing  of  a  given  laser  field.  They  are  important  effects  for  a  single  laser  field  as  well 


(a)  £  =  0,  m  s  0; 

u_  =  0.01  Wq  (unit 

V20  =  6.4  x^l03  cm/sec). 


(b)  i  =  0,  m  =  0; 

=  ui„  (unit 

V20  =3.5  X  103  cm/sec). 


(c)  1  =  0,  m  =  Oj 

=  2«i^  -  0.01  oIq 

(unit  V20  =  1*7  xTO**  cm/sec). 


Figure  13. 


Second-order  dimensionless  radial  velocity  versus  r/a, 

as  100  A.  Same  configuration  treated  in  Figure  2. 
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(d)  a  =  2,  m  =  0; 

io_  =  0.01  w_  (unit 

V20  =  6.4  x^l03  cm/sec). 


(e)  a  =  2,  m  =  0; 

u)^  =  u_  (unit 

V20  X  103  cm/sec). 


(f)  4  =  2,  m  =  0; 

=  2<i)e  -  0.01 

(unit  V20  =  1«7  x^io**  cm/sec). 


Figure  13.  (Concluded) 

as  for  phase  conjugation  and  coherent  two-wave  mixing.  The  cross-effect  is  described  by 
the  cross-intensity-dependent  nonlinear  dielectric  function  and  coherent  coupling 
coefficient  which  are  responsible  for  the  cross-phase  modulation  and  coherent  coupling 
between  the  laser  fields.  They  are  important  effects  in  the  phase  conjugation  and 
coherent  two-wave  mixing.  In  particular,  the  coherent  coupling  coefficient  is  the  main 
mechanism  in  the  coherent  two-wave  mixing. 

The  electrodynamic  problem  has  been  treated  globally,  calculating  the  dioole 
moment  induced  upon  the  sphere  by  the  applied  field,  and  thus  obtaining  directly  the 
surface-plasma  resonance  within  the  quasi-hydrodynamic  formalism.  Therefore,  if  the 
incident  frequency  is  near  the  surface  plasma  resonant  frequency  of  the  microparticles, 
there  will  be  local  field  enhancement  of  the  third-order  susceptibility  just  as  the  theory 
states. 
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(a)  u)^  =  2u)g  -  0.05  Up  =  0.2225  pm. 


(b)  Uq  =  0.39  -  0.41  pm. 


(c)  u  =  0.05  Up  =  6.54  pm. 
f  =  lO"**,  L  =  I  mm,  and  a  =  100  A. 


thto,) 


Figure  14.  The  incident  62  and  coherent  angles  versus  e^(u^). 

3.3.1  Third-Order  Suscqptibility  in  the  Frequency  Doniain 

3.3.1. 1  Self-Action  and  Cross  Action  Effects.  For  self-action,  consider  a  metallic 
sphere  composite,  which  is  irradiated  by  an  intense  laser  beam  of  frquency  uj,  amplituoe 
A|,  propagating  along  the  x-axis  and  polarized  in  the  z-direction.  The  equation,  which 
governs  the  induced  third-order  density  fluctuation  at  u|(r,uj)  is  Equation  12. 
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£h  (“+I 

(a)  u)^  =  2uig  -  0.05  Up.  The 
conversion  factors  are 
29.19  W/cm2  and  0.12  W/cm2 
for  coherent  and  incoherent 
intensities. 


(b)  u^  =  Up,  the  conversion 
factors  are  0.23  W/cm2  and 
0.30  mW/cm2  for  coherent 
and  incoherent  intensities. 


13.76  yW/cm2  and 
6.65  nW/cm2  for  the 
coherent  and  incoherent 
intensity,  f  =  lO"**, 

L  =  1  mm,  and  a  =  100  A. 


Figure  15.  The  reduced  coherent  and  incoherent  intensities  versus  e^(u+). 


The  induced  third-order  fields  are  given  as  follows.  Inside  the  sphere,  the 
density  and  electrostatic  potential  are  given  by 
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•C441U 


•no-®  cm) 

(b)  For  second-harmonic  part,  Q(2u)/a^. 

Figure  16.  The  dimensionless  function  Ujq  versus  particle  sizes. 


47 


Figure  16.  (Concluded) 


CXitside  the  sphere,  the  electrostatic  potential  is  given  by 


o(«J 

♦3  - cose  +  --J—  YjqCb)  (50) 

where  u^^^wj)  and  o(«»})  are  the  induced  third-order  dipole  and  octupole  nrtoments  of  the 
sphere.  The  induced  third-order  dipole  moment  is  given  in  Appendix  C  (Equation  C.IO). 
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The  third-order  octup>ole  moment  (t  =  3)  is  much  smaller  compared  to  the 
third-order  dipole  moment  (4  =1)  for  spheres  in  the  Rayleigh  limit.  Thus,  the  third-order 
octupole  moment  contribution  to  the  third-order  susceptibility  is  negligible. 

For  a  very  dilute  metallic  sphere  composite  slab,  the  third-order  coherent  field 
at  point  r  outside  the  slab  generated  by  microparticle  third-order  polarization  at  uj  is 
given  by 


=  fL2irik.{-0sineP<^^(a.j))  (51) 

where  f  is  the  volume,  L  is  the  thickness  of  the  slab,  and  (u»^ )  is  given  in 
Appendix  C  (Equation  C.ll). 

The  corresponding  third-order  susceptibility  of  the  composite  is  given  by 


f  -  ■  -I-  ■  ..2 


(52) 


where  the  dimensionless  function  Pj  |jn(uj)  is  given  in  Appendix  C  (Equation  C.IO),  and 
the  local  field  enhancement  factors  are  clearly  given. 

For  a  silver  sphere  of  100  A  in  radius  embedded  in  a  host  with  =  1,  one  has 


3.59  X  10"^° 


(53) 


After  lengthy  numerical  integration,  one  obtains  at  the  surface  plasma  frequency 


x(3>(«5) 


(-  5.3  +  1.81)  X  10”®  esu/sphere 


(54a) 
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(54b) 


for  V  =  0.02  Up,  and 


*  (-  8.3  +1  4.1)  X  10*^^  esu/sphere 


for  V  =  0.1  Ujj,  respectively.  The  experimental  value  (Ref.  18)  for  single  silver  sphere  is 
2.4  X  10  esu,  which  was  measured  in  the  degenerate  four-wave  mixing  configuration, 
and  is  fairly  close  to  the  value  given  in  Equation  54a. 

The  degenerate  third-order  susceptibility  for  a  very  dilute  silver  sphere 

(3) 

composite  with  volume  fraction  f  is  given  by  f  x'  '(uj,  where  the  numerical  values  of 

(3)  * 

x'  '(ug)  are  given  by  Equations  54. 

For  cross-action  effects,  consider  the  metallic  microsphere  composite,  which 
is  irradiated  by  two  intense  laser  beams  of  frequencies  u^  and  U2»  electric  field 
amplitudes  and  ^2*  both  propagating  along  the  x-axis  and  linearly  polarized  in  the 
z-direction.  Then  the  third-order  susceptibility  under  this  configuration  is  given  by 


(3) 


y  (»  1 


OjC.,)  |0j(.j)|  ,55, 


where  the  i  =  j  term  is  the  self-action  contribution  as  given  by  Equation  52  and  the  i  *  j 
term  is  the  cross-action  contribution. 

According  to  the  classification  discussed  in  Section  2.3,  the  third-order 
susceptibility  in  Equation  55  for  i  =  1  can  be  cast  into  the  following  groups. 


(3)  /  X  ^  (30/  .  X  ^  (3X)/  .  V 

*  (3S)/  .  X 


(56) 


which  are  the  coherent  coupling,  cross-phase-modulation  and  self-phase  modulation  parts 
of  the  susceptibility,  respectively. 
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In  Equation  56,  the  coherent  coupling  part  is  given  by 


"  4zzz^“1*“2  4zzz^“l*’“2'*’“+^ 


(57a) 


which  comes  from  the  cascading  sum  and  difference  frequency  effects. 

In  Equation  56,  the  cross-phase  modulation  part  is  given  by 

’‘zzzz^“l*“l*“2*'“2^*4zzz^“r*“l'*'  (“2"“2^^‘^4zL^“1*“1*“2*"“2^ 


(57b) 


which  comes  from  the  cascading  cross-rectified  effect  and  direct  third-order  contribu¬ 
tion. 

In  Equation  56,  the  self-phase  modulation  part  is  given  by 


4zzz^“l*“l*“l**“l^  '  ’‘zzzzW*"!'^ 

(3) 


which  comes  from  the  cascading  self-rectified  and  self-second-harmonic  effect  and  the 
direct  third-order  contribution. 


3.3.1.2  Polarization  Effect.  For  polarization  effect,  consider  the  more  general 
situation  discussed  in  Section  2.0  (Fig.  1).  There  are  an  elliptically  polarized  wj  wave 
with  polarization  vector  ej^  =  y+nz  and  a  linearly  polarized  «d2  wave  with  polarization 
vector  ©2  “  ^  incident  normal  to  the  microparticle  composite,  i.e.,  02  =  it/2.  The 
normal  direction  of  the  composite  or  the  propagation  direction  of  both  waves  is  in  the 
x-direction. 
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For  a  very  dilute  metallic  sphere  composite  slab,  the  third-order  coherent  field 
at  point  r  outside  the  slab  generated  by  microparticle  third-order  polarization  at  is 
given  by 


=  fL2ir1kj^{-eiis1neP^^^(a»j)  +  (ecosesin4>-4>cos*)Pp^(uj)} 

(58a) 

=  fL2ii1k2  {-esin0P^^)(a.2)}  (58b) 

(3) 

where  f  is  the  volume,  L  is  the  thickness  of  the  slab,  and  P||j  '  (u^)  is  given  in 
Appendix  C  (by  Equation  C.ll). 

Then  the  nonvanishing  components  of  the  third-order  susceptibility  are  given 
by 


1 _ 


(59a) 


for  self-action  contribution  and 


X  - i - 5 

(3)  (3) 


(59b) 


‘^l("i)  (59c) 


with  j  *  i  for  cross-action  contribution,  and  i  =  1  or  2. 

The  third-order  susceptibility  in  Equation  58  can  also  be  grouped  into  the 
coherent  coupling,  cross-phase-modulation  and  self-phase  modulation  parts  as  was  done 
for  Equation  55. 
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3.3.2  Time-Dependence  of  Third-Order  Susceptibility  and  Polarization 


In  Section  3.3.1,  the  third-order  susceptibility  of  metallic  microsphere 
composites  in  the  frequency  domain,  which  is  convenient  to  use  in  the  case  of 
monochromatic  or  quasi-monochromatic  light  beams  was  calculated.  However,  in  the 
case  of  ultrashort  light  pulses,  the  third-order  susceptibility  in  the  time  domain  will  be 
more  useful. 

The  more  general  time-dependent  third-order  nonlinear  polarization  is  given  by 
p(^>(t)  .  f  dtj  f  dtj  dt,  E(t-tj)  E(t-t2)  E(t-t3) 

0  0  0  (JO) 

According  to  Equation  60  the  time-dependent  third-order  susceptibility  for 
metallic  microsphere  composites  can  be  obtained  from  the  Fourier  transform  of  the 
frequency-dependent  third-order  susceptibility  in  Section  3.3. 1  and  is  given  by 


/4\  dul  dui  d<u4  wit  't'wito) 

X<3>(ti,t2.t3)  •  /  5(1  /  5^2  J  e  '  1  2  3  3' 


(61) 


where 

W  *  wj^  -t-  W2  <i>2* 

In  Section  2.3,  the  different  classification  of  the  third-order  responses  was 
discussed.  Here  only  the  time-dependent  third-order  polarization  for  coherent  coupling 
terms  is  worked  out  explicitly.  If  the  metallic  microsphere  composite  is  irradiated  by 
the  two  Gaussian  light  pulses  at  carrier  frequencies  wj  and  W2  with  pulse  full  widths  T^ 
and  T2,  and  probe  delay  time  tq,  the  time-dependent  third-order  susceptibility  for  .the 
metallic  microsphere  composite  can  be  obtained  from  Fourier  transformation  of 
Equations  39. 

For  a  linearly  polarized  probe  wave  at  W2»  the  coherent  coupling  parts  of  the 
time-dependent  third-order  susceptibility  are  given  by 
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where  6(12)  is  a  step  function,  the  demominator  is  given  by 


di(«s)=Re  {Di(ai^)} 


(63a) 


the  inverse  of  intrinsic  decay  time  is  given  by 

I«Dl(«s) 


(63b) 


and  t|  and  t2  are,  respectively,  the  time  coordinate  for  the  pump  and  probe  waves.  In 
Equation  61,  the  P-f unctions  are  given  by 


P'l^lii(®2)  “  ^  ®2^2^*^2®^  (64a) 

which  are  the  same  homogeneous  and  inhomogeneous  factors  shown  in  Appendix  C 
(Equations  C.IO). 

For  simplicity,  it  is  assumed  that  both  pump  (u^)  and  probe  (W2)  waves  are 
propagating  along  the  x-axis,  i.e.,  62  s  «/2.  The  envelopes  of  two  pulses  are  chosen  to  be 


•ST-  '*■ 

I4  y  2  p 

A^Cx.tp  =  2  [— J— --1  exp  I-  (t^-  /  4T‘l 


where  Ij  and  vj  are  the  light  intensity  and  phase  velocity  in  the  host  medium,  and  Tj  is 
the  pulse  full  width,  i  =  1  or  2. 

Then  the  total  electric  field  for  the  two-wave  mixing  can  be  written  as 

*  -iu-lt  ^  -iWj(t-Trj) 

ECr.t)  =  ejAj(x,t)  e  +  e2A2(x,t  -  Tp)  e  ^  +  c.c. 


where  tq  is  the  probe  delay  time. 

Substituting  Equations  62,  63  and  66  into  Equation  60,  the  coherent  parts  of 
the  time-dependent  third-order  polarization  for  the  probe  wave  at  u)2  are  given  by 
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where  the  inverse  of  effective  decay  time  for  short  pulses  with  relative  delay  is  given  by 


'’e  •  '’s  - 

'2 

the  overlap  function  is  given  by 


(68a) 


0(t,Tp)  -  exp[-(t-Tp-  /  2t|  -  (t-x/vj)2  /4^2j 


(68b) 


and  the  dimensionless  polarization  prefactors  for  the  homogeneous  and  inhomogeneous 
contributions  are  given  by 
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(68c) 


H(vJe)  =  2T^  [1  -  VgP(vJ^)l 


erfc(vgTg) 


(68d) 


with  the  effective  pulse  width  T_,  which  is  defined  by 


-2 

'e  '2  'l 


(68e) 


In  Equation  66,  the  subscript  m  =  1  or  0  corresponds  to  the  y-  or  .z-component  of  the 
time-dependent  third-order  polarization  of  the  probe  wave  at  U2» 


3.4  SPECTRAL  RANGE  AND  RESPONSE  TIME 

The  spectral  range  of  the  nonlinear  third-order  susceptibility  of  the  micro¬ 
particle  composites  is  determined  by  the  influential  range  of  the  local  field  enhancement 
factor.  This  can  be  easily  determined  by  the  FWHM  of  the  enhancement  factor.  For 
silver  microsphere  composites  from  Figures  7  and  8,  the  spectral  range  is  between 
0.40  »inn  and  0.44  ym.  For  gold  microsphere  composites,  there  is  a  smaller  enhancement 
factor  and  larger  spectral  range  because  of  gold's  larger  material  intrinsic  damping 
fitting  parameter  as  compared  to  silver's. 

The  response  times  of  the  microparticle  composite  can  be  read  out  from  their 
time-dependent  polarization  given  in  Equation  67.  The  response  time  is  the  effective 
decay  time  (l/v^)  in  the  equation,  which  is  determined  by  the  intrinsic  decay  time  (1/vg) 
of  the  microparticles  as  well  as  the  ratio  of  the  probe  delay  time  (tq)  to  its  pulse  width 
The  intrinsic  decay  times  for  the  parameters  used  in  Figures  3  through  10,  are  20  fs 
and  3  fs  for  v  =  0.02  Wp  and  v  =  0.1  Up,  respectively.  These  intrinsic  decay  times  are 
typical  of  the  electronic  nonlinearity. 
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SATURATION  CHARACTERISTICS 


In  a  real  material  the  optical  nonlinearity  cannot  grow  without  limit  and  many 
materials  will  permit  nonlinear  changes  only  in  a  refractive  index  that  saturates  quite 
rapidly  with  power.  Saturation  often  limits  refractive  index  changes  to  an  upper  limit  of 
0.1  but  apparently  0.001  is  fairly  common  with  a  lower  limit  being  typically  0.0001. 

In  this  theory  the  nonlinear  electrodynamics  of  microparticles  is  developed  in 
terms  of  perturbation  expansion.  The  saturation  will  set  in  if  the  expansion  parameter  is 
approaching  unity. 

From  the  calculated  perturbed  density,  electric  field  or  other  fields,  the 
expansion  parameter  can  be  easily  extracted.  Then,  by  equating  the  expansion  parameter 
to  unity,  the  saturation  field  and  intensity  for  the  nonlinear  electrodynamics  of  the 
microparticles  can  be  obtained. 

Comparison  of  the  perturbed  density  given  in  Appendices  A,  B  and  C  (by 
Equations  A.2,  B.7  and  C.6)  shows  that  the  following  order  of  magnitude  relationship 
exists  between  the  second-order  nonlinear  density  and  the  first-order  linear  density  and 
the  nonlinear  density  of  successive  orders 
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where 


r  =  ing-  DWuJ  = 
sat  ea 


3eqyp 


(70) 


with  A^  and  D|(u|)  are  the  electric  field  amplitude  and  the  resonant  denominator  of  the 
incident  light.  In  Equation  70,  r^  is  a  widely  used  measure  of  the  electron  density,  and 
defined  as  the  radius  of  a  sphere  whose  volume  is  equal  to  the  volume  per  free  (or 
conduction)  electron.  Note  that  the  saturation  electric  field  given  in  Equation  70  is 
reduced  by  the  surface  plasma  resonant  effect,  viz.  Dj((>ip,  as  well  as  the  background 
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dielectric  constant  ejj  and  is  enhanced  by  the  Thomas-Fermi  screening  effect,  viz.  the 
inverse  of  the  Thomas-Fermi  screening  length  qjp. 

The  saturation  fields  for  the  parameters  used  in  Figures  3  through  10,  are  5.1 
and  U9.6  k  statvolts  per  centimeter,  which  correspond  to  the  saturation  intensities,  being 
3.0  and  200  GW/cm  for  silver  spheres  with  v  =  0.02  Up  and  0.1  Up,  respectively. 
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4.0  EXPERIMENTAL  RESULTS 


4.1  PHASE  CONOUGATION  EXPERIMENTS 

4.1.1  Phase  Conjugation  in  CS2 

4.1.1. 1  Experimental  Apparatus.  Phase  conjugation  in  CS2  was  studied  experimentally 
utilizing  the  four  wave  mixing  experimental  system  shown  in  Figure  17.  A  passively 
mode  locked  Nd:YAG  oscillator  produced  a  train  of  S  to  10  pulses,  with  a  pulse  duration 
of  25  ps  in  a  TEMqq  transverse  mode.  A  single  pulse  was  selected  out  of  the  pulse  train 
and  amplified  through  two  single  pass  amplification  stages.  Between  the  two  amplifiers, 
the  beam  was  expanded  by  a  factor  of  3  by  means  of  an  inverted  telescope,  to  increase 
the  amplification  and  to  reduce  the  divergence  of  the  beam.  The  amplified  pulse  was 
frequency  doubled  to  532  nm,  resulting  in  a  single  picosecond  pulse  with  a  maximum 
energy  of  15  m3,  which  corresponds  to  a  peak  power  of  ~  600  MW.  Vertical  polarization 
of  the  frequency  doubled  pulses  was  ensured  by  a  polarizing  beam  splitter.  Beam  splitter 
BSj  (R  =  0.15)  was  used  to  create  the  probe  beam  while  beam  splitter  BS2  (R  =  0.50) 
reflects  a  small  portion  of  the  probe  pulse  into  a  reference  pulse,  Ij..  The  reference  pulse 
is  reflected  back  to  BS2  by  mirror  Mg,  whose  reflectivity  was  chosen  to  meet 
experimental  requirements.  The  probe  beam  traverses  the  sample  at  an  angle  of  6**  with 
the  pump  beam  I^.  The  temporal  overlap  was  adjusted  by  the  variable  delay  line 
consisting  of  mirrors  M3  through  M3.  The  backward  pump  beam  I|,  was  created  by 
retroref lection  of  I|  by  mirror  Mj  in  close  contact  with  the  sample  cell  to  ensure  a 
maximum  temporal  overlap  of  the  pulses  within  the  nonlinear  medium.  The  0.55-mm 
thickness  of  the  cell  wall  added  to  the  retrorefiection  setup  an  inherent  delay  of  3  ps 
between  the  two  pump  pulses.  This  delay  is  much  smaller  than  the  pulse  width  and  the 
experimentally  determined  coherence  length  of  the  pulses,  which  was  14  ps.  In  some 
experiments  a  special  cell  was  used,  the  rear  end  of  which  was  a  mirror.  This  cell  was, 
however,  more  difficult  to  handle  and  did  not  improve  the  maxiumum  attainable  phase- 
conjugate  reflectivity.  The  €$2  was  HPLC-grade  (glass  distilled  and  filtered)  and  was 
contained  in  quartz  spectroscopic  cells  with  1-,  2-  or  lO-mm  optical  path  lengths.  The 
diameter  of  the  pump  beams  in  the  sample  cell  was  3  mm.  Some  of  the  experiments 
were  performed  focusing  the  prdbe  beam  into  the  sample  cell.  The  phase  conjugated 
beam  I^  retraced  the  path  of  the  probe  and  was  detected  by  a  fast  photodiode  (ITT 
F4000).  The  intensity  was  attenuated  by  three  glass  slides  (G)  which  could  be  replaced  by 
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mirrors  when  necessary.  The  reference  pulse  is  detected  by  the  same  photo  diode  and 
precedes  the  conjugate  pulse  by  7  ns.  Both  pulses  were  measured  using  a  Tektronix 
500  MHz  storage  oscilloscope.  The  detection  system  was  calibrated  by  placing  mirror 
M[^  in  the  path  of  the  probe  beam.  Pulse  energies  were  measured  by  a  pyroelectric 
energymeter  (Laser  Precision  R37200). 


mOTODIODE 


Figure  17.  Experimental  setup  for  the  production  of  phase  conjugation  by 

degenerate  four-wave  mixing  of  picosecond  pulses.  The  backward  pump 
pulse  is  produced  by  retroref lection  of  the  forward  pump  beam.  M; 
mirror;  D:  diaphragm;  SHG:  second  harmonic  generator;  1^:  forward  pump 
beam;  IpZ  proi^  beam;  1^:  phase  conjugated  beam;  1^:  reference  beam. 


4.1. 1.2  Experimental  Results  and  Discussion.  Using  the  experimental  system 
described  in  the  previous  section,  phase  conjugate  reflectivity  was  measured  as  a 
function  of  pulse  energy  for  different  values  of  probe  to  pump  pulse  fluence  Fpj./Ff,  for 
different  ratios  of  the  two  pump  pulse  energies  Ef/E|,  and  for  different  cell  lengths  L. 
Also,  the  presence  of  stimulated  Brillouin  scattering  and  whole  beam  self-focusing  of  the 
pump  pulses  were  checked.  Both  of  these  nonlinear  effects  were  found  to  be  negligible. 
In  an  attempt  to  detect  two-photon  absorption  in  CS2f  the  transmission  of  picosecond 
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pulses  was  measured  as  a  function  of  their  energy.  However  no  significant  energy 
dependence  of  the  transmission  could  be  detected  up  to  an  energy  of  5  m3  per  pulse. 

In  the  retroreflection  setup,  the  maximum  interaction  length  of  the  two  pump 
beams  is  limited  to  one  half  of  the  pulse  length,  so  it  is  advantageous  to  minimize  wasted 
overlap. 

In  the  first  experiment  a  specially  designed  cell,  2-mm  long,  with  the  rear 
surface  consisting  of  a  dielectric  mirror  was  used.  Using  this  cell,  the  region  of  temporal 
overlapping  of  the  two  pump  beams  in  the  CS2  was  increased.  Figure  18  shows  a  typical 
result  for  reflectivity  versus  pump  pulse  energy,  using  the  special  designed  cell.  The 
ratio  of  pump  to  probe  pulse  energy  was  14.  The  decrease  in  the  slope  of  the  reflectivity 
versus  energy  was  accompanied  by  a  decrease  in  the  quality  of  the  phase  conjugated 
beam,  which  was  attributed  to  small-scale  self-focusing;  the  pump  pulse  seems  to  break 
up  in  bright  and  dark  regions,  and  this  same  distorted  profile  is  transferred  to  the  phase 
conjugate  beam.  The  power  density  at  which  saturation  takes  place  is  about  1  GW/cm  , 
although  theory  predicts  a  critical  power  threshold  for  macroscopic  self-focusing  in  CS2 
at  10  W.  The  short  optical  path  length  of  the  sample  cell  causes  the  self-focusing  to 
occur  outside  the  cell.  In  this  case  there  is  no  change  into  the  amplitude  profile  of  the 
beam  inside  the  cell,  only  a  change  in  the  phase  profile,  i.e.,  the  shape  of  the  wavefronts. 
According  to  Akhmanov  a  1-mm  cell  containing  a  Kerr  like  medium  acts  as  a  thin  lens 
with  focal  length  Zj,  given  by  Reference  16; 


4 


(71) 


where  Wq  is  the  1/e  beam  radius  at  the  cell  entrance,  Iq  is  the  intensity  at  the  beam 

center,  n2  the  nonlinear  refractive  index  and  L  the  cell  length.  At  a  power  level 
2 

Iq  =  1  GW/cm  ,  Wq  =  3  mm  and  L  =  1  mm,  Zj  equals  5  m  in  the  case  of  CS2.  This  shows 
that  the  short  interaction  region  in  picosecond  FWHM,  permitting  the  use  of  a  short 
sample  cell,  reduces  the  effects  of  macroscopic  or  whole  beam  self-focusing. 
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PUMP  PULSE  ENERGY,  mJ 

Figure  18.  Phase-conjugate  reflectivity  of  CS2  in  a  2-mm  cell  versus  the  pump 
energy  in  the  forward  pump  beam,  using  a  sample  cell  with 
retroref  lection  mirror  within  the  sample.  The  pump/probe  pulse  energy 
ratio  is  14:1.  The  curve  clearly  indicates  the  negative  deviation  from  a 
square  dependence  of  the  reflectivity  on  the  pump  energy  at  high 
energies  and  is  typical  for  most  experiments  we  performed. 


In  order  to  determine  the  influence  of  a  difference  in  energy  of  the  forward 
and  backward  pump  pulses,  the  conjugate  reflectivity  versus  pump  pulse  energy  was 
measured  using  a  retroref  lector  with  a  reflectivity  less  than  unity.  Figure  19  shows 
reflectivity  versus  reduced  energy  curves  for  three  different  values  of  the  backward  to 
forward  punip  pulse  energy  ratio  Taking  into  account  the  Fresnel  losses  at  the 

sample  cell  window,  the  following  values  of  E||,/E|  were  measured:  0.91,  0.71  and  0.30. 
Reduced  pump  energy  E^  is  defined  as  the  geometric  mean  of  the  forward  and  backward 
pump  pulse  energies  E^  and  E^:  Examination  of  the  plot  in  Figure  19  makes 

evident  the  sharp  decrease  in  efficiency  with  the  increasing  energy  difference  between 
forward  and  backward  pump  pulses.  The  solid  lines  represent  the  theoretical  curves 
obtained  using  the  model  described  in  Reference  16.  The  wave  vector  k  of  an 
electromagnetic  field  depends  on  the  intensity  in  a  nonlinear  medium  since  the  refractive 
index  is  dependent  on  the  electromagnetic  field  strength.  Therefore,  for  two 
counterpropagating  fields  with  identical  vacuum  wave  vectors,  kf  and  k^,  within  a 
nonlinear  medium,  kf  k|j  will  not  be  equal  to  zero,  resulting  in  a  phase  mismatch  of  the 
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Figure  19.  Phase  conjugate  reflectivity  of  CS2  in  a  2-mm  sample  cell  versus  the 
reduced  pump  energy  for  three  different  values  of  the  ratio  of 
backward  to  forward  pump  pulse  energy  E^/EfS  (•)  0.91,  (o)  0.71, 

(a)  0.30.  The  solid  lines  are  calculated  curves,  assuming  a  1/e  beam 
radius  of  1.3  mm  (see  text).  The  ratio  of  forward  pump  pulse  energy  to 
probe  pulse  energy  was  14:1. 


four  wave  mixing  process.  If,  however,  the  vacuum  intensities  are  equal,  then 
k|  +  k|j  =  0,  and  there  should  be  no  dephasing.  The  mismatch  can  be  calculated  as  a 
function  of  the  energy  ratio  r  and  the  forward  pump  pulse  intensity  I^,  analogous  to 
Reference  17: 


•  |k,+  '‘bl  •  J  “  ,  ^  If 


(72) 


Using  the  nonlinear  coupling  coefficient  k: 


(73) 
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and  being  the  intensities  of  forward  and  backward  pump  beam,  the  expression 
(Equation  72)  can  be  written  as  following: 


|Ak|  *  9<  [^]  K  (74) 

/r 

Using  Equation  74,  the  phase  conjugate  reflectivity,  R,  with  asymmetric  pump  beams  can 
be  calculated  as  follows: 


R  =  I  s1n^»L/lV)  _  ,  (^5) 

cos^(icL/l+p^)  +  p^ 

A  ^ 

where  p  =  -g-  (l-R//r)  and  L  is  the  interaction  length  of  the  pump  and  probe  pulses 
inside  the  nonlinear  medium.  The  theoretical  curves  in  Figure  19  were  calculated  using 
the  value  of  =  2  x  10"^^  (SI  units),  leading  to  the  following  expression  for  the 
coupling  constant  in  €$2: 

K  =  0.21  If  (MW/cm^m)  (76) 


The  curves  were  calculated  assuming  a  beam  diameter  of  3  mm.  For  the  ratio 
Eb/Ef  =  0.53  the  theory  predicts  the  maximum  reflectivity  very  well.  The  different 
maximum  positions  of  the  theoretical  and  experimental  curves  are  attributed  to  the  fact 
that  the  beam  diameter  was  an  estimate.  These  results  show  that  for  ratios  Ejj/Ef  of 
0.71  and  0.91  and  at  pump  energies  higher  than  0.75  m3,  the  experimental  data  for  the 
phase-conjugate  reflectivity  fall  below  the  theoretically  expected  values,  even  when  the 
dephasing  due  to  pump  pulse  energy  asymmetry  is  taken  into  account.  The  experimental 
data  shown  in  Figure  19  lead  to  the  same  conclusion,  because  in  this  experiment  the 
pump  pulse  difference  was  minimized.  TTie  conclusion  is  that  pump  pulse  asymmetry  is 
not  the  limiting  factor  for  the  phase  conjugate  efficiency  in  these  experiments. 
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A  second  important  factor  which  limits  the  efficiency  of  phase  conjugation  is 
the  depletion  of  the  pump  beams.  It  is  obvious  that  the  four -wave  mixing  process  can 
never  produce  a  conjugate  pulse  which  has  a  higher  energy  density  or  fluence  than  that  of 
the  pump  pulses.  However,  this  effect  has  only  seldom  been  taken  into  account  when 
solving  the  coupled-wave  equations  for  the  optical  phase  conjugation  (OPC)  by 
degenerate  four-wave  mixing  (DFWM)  (Refs.  18,19).  Numerical  solutions  to  these 
equations,  incorporating  pump  depletion,  indicate  no  significant  depletion  effects  when 
the  pump  to  probe  intensity  ratio  is  higher  than  10  (Ref.  19).  In  Figure  20  the  results  of 
six  different  experiments,  each  having  a  different  ratio  of  probe  to  pump  fluence  Fpj./Ff 
are  compiled.  The  maximum  reflectivity  obtained  in  each  experiment  is  plotted  as  a 
function  of  the  Fpj./Ff.  In  all  experiments,  except  the  one  with  Fpf/Ff  =  5,  this 
maximum  was  the  saturation  level  of  the  reflectivity  occuring  at  pump  pulse  energy 
levels  above  2  m3.  Note  that  the  power  dependence  in  each  of  these  experiments  was 
measured  by  varying  the  intensity  of  the  source  beam  so  as  to  obtain  reflectivity  data  at 
a  constant  pump  to  probe  fluence  ratio.  In  the  normal  setup,  Fp^/Fj  was  equal  to  0.05. 

In  order  to  investigate  the  influence  of  variations  in  Fp^/Fj,  this  ratio  was  decreased  by 
neutral  density  filters  in  the  probe  beam.  At  the  lowest  value  of  Fp^/Fj  a  reflectivity  of 
more  than  300  percent  was  obtained,  which  means  phase  conjugated  amplification  of  the 
probe  pulse.  The  reflectivity  versus  energy  curve  for  this  experiment  is  shown  in 
Figure  21.  Note  that  this  curve  shows  very  little  saturation  effect,  as  compared  to 
Figure  18  for  example.  The  cell  length  was  2  mm  in  this  case.  In  the  other  extreme, 
where  the  ratio  Fp|./Fj  was  increased  to  5  by  focusing  the  probe  beam  into  the  sample 
cell,  a  reflectivity  of  only  0.1  could  be  obtained.  The  probe  beam  was  focused  into  the 
center  of  the  collimated  pump  beams,  reducing  its  radius  by  a  factor  of  10,  causing  the 
fluence  of  the  probe  beam  to  increase  by  a  factor  100.  However,  focusing  of  the  probe 
beam  has  been  shown  to  increase  the  quality  of  the  phase  conjugation  (Refs.  20,21): 
under  certain  conditions  a  lens  transforms  the  probe  beam  to  its  spatial  Fourier 
transform  in  the  focal  plane,  located  in  the  interaction  region  of  the  pump  beams.  In  this 
case,  the  Gaussian  reflectivity  profile  determined  by  the  product  of  the  pump  pulse 
transverse  intensity  profiles,  acts  as  a  spatial  filter,  because  the  reflectivity  is  highest  in 
the  center  of  the  Fourier  plane,  where  the  low  spatial  frequency  components  of  the  probe 
beam  are  concentrated.  In  this  case  the  phase  conjugated  pump  beam  had  the  same 
diameter  as  the  probe  beam,  whereas  a  phase  conjugated  reflection  of  a  collimated  probe 
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Figure  20.  Maximum  value  of  the  phase  conjugate  reflectivity,  as  obtained  from 

different  experiments,  plotted  versus  the  ratio  of  probe  pulse  fluence  to 
forward  pump  beam  fluence,  Fp^/F|. 

beam  reduces  its  beam  radius.  The  results,  however,  indicate  that  focusing  of  the  probe 
beam  leads  to  a  stronger  dependence  of  the  PC-efficiency  on  the  energy  of  the  probe 
pulse,  reducing  the  fidelity  of  the  phase  conjugation. 

Self-focusing  of  the  pump  beams  is  very  important  as  a  limiting  factor  for  the 
fidelity.  Although  whole  beam  self-focusing  was  never  observed,  the  phase-conjugated 
beam  at  high  pump  pulse  energy  was  degraded  by  small  scale  self-focusing,  leading  to 
beam  breakup.  Experiments  using  a  10-mm  sample  cell  length,  indicate  a  lower 
maximum  reflectivity  than  2-mm  long  cells.  Normally  one  would  expect  a  slight  increase 
due  to  the  increased  overlap  region.  This  effect  is  attributed  to  an  increase  of  the 
wavefront  distortions  induced  by  small  scale  self-focusing  in  a  longer  sample  cell. 

In  this  discussion,  generally,  pump  energy  rather  than  pump  intensity  has  been 
used  as  the  independent  parameter.  It  was  not  possible  to  measure  accurately  the 
beam  diameters  used  in  these  experiments.  Variation  of  the  beam  diameters  is 
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Figure  21.  Phase  conjugate  reflectivity  of  CS2  in  a  2-mm  cell,  at  a  pump  to  probe 
energy  ratio  of  500:1,  plotted  versus  forward  pump  pulse  energy.  A 
maximum  reflectivity  of  3  has  been  observed,  i.e.,  an  amplification  of 
the  phase  conjugated  reflection  of  the  probe  beam. 


estimated  to  be  20  percent  from  experiment  to  experiment  since  realignment  was 
involved.  In  addition  to  this,  slight  variations  of  the  pulse  width  could  occur  due  to 
changing  concentration  of  the  passive  modelocker  dye. 


4.1. 1.3  Conclusion.  The  efficiency  of  the  phase  conjugate  reflection  of  picosecond 

pulses  in  CS2»  a  transparent  optical  Kerr  medium  has  been  measured.  Very  high 
reflection  efficiencies  were  obtained  although  some  saturation  of  the  reflectivity  at  high 
fluences  was  observed.  This  was  attributed  to  "small  scale"  self-focusing,  which  is  a 
typical  problem  in  high  energy  optical  systems.  The  increased  efficiency  for  the  phase 
conjugation  of  short  pulses  in  Kerr  media  is  explained  by  the  fact  that  the  influence  of 
the  processes  such  as  stimulated  Brillouin  scattering  (SBS)  and  whole  beam  self-focusing, 
which  limit  the  phase  conjugate  efficiency  of  nanosecond  pulses,  are  reduced  when 
picosecond  pulses  are  used.  In  addition,  a  strong  dependence  of  the  phase  conjugation 
efficiency  on  the  ratio  of  the  pump  to  probe  fluences  was  observed. 
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4.1.2  Phase-Conjugation  in  Gold  Colloicb 


4.1.2.1  Experimental  Setup.  Two  types  of  gold  colloid  samples  were  used,  a  water- 
based  colloid,  prepared  following  the  method  of  Turkevich  (Ref.  22)  and  a  solid  glass- 
based  sample  obtained  from  Schott  glass  (sample  type  RG-6,lmm  thickness).  A  water- 
based  colloid  solution  contains  approximately  74  mg  of  gold  per  liter  of  water.  Since  the 
density  of  gold  is  19.3  g/cm  ,  this  corresponds  to  a  total  volume  of  4  x  10  cm  or  a 
volume  fraction  (Volg^y/Vol^^^^l)  of  ~  4  x  10"^.  Typical  values  for  the  radius  of  the  gold 
spheres  in  such  solutions  range  from  -  40  A  to  ~  120  A,  corresponding  to  sphere  volumes 
between  2.7  x  10  and  7.2  x  10  cm  .  Electron  micrographs  indicated  that  particles  in 
the  solutions  had  an  average  radius  of  120  A.  This  is  equivalent  to  a  number  of  spheres  per 
liter  ranging  from  1.5  x  lO^^/l  to  5.6  x  10 ^^/l.  The  water-based  colloid  had  a  local 
absorption  maximum  due  to  surface  plasmons  on  the  gold  particles,  measured  by  a  Cary 
219  spectrophotometer,  of  2.76  cm~^  at  519  nm  for  a  sample  thickness  of  2  mm. 

Similarly,  the  absorption  of  the  Schott  RG-6  sample  had  a  local  maximum  of  9.21  cm~^ 
at  535  nm.  This  was  superimposed  on  a  background  absorption,  due  to  interband 
transitions,  which  becomes  noticeable  at  -  620  nm  and  increases  in  magnitude  as 
wavelength  decreases.  The  width  of  the  peak  for  the  local  maximum  was  approximately 
100  nm  for  both  samples  before  bleaching.  After  bleaching,  the  width  could  not  be 
estimated  accurately. 

Optical  phase  conjugation  measurements  were  made  in  two  different  pump 

intensity  regions.  The  phase  of  the  nonlinear  susceptibility  was  measured  by  means  of  a 

nonlinear  interferometer,  described  elsewhere  (Ref.  23).  A  semiconductor  colloid  glass, 
(3) 

whose  yc  '  (both  magnitude  and  phase)  was  previously  determined,  was  used  in  one  arm  as 
a  reference  and  a  gold  colloid  sample,  whose  phase  was  to  be  determined,  was  placed  in 
the  other  arm.  The  interferometer  operated  at  a  wavelength  of  532  nm  via  a  CW  mode- 
locked  laser  providing  1200  mW  of  average  power  at  a  repetition  rate  of  82  MHz.  The 
peak  intensities  at  the  sample  were  -  1  to2  kW/cm  . 

Since  the  phase-conjugate  reflectivities  of  most  high  speed  (i.e.  picosecond 
response)  materials  is  very  low  at  the  low  peak  intensity  levels  of  CW  mode  locked  lasers 

(a  few  kilowatts  per  square  centimeter),  scattering  is  comparable  to  the  P.C. 

-10  -9 

reflectivities  of  ~  10  -10  .  For  this  reason,  a  lock-in  amplifier  was  used  in  addition  to 

a  series  of  diaphragms  to  reject  scattered  light. 
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A  pulsed  Nd:YAG  laser  operating  at  2  Hz  and  332  nm  was  used  for  the  high 
intensity  phase  conjugate  (P.C.)  reflectivity  measurements  (see  Figure  22).  A  passively 
mode-locked  oscillator  was  followed  by  a  pulse  selector,  two  amplifier  stages,  and  a 
frequency  doubling  crystal.  This  provides  a  23  ps,  332  nm  pulse  with  a  diameter  of 
~  2  mm  and  energies  ranging  from  a  few  tens  of  microjoules  to  a  few  millijoules.  This 
pulse  is  split  into  four  components.  Two  counterpropagating  pump  beams  and  one  weaker 
probe  beam  are  arranged  in  the  counterpropagating  DFWM  geometry.  The  fourth 
component  is  used  as  a  reference  signal.  The  reference  signal  as  well  as  the  phase- 
conjugate  signal  are  detected  by  high-voltage  vacuum  photodiodes.  These  are  scanned  by 
gated  integrators  at  a  4  ns  gate  width  (long  enough  to  allow  for  trigger  jitter  and  short 
enough  to  minimize  detection  of  stray  light).  These  integrators,  in  turn,  are  monitored 
with  a  computer  interface  (A/D  conversion)  and  ultimately  the  data  are  manipulated  by  a 
DEC  microvax  computer. 


Figure  22.  Laser  system  for  high  intensity,  picosecond  DFWM  experiments. 
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The  water-based  colloid  suspension  was  studied  both  in  a  static,  closed  cell  and 
in  a  flowing  cell  in  circuit  with  a  200-cc  reservoir.  The  flowing  cell  arrangement  was 
used  to  minimize  the  influence  of  sample  deterioration  during  reflectivity  measurements, 
while  the  static  cell  was  used  for  absorption  band  measurements  before  and  after 
damage. 

During  some  experiments,  the  bleaching  of  the  gold  colloid  glass  was 
monitored  by  measuring  the  absorption  in  the  sample.  The  intensity  transmitted  through 
the  four  wave  mixing  region  in  the  sample  was  measured  by  a  photodiode,  as  well  as  a 
reference  signal.  Thus,  this  system  could  monitor  the  influence  of  each  successive  laser 
pulse  on  the  absorption  of  the  gold  colloid  sample.  The  absorption  monitor  at  ^43  nm  was 
assumed  to  be  indicative  of  absorption  changes  at  532  nm,  the  wavelength  at  which  the 
phase  conjugation  experiments  were  carried  out. 

4. 1.2.2  Results.  Data  from  the  high-energy,  pulsed  YAG  laser  for  gold  colloid  glass 

and  water  suspensions  provided  a  power  dependence  reference  for  the  phase  conjugate 
reflectivity  for  intensities  ranging  up  to  a  few  gigawatts  per  square  centimeter 
(Figures  23,  24a).  In  the  initial  studies,  the  probe  beam  power  was  ~  10  percent  of  the 
pumps  and  in  later  experiments  it  was  adjusted  to  be  50  percent  of  the  pump  power.  No 
qualitative  difference  (due  to  the  probe/pump  ratio  change)  was  detected  in  the  structure 
of  the  measured  reflectivity  power  dependence.  At  these  low  overall  reflectivities,  no 
influence  from  pump  depletion  was  expected  since  a  probe  intensity  comparable  to  a 
pump  still  would  drain  the  pump  by  only  0.001  percent. 

Data  for  the  water-based  suspension  indicate  that  the  reflectivity  does  not 
vary  significantly  with  pump  intensity  in  the  region  of  ~  500  MW/cm^.  This  does  not  mean 
that  the  P.C.  reflectivity  is  constant  but  rather  that  it  does  not  vary  strongly  and 
certainly  does  not  vary  quadratically  with  the  pump  intensity  as  would  be  expected  for  a 
P.C.  process. 

For  the  water  suspension,  some  experiments  indicated  a  "flat"  spot  or  local 
maximum  in  the  intensity  dependence  at  a  pump  energy  flux  of  a  few  hundred  micro¬ 
joules.  At  higher  pump  intensities,  the  reflectivity  begins  to  rise  again  as  a  function  of 
pump  intensity.  It  is  difficult  to  specify  this  power  dependence  analytically  since,  in  the 
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Figure  23.  Phase  conjugate  reflectivity  of  Schott  RG-6  gold  colloid  glass  for  high 
pump  intensity,  probe/pump  ratio  =  0.1,  low  repetition  rate,  25  ps  pulses 
for  all  beams  incident  (■);  pump  1  blocked  (O);  pump  2  blocked  (a)  and 
all  beams  incident  on  an  undoped  glass  sample  (+).  The  flat  response  for 
gold  glass  is  largely  due  to  the  formation  of  permanent  gratings  in  the 
sample.  All  beams  are  linearly  copolarized. 


pump  energy  range,  the  slope  of  the  curves  varies  continuously  on  a  logarithmic  scale  but 
the  range  is  not  large  enough  to  assign  an  exponential  power  dependence  (Ref.  19).  At 
first  glance,  these  power  dependence  curves  do  not  display  characteristics  that  would  be 
associated  with  the  usual  approximations  (i.e.  'Hangent-squared*  power  dependence  or 
saturation  of  reflectivity).  These  results  will  be  discussed  below. 

During  exposure  to  the  high  power  pulsed  laser,  the  Schott  RG-6  glass  gold 
colloidal  suspension  developed  permanent  gratings  which  caused  the  '^hase  conjugate” 
reflectivity  to  remain  approximately  constant  over  a  wide  pump  beam  energy  range.  One 
possible  (partial)  explanation  for  this  result  is  the  possibility  of  total  particle 
decomposition  in  the  regions  of  the  antinodes  of  the  optical  standing  waves.  Thus,  the 
total  diffraction  grating  would  be  the  sum  of  the  linear  gratings  formed  by  sample 
damage  and  the  nonlinear  gratings  which  still  exist  at  the  nodes  of  the  optical  standing 
waves.  If  the  linear  contribution  is  assumed  to  dominate  in  the  experiment,  constant 
reflectivity  would  be  expected.  Interestingly,  the  total  phase  conjugate  signal  is 
significantly  greater  than  the  sum  of  the  two  'Operate”  permanent  gratings,  i.e.  the 
signals  detected  when  one  or  the  other  of  the  pump  beams  was  blocked.  This  implies  the 
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P.C.  Reflectivity 


Pump  Energy  (ml) 

(a)  Phase  conjugate  signal  for 
water-based  gold  colloid  (+) 
and  plain  water  (■)  in  a  2-mm 
flow-through  cell. 


Pump  Energy  (mJ) 

(b)  Phase  conjugate  reflectivity  if 
the  signal  from  the  water  and  cell 
is  assumed  in  phase  (+)  and  antiphase 
<■)  with  that  of  the  gold  particles. 


Pump  Energy  (mJ) 

(c)  Phase  conjugate  reflectivity  when  water  and  cell  nonlinearites  are  assumed  to  be 
real  and  that  of  the  gold  particles  are  assumed  to  be  imaginary.  All  beams  are 
linearly  copolarized. 

Figure  24.  Phase  conjugation  signal  for  water-based  gold  colloid  and  plain  water  in  a 
2-mm  flow-through  cell. 
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dominance  of  the  nonlinear  contribution,  in  direct  contradiction  to  the  initial  premise.  A 
hand-waving  solution  is  to  say  that  the  nonlinear  solution  dominates  with  saturation  of 
the  nonlinearity  being  the  crucial  factor  in  the  case  of  all  beams  incident  on  the 
sample. 

There  was  also  a  permanent  change  in  the  absorption  of  the  colloids,  both 
liquid  and  solid  (Figures  25,  26),  upon  exposure  to  the  high-power  pump  pulses.  This 
appeared  as  a  "bleaching"  of  the  colloid  or  occasionally  the  appearance  of  a  pale  bluish 
spot.  These  changes  appeared  to  be  permanent.  Data  from  the  spectrophotometer 
confirmed  this;  absorption  dropped  in  the  blue  and  increased  in  the  red  wavelength 
region. 
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Figure  25.  Change  in  optical  absorption  at  x  =  543  nm  of  Schott  RG-6  as  a  function 
of  total  energy  flux  from  pulsed  Nd:YAG  laser  (x  =  532  nm)  for  average 
pulse  energies  of  0.1  m3/pulse  (-•■)  and  1.0  m3/pulse  (□). 


The  change  in  absorption  in  a  Schott  //RG-6  filter,  as  a  function  of 
accumulated  laser  energy  incident  upon  the  sample  spot  was  measured  for  two  different 
average  laser  pulse  energies,  0.1  m3  and  1.0  m3  (Fig.  25).  Note  that,  upon  interrogation 
by  1  m3  pump  pulses,  the  absorption  at  543  nm  has  permanently  dropped  from  0.41  to  0.3 
after  only  a  few  millijoules  of  accumulated  energy.  For  the  same  total  accumulated 
energy  flux,  the  lower  energy  pulses  have  done  far  less  damage.  The  most  extreme 
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Wavelength  (nm) 

Figure  26.  Continuous  absorption  of  water-based  gold  colloid  in  a  2-mm  cell  before 
(a)  and  after  (b)  sample  damage. 

case,  namely  irradiation  with  a  few  hundred  millijoules  per  second  from  the  CW  mode- 
locked  Nd:YAG  laser,  did  not  result  in  any  change  in  absorbance  although  the 
accumulated  energy  was  many  orders  of  magnitude  larger  than  the  pulsed  laser  provided. 

At  the  highest  pulse  energies,  there  was  also  an  increase  in  the  far  field  spot 
size  of  the  phase  conjugate;  the  transmitted  probe  pulse  showed  a  similar  phenomenon. 
Since  the  transmitted  probe  pulse  shares  this  "bloomed"  appearance  (generally  appearing 
as  a  somewhat  distorted  series  of  concentric  rings)  in  the  far  field,  the  most  probable 
cause  of  the  phase-conjugate  "bloom"  structure  is  that  it  is  an  accurate  replication  of  the 
self-focusing  or  self-defocusing  of  the  probe  signal.  That  this  phenomenon  occurs  at  high 
intensities  is  expected;  the  self-focusing  or  self-defocusing  phenomenon  is  also  due  to  the 
third  order  nonlinear  susceptibility. 

Figure  27  shows  phase  conjugate  signal  interference  and  pump  signal  inter¬ 
ference  from  Schott  RG~6  gold  colloid  glass  and  Schott  OG-530  semiconductor  colloid 
glass  as  a  function  of  the  delay  of  one  arm  of  the  interferometer.  Since  the  maxima  and 
minima  of  the  two  signals  coincide  and  the  phase  of  x'  of  OG-530  is  known  to  be 
negative  and  real,  it  was  determined  that  for  the  high  repetition  rate,  low  peak  intensity 
(2  kW/cm^)  laser  pulses,  the  phase  of  the  total  third  order  nonlinear  susceptibility  for  the 


Figure  27.  A  comparison  of  the  interference  of  conjugate  signals  (•)  and  reference 
signals  from  Schott  RG-6  gold  colloid  glass  and  Schott  OG-530 
semiconductor  colloid  glass  as  a  function  of  delay  in  the  arm  of  the 
nonlinear  interferometer.  Since  the  phase  of  the  P.C.  signal  from  CXi-530 
is  known  to  be  negative  and  real  in  this  region,  the  phase  of  RG-6  must 
be  the  same. 

(3) 

gold  colloid  sample  is  negative  and  real.  In  other  words,  for  x  -Xreai+ixjmaginary* 
was  determined  that  the  effective  Xi^aginary  negligible  compared  to  Xreal 
Xreal  **  negative.  Since  the  colloid  sample  is  actually  a  composite  system,  this 
cannot  be  attributed  to  the  particles  without  further  assuming  that  the  nonlinear 
susceptibility  of  the  glass  matrix  is  insignificant  compared  with  that  of  the  gold  particles 
at  this  power  level.  Evidence  supporting  this  claim  is  that  a  plain  glass  sample  of  1-mm 
thickness  did  not  provide  a  detectable  phase-conjugate  signal  when  pumped  by  the  CW 
mode-locked  Iciser  of  the  Twyman-Green  interferometer.  Therefore,  the  negative  and 
real  nonlinear  susceptibility  is  attributed  to  the  gold  particles. 

4.1.2.3  Conclusions.  The  intensity  dependence  of  the  picosecond  phase-conjugate 
reflectivity  in  gold  colloids  for  pump  intensities  on  the  order  of  1  GW/cm  has  been 
determined  experimentally;,  the  intensity  dependence  of  sample  damage  and  the  phase  of 
the  nonlinear  susceptibility  of  gold  colloids  has  been  investigated  via  interferometric 
methods  using  a  CW-mode-locked  Nd:YAG  laser;  and  also  the  significance  of  the  phase 
relationship  of  the  nonlinear  susceptibilities  in  a  multi-component  system  has  been 
recognized.  This  is  significant  in  an  applications  environment  since  the  functional 
intensity  dependence  could  have  its  shape  ^tailcM’-made”  by  adjusting  the  stoichiometry  of 
the  constituents. 
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4.2 


COHERENT  TWO-WAVE  MIXING  EXPERIMENTS 


4.2.1  Experimental  Setup 

A  pulsed  Nd;YAG  laser  operating  at  4  Hz  was  used  for  the  two-wave  mixing 
experiments.  The  laser  system  consisted  of  a  passively  mode-locked  oscillator  which 
provided  a  pulse  train  of  ~  10  35-ps  pulses  at  x  =  1.064  pm.  A  pulse  selector  directed  one 
of  the  most  intense  pulses  of  the  train  to  a  double-pass  NdrYAG  amplifier.  The  selected 
pulse  was  chosen  from  the  leading  half  of  the  pulse  train  to  minimize  the  influence  of 
pulse  chirp  on  the  two-wave  mixing.  This  pulse  was  subsequently  frequency  doubled  in  a 
KDP  (KH2P0^)  crystal  and  residual  1.064  pm  radiation  was  then  eliminated  by  a  harmonic 
beam  splitter.  Pulse  energies  at  532  nm  were  on  the  order  of  1  to  5  m3,  depending  on  the 
experiment,  with  a  standard  deviation  of  25  f>ercent  within  any  one  experiment.  In  all 
cases,  the  beam  diameter  was  -  2  mm  and  was  approximately  Gaussian. 

The  532  nm  pulse  is  split  into  two  components,  corresponding  to  the  pump  and 
probe  pulses.  Another  beam  splitter  deflected  roughly  8  percent  of  the  pump  pulse  into  a 
detector  for  normalization  purposes.  The  final  input  intensity  (and  energy)  ratio  was 
variable  with  1  <  ^purnp/'^probe  ^  depending  on  the  beam  splitter  selected.  A  delay 
stage  in  the  pump  beam  path  controlled  the  relative  arrival  times  of  the  pump  and  probe 
pulses  in  the  sample.  The  angle  between  the  pump  and  probe  wavevectors,  outside  the 
sample,  varied  depending  on  the  experiment,  from  <  0.3®  (~  4  mrad)  to  -  4.5® 

(~  78  mrad).  Since  these  pulses  provided  intensities  on  the  order  of  10^  W/cm^,  focusing 
the  beams  was  unnecessary.  The  initial  pump  reference  and  final  probe  energies  were 
measured  as  a  function  of  their  relative  arrival  times  in  the  sample  by  pyroelectric 
detectors.  The  analog  output  of  these  meters  were  sent  to  an  A/D  converter  which  was 
interfaced  with  a  computer.  The  apparatus  is  shown  schematically  in  Figure  28. 

In  some  experiments  both  the  pump  and  probe  were  linearly  polarized,  in 
others  the  pump  was  circularly  polarized  and  the  probe  was  linearly  polarized.  For  the 
experiments  in  which  the  pump  pulse  had  circular  polarization,  the  probe  pulse  had 
polarizing  beam  splitters  before  and  after  the  sample.  One  of  the  motivations  ior  an 
interaction  of  a  circularly  polarized  pump  and  a  linearly  polarized  probe  was  to  minimize 
the  influence  of  the  "Kerr  Shutter"  effect.  It  is  well  known  that  a  strong,  linearly 
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Figure  28.  Experimental  apparatus  and  beam  geometry  for  degenerate  two-wave 
mixing.  PL  =  polarizer;  DET  =  detector;  DS  =  delay  stage. 

polarized  pump  causes  a  transient  linear  birefringence  in  Kerr  media;  this  effect  is  often 
used  to  rotate  the  polarization  of  a  probe  pulse  and  thus  "time"  two  picosecond  pulses  via 
the  "Kerr  shutter."  A  circularly  polarized  pump  causes  a  transient  circular  birefringence 
in  a  material,  but  this  effect  is  expected  to  be  insignificant  in  comparison  to  the  linear 
effect.  For  this  reason,  it  is  expected  that  any  net  change  in  the  measured  probe  energy 
is  a  result  of  energy  transferred  from  the  pump  rather  than  polarization  rotation  of  the 
probe  energy. 

In  a  third  group  of  experiments,  both  pump  and  probe  pulses  had  circular 
polarizations;  a  comparison  was  made  between  an  interaction  of  pulses  with  identical 
polarization  helicity  and  an  interaction  of  pulses  with  opposite  polarization  helicity. 

The  samples  included  carbon  disulfide,  nitrobenzene,  semiconductor  colloid 
glass  and  a  water-based  gold  particle  colloid.  The  liquid  samples  were  contained  in 
quartz  cells  of  various  thicknesses.  All  of  these  materials  should  have  nonresonant 
electronic  polarizabilities  with  response  times  on  the  order  of  10~^^  s;  they  should  have 
additional  response  times  corresponding  to  molecular  reorientation,  and  various  decays 
from  other  excited  states.  Carbon  disulfide  has  a  characteristic  response  time  of  ~  1  ps 
(Ref.  28)  and  nitrobenzene  has  two  characteristic  response  times,  ~  8  ps  and  a  few 
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nanoseconds  (Refs.  24,23).  Likewise  the  semiconductor  glasses  have  two:  ~  30  ps  and  a 
few  nanoseconds  (Ref.  26).  Finally,  the  water-based  gold  colloid  should  have  a  very  fast 
response  corresponding  to  a  surface  plasmon  decay  (Refs.  27,28,29)  as  well  as  a  much 
slower  response  due  to  thermal  effects. 

4.2.2  Experimental  Results 

The  initial  experiments  were  designed  to  produce  visible  self  diffraction;  this 
phenomena  is  fairly  well  understood  (Ref.  30)  and  the  experiment  is  facilitated  by  the 
fact  that  the  propagation  directions  of  the  diffracted  orders  allows  detection  of  the 
process  by  eye.  At  small  angles,  several  orders  of  self  diffraction  are  visible  when  using 
CS2  as  the  active  medium.  A  sample  thickness  of  ~  5  mm  showed  the  best  compromise 
between  the  intensity  of  each  diffraction  order  and  the  number  of  diffraction  orders. 
Sample  thicknesses  of  1,2,  and  10  mm  also  provided  visible  self  diffraction,  but  the  1- 
and  2-mm  cells  did  not  provide  as  much  diffracted  beam  intensity  as  the  5-mm  sample 
while  the  10-mm  sample  gave  more  intense  diffracted  beams  but  fewer  visible  orders  of 
diffraction.  All  visible  self  diffraction  disappeared  for  angles  larger  than  -  2* 

(35  mrad).  Self  diffraction  was  present  both  in  the  case  of  a  circularly  polarized  pump 
and  a  linearly  polarized  pump  but  the  intensity  was  visibly  higher  in  the  case  of  parallel 
linearly  polarized  pump  and  probe  pulses.  Semiconductor  colloid  glasses  and 
nitrobenzene  also  displayed  self  diffraction  effects,  though  not  to  extent  of  the  CS2 
sample. 

Once  the  influence  of  self  diffraction  had  been  verified,  the  energy  of  the 
transmitted  probe  pulse  relative  to  the  pump  reference  signal  (before  the  sample)  was 
monitored.  In  these  first  experiments,  an  increase  in  the  probe  energy  of  as  much  as 
~  10  percent  was  measured,  depending  on  the  relative  positions  of  the  pump  and  probe 
pulses.  At  some  positions,  there  was  an  apparent  slight  loss  in  the  probe  energy,  but  this 
loss  was  of  the  same  magnitude  as  the  noise  of  instrumentation  so  its  presence  was 
neglected  at  that  time.  The  angle  between  the  pump  and  probe  pulses  was  then  increased 
to  -  4.5*  where  the  higher  diffraction  orders  were  weak  enough  to  be  considered 
insignificant.  It  was  found  that  the  probe  gains  and  losses  were  present  even  without  the 
influence  of  self  diffraction* 

By  further  experimenution,  the  losses  and  gains  of  the  probe  were  successfully 
quantified  as  functions  of  several  parameters.  In  all  of  the  experiments,  the  laser 
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intensities  were  typically  on  the  order  of  1  GW/cm  .  The  time  axes  shown  on  Figure  29 
are  not  absolute.  Although  the  scales  are  accurate,  the  origin  has  been  arbitrarily 
assigned.  The  actual  origin  (location  of  zero  delay  between  the  pump  cuid  probe)  is 
located  in  the  region  of  pump  probe  overlap.  For  an  experiment  in  which  the  pump  and 
probe  are  initially  identical,  it  may  be  easily  deduced  that  the  origin  occurs  where  the 
net  energy  transfer  is  zero  (at  the  crossing  point).  For  other  configurations,  the  origin  is 
shifted  slightly  with  respect  to  the  crossing  point. 

When  the  semiconductor  colloid  sample  was  tested  (Figure  29),  the  results 
were  those  of  a  typical  pump-probe  transient  absorption  experiment:  when  the  output 
probe  polarization  parallel  to  the  input  probe  polarization  was  monitored,  there  was 
increased  transmittance  (or  decreased  absorption)  when  the  pump  pulse  passed  through 
the  sample  first.  The  decay  of  the  transmittance  corresponds  to  the  excited  state  free 
carrier  decay  in  the  semiconductor.  There  was  also  a  coherence  spike  (Ref.  34) 
superimposed  on  the  transient  absorption  background.  When  the  output  probe 
polarization  perpendicular  to  the  input  probe  polarization  was  monitored,  there  was  only 
a  Coherence  spike"  present.  The  width  of  this  spike  was  -  15  ps. 

The  gold  colloid  also  displayed  (Figure  30)  a  coherence  spike  in  the 
perpendicular  polarization,  but  the  excessive  noise  in  the  parallel  polarization 
experiment  precluded  the  detection  of  any  coherent  phenomena  here.  Interestingly,  the 
gold  colloid  showed  an  increased  probe  transmittance  for  all  times  when  the  pump  passed 
through  the  sample  first.  It  is  believed  that  particle  damage  is  an  important  process  at 
these  high  intensities  and  can  explain  the  apparently  permanent  change  in 
transmittance.  The  nonlinear  coupling  in  the  gold  colloid  proved  to  be  the  weakest 
coupling  in  the  materials  studied. 

The  most  interesting  results  (Figures  31,  32)  appeared  for  the  Kerr  media,  CS2 
and  nitrobenzene.  Both  of  these  materials  showed  a  'kroherence  spike"  representing 
energy  transfer  from  the  circular  pump  pulse  to  the  probe  pulse  in  the  polarization 
perpendicular  to  the  original  probe  polarization.  The  width  of  the  spike  in  the  CS2  sample 
was  apqproximately  40  ps  and  the  peak  energy  transfer  corresponded  to  -  30  percent  of 
the  probe  input  energy.  The  corresponding  peak  energy  transfer  in  the  nitrobenzene  was 
>  15  percent.  The  width  of  the  spike  in  nitrobenzene  could  not  be  easily  defined  as  it 
apparently  had  a  narrow  spike  (of  width  >  15  ps)  superimposed  on  a  broader  (30-40  ps) 
coherent  background.  Similar  phenomena  (i.e.,  probe  diffraction  experiments)  have  been 
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a)  the  polarization  perpendicular 
to  that  of  the  initial  probe 
polarization. 


0^ 


b)  the  polarization  parallel  to 
the  initial  polarization. 


Relative  Probe  Delay  (psec) 


Figure  29.  Change  in  probe  transmission  in  a  semiconductor  doped  glass  sample  as  a 
function  of  relative  probe  delay.  Synchronization  of  the  pump  and  probe 
occurs  at  ~  200  ps  on  the  graphs  (i.e.y  axis  origin  is  at  bitrary).  Each  point 
represents  one  laser  shot. 

discussed  by  several  researchers  (Refs.  30,31,32).  The  shape  (intensity)  of  the  diffracted 
signal  as  a  function  of  probe  delay  is  explained  by  considering  the  pulse  durations 
(causing  the  broad  background)  and  the  coherence  time  of  each  of  the  pulses  (causing  the 
narrower  spike  centered  on  the  broad  background)  in  materials  with  transient  grating 
decay  times  which  are  large  compared  to  the  ct^rence  time  and  pulse  length,  so  that  an 
"integrated-intensity"  grating  may  be  assumed. 

The  13  to  20  ps  coherence  spikes  in  both  the  nitrobenzene  and  semiconductor- 
doped  glass  samples  suggest,  therefore,  that  the  laser  pulses  are  not  transform-limited. 
This  is  expected  since  passively  mode-locked  lasers  typically  have  incomplete  mode¬ 
locking. 
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a)  the  polarization  perpendicular 
to  that  of  the  initial  probe 
polarization. 
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b)  the  polarization  parallel  to 
the  initial  polarization. 
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Figure  30.  Change  in  probe  transmission  in  a  liquid  gold  colloid  sample  as  a  function 
of  relative  probe  delay  for  monitoring:  Synchronization  of  the  pump  and 
prcriie  occurs  at  ~  200  ps  on  the  graphs  (i.e.>  axis  origin  is  arbitrary).  Each 
point  represents  20  laser  shots. 


The  intensity  of  the  coherence  spike  was  measured  as  a  function  of  pump 
intensity  in  CS2<  For  normalization  purposes,  the  energy  of  the  probe  coherence  spike 
was  divided  by  the  incident  laser  energy.  Figure  33  shows  the  expected  square  law 
dependence  on  the  incident  beam  intensity. 

In  the  case  of  the  output  probe  polarization  being  parallel  to  the  input 
polarization,  the  coherent  signal  from  CS2  was  antisymmetric  with  respect  to  pump  and 
probe  relative  delay  times.  When  the  pump  pulse  arrived  in  the  sample  slightly  before 
the  probe  pulse,  there  was  a  net  energy  transfer  from  the  pump  into  the  probe.  This 
energy  transfer  was  on  the  order  of  20  percent  of  the  initial  probe  energy.  When  the 
probe  arrived  in  the  sample  slightly  before  the  pump  pulse,  it  lost  >  20  percent  of  its 
initial  energy.  This  result  is  in  agreement  with  the  set  of  integro-differential  equations 
presented  earlier.  Examination  of  the  factors  (A  -B  )  in  the  coupled  field  amplitude 


a)  no  polarization  discrimination. 


b)  polarization  parallel  to  the 
initial  probe  polarization. 


c)  polarization  perpendicular  to 
the  initial  polarization. 
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Figure  31.  Change  in  probe  transmission  in  a  CS2  sample  as  a  function  of  relative 
probe  delay.  Synchronization  of  the  pump  and  probe  occurs  at  ~  200  ps 
on  the  graphs  (i.e.,  axis  origin  is  arbitrary).  Each  point  represents  20 
laser  shots. 


equations  indicates  that  a  reversal  of  leading  and  trailing  pulses  is  expected  to  result  in  a 
reversal  of  the  direction  of  energy  transfer.  Experimental  results  were  inconclusive  for 
nitrobenzene;  although  the  nonlinear  optical  susceptibility  of  nitrobenzene  has  a 
component  with  a  relaxation  time  comparable  to  the  pulse  width,  the  excitation  of  the 
nonlinearity  is  not  instantaneous  and  therefore  the  effective  nonlinear  excitation  by 
these  laser  pulses  is  diminished.  The  magnitude  of  the  transfer  was  apparently 
<10  percent  of  the  probe  energy.  However,  the  noise  levels  of  the  instrumentation 
prevented  quantification. 
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a)  the  polarization  parallel 
to  that  of  the  initial  probe 
polarization. 


b)  the  polarization  perpendicular 
to  the  initial  polarization. 
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Figure  32.  Change  in  probe  transmission  in  a  nitrobenzene  sample  as  a  function  of 
relative  probe  delay.  Synchronization  of  the  pump  and  probe  occurs  at 
~  200  ps  on  the  graphs  (i.e.,  axis  origin  is  arbitrary).  Each  point 
represents  a  laser  shot. 


Figure  33.  Change  in  pr<^  transmission  in  the  polarization  orthogonal  to  that  of 
the  imtial  probe  polarization  as  a  function  of  pump  energy  in  the  region 
of  >  1  m3  per  pulse. 
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Experiments  were  also  carried  out  in  CS2  as  a  function  of  different  initial 
energy  ratios  between  the  pump  and  probe  pulses  for  two  polarization  combinations: 
linear-pump/linear-probe  and  circular-pump/linear-probe. 

In  the  experiments  with  pump  and  probe  beams  having  parallel  linear 
polarizations  (Figure  34),  the  results  were  generally  not  strictly  antisymmetric  as  a 
function  of  pump-probe  relative  delay.  As  the  initial  pump/probe  intensity  ratio 
increased,  the  probe's  maximum  gains  further  exceeded  its  maximum  losses.  For 
example,  imagine  a  case  in  which  the  pump  intensity  is  initially  three  or  four  times  that 
of  the  probe.  If  the  probe  preceeded  the  pump,  it  might  lose  2  percent  of  its  energy  to 
the  pump  as  they  passed  through  the  sample.  If  instead  the  pump  preceeded  the  probe, 
the  probe  might  experience  gains  of  6  to  7  percent.  Even  in  the  case  of  initially 
equivalent  pump  and  probe  intensities,  the  antisymmetric  coherent  signal  is  present.  In 
this  case,  the  probe's  (or  pump's)  maximum  gains  should  be  the  same  as  its  maximum 
losses.  However,  the  situation  is  complicated  experimentally  by  the  simultaneous 
occurence  of  self  diffraction.  The  energy  exchange  increases  with  increasing  interaction 
length,  while  the  finite  beam  diameter  requires  small  intersection  angles  for  a  long 
overlap  length.  Unfortunately,  self-diffraction  increases  with  decreasing  intersection 
angles.  The  result  is  that  energy  is  drained  from  both  the  pump  and  the  probe  samples  as 
they  pass  through  the  sample.  This  means  that  the  maximum  net  gain  for  the  probe  (and 
pump)  is  decreased,  while  the  maximum  net  loss  is  increased  due  to  the  presence  of  the 
self-diffraction. 

In  order  to  verify  this,  measurements  were  made  (Figure  33)  as  a  function  of 
the  intersection  angle  for  an  initial  pump/probe  intensity  ratio  of  unity.  At  small  angles, 
the  probes  losses  were  indeed  greater  than  its  gains;  as  the  intersection  angle  increased, 
the  probe's  gain  indeed  became  antisymmetric  as  a  function  of  pump/probe  delay. 

Analogous  experiments  were  performed  using  a  circularly  polarized  pump  pulse 
with  a  linearly  polarized  probe.  In  this  case,  higher  initial  pump/probe  intensity  ratios 
resulted  in  lower  maximum  probe  gains  and  higher  maximum  probe  losses.  Again,  as  the 
initial  pump/probe  intensity  ratios  approached  unity,  the  probe's  maximum  losses  became 
comparable  to  its  gains.  The  influence  of  the  seif  diffraction  was  similar  to  that  in  the 
linear  polarization  experiments. 
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a)  Initial  pump  probe  intensity 
ratio  of  3/1. 


b)  Initial  pump  probe  intensity 
ratio  of  3/1. 


c)  Initial  pump  probe  intensity 
ratio  of  2.3/1. 


d)  Initial  pump  probe  intensity 
ratio  of  2/1. 
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Figure  34.  Change  in  probe  transmission  in  a  CS2  sample  as  a  function  of  relative 
probe  delay  for  monitoring  the  polarization  parallel  to  the  initial 
polarization  for  four  initial  pump  probe  intensity  ratios  for  an 
intersection  angle  of  ~  0.3*.  Synchronization  of  the  pump  and  probe 
occurs  at  ~  130  ps  on  the  graphs  (i.e.»  axis  origin  is  arbitrary).  Each  point 
reixresents  20  la^r  shots. 
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Transmission  Transmission  Transmission 


a)  Intersection  angle  of  0.24*. 


Relative  Probe  Delay  (psec) 


c)  Intersection  angle  of  0.68*. 


d)  Intersection  angle  of  0.9*. 
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Figure  33.  Change  in  probe  transmission  intensity  in  a  CSj  sample  as  a  function  of 
relative  probe  delay  for  monitoring  the  polarization  parallel  to  the  initial 
polarization  for  five  intersection  angles.  Synchronization  of  the  pump 
and  probe  occurs  at  ~  130  ps  on  the  graphs  (i.e.,  axis  origin  is  arbitrary). 
Each  point  represents  20  laser  shots. 


A  final  experiment  in  CS2  demonstrated  that  the  nonlinear  coupling  of  two 
pulses  with  identical  circular  polarization  was  significantly  weaker  than  the  coupling  of 
pulses  with  opposite  circular  polarization  (Fig.  36).  This  is  explained  by  noting  that  the 
sum  of  two  corotating  circular  polarizations  is  itself  circular;  a  circular  polarization 
allows  only  weak  molecular  alignment  (Ref.  33).  However,  the  sum  of  two  counter¬ 
rotating  circular  polarizations  has  a  linear  component  which  can  induce  significant 
molecular  reorientation. 

The  results  of  most  of  these  experimental  configurations  discussed  are 
summarized  in  Table  1. 
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b)  two  identical  pulses  with  corotating  circular  polarizations. 

Figure  36.  Change  in  probe  transmission  in  a  CS2  sample  as  a  function  of  relative 
probe  delay.  One  can  deduce  that  zero  delay  (temporal  synchronization 
of  the  pump  and  probe)  occurs  a+  ~  170  ps  on  the  time  axis. 

Table  1 


Adjustable  parameters  of  degenerate  two-wave  mixing 


^^^olarizations 

Circular  Pump  /  Linear  Probe 

Circular  Pump  /  Circular  Probe 

Parameter^^^ 

Linear  and 
Parallel  Pump 
and  Probe 

Madtoi  Parallel 
to  incident  Probe 

Monitor 

Perpendicular  to 
incident  Probe 

Co-Rotating 

polarizations 

Counter-Rotating 

polarizations 

Initial  ratio  of 
Pump/Probe 
Intensity  Increases: 

Probe  gains 
exceed  probe 
losses 

Probe  losses 
greater  than 
probe  gains 

— 

Intersection  angle, 
a,  of  pump  and 
probe  decreases'. 

Probe  losses  typically  exceed 
probe  gains  due  to  tlw  influence 
of  self  diffiraction  losses. 

— 

Negligible 

Influence 

losses  to 
diifeacted 
beams  increase 

Increasing  laser 
intensities 

Increase  in  coupling  and 
energy  transfer  between  pump 
and  probe  pulses 

Probe  intensity 
increases  as  the 
square  of  punq> 

increase  in  coupling  and  energy 
transfer.  Note  that  corouting  coupling 
is  «  counterrotating  coupling. 

Increasing  optical 
path  length  in 
nonlinear  medium 

Energy  transfer  in  Caibon  Disulflde  showed  a  small 
incr^  with  path  length  as  optical  path  was 
increased  fimn  1mm  to  2mm  to  5mm  to  10mm 

— 

Variation  of 
Pump-Probe 
Relative  Delay: 

As  the  pulses  pass  through  the 
sample,  probe  loses  energy  if  it 
leads  pump;  gains  energy  if  it 
lags  pump 

I’robegainis 
ipprox.  symmetric 
with  delay.  Loss 
lot  possible 

Probe  gains  energy  when  it 
lags  pump;  loses  energy  when 
it  leads  pump. 

Correlation  btwn. 
pulse  thiraiion  and 
material  relaxation 
time 

The  nonlinear  coupling  which  allows  energy  transfer  is  strongest  when  the  relaxation  time  of 
the  optical  nonlineuity  is  comparable  to  the  pulses  length.  For  cases  in  which  the  nonlinear 
relaxinion  time  is  significantly  smaller  or  larger  than  the  pulse  width,  energy  transfer  is 
negligible.  Reqxmse  time  for  excitation  of  the  optical  nonlinearity  must  also  be  considered 
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3.0  COMPARISON 


5.1  LINEAR  ABSORPTION  SPECTRUM 

The  absorption  of  water-based  gold  colloids  in  a  2  mm  cell  was  given  by 
Figure  26.  Two  fitting  parameters,  ejj  =  11.6  and  v  =  0.09  Up,  were  used  to  interpret  the 
linear  absorption  spectrum  for  gold  particles.  The  two  parameters  appear  in  the  resonant 
denominator,  which  is  given  in  Appendix  A  (Equation  A.2a),  and  calculated  from  the 
measured  absorption  coefficient  at  the  surface  plasma  frequency  both  for  water-based 
gold  colloids  and  Schott  PG-6  glass  in  Section  4.1.2.  The  theoretical  absorption  spectra 
for  these  two  parameters  are  given  in  Figure  37.  There  is  indeed  a  reasonable  fit  for  the 
low  energy  wing  of  the  absorption  spectra,  where  the  interband  transition  is  not  the 
dominant  contribution. 


Figure  37.  The  linear  absorption  spectra  of  gold  colloids  within  the  quasi¬ 
hydrodynamic  theory.  The  fitting  parameters  are  =  1 1.6  and 
V  =  0.09  Up. 


5.2  PHASE  CONJUGATION 

Starting  from  one  of  the  simplest  equations  for  phase-conjugate  reflectivity 
(R): 

R  ~  tan2(KL)  (77) 

(3) 

where  K  is  defined  as  the  nonlinear  coupling  coefficient  (K  is  proportional  to  x'  '  times 
the  field  intensity)  and  L  is  the  interaction  length,  consider  the  low-conversion  (low 
reflectivity)  limit.  In  this  region,  the  phase  conjugation  reflectivity  is  proportional  to 
the  square  of  K  and  the  approximations  for  CS2  and  gold  colloid  samples  can  be  equated, 
correcting  for  the  reflectivity  difference  between  the  two  and  the  actual  volume 
fraction  (~  4  x  10"^)  for  the  gold  colloid: 

^Au  Colloid  =  *^CS2  * 
or 

<*^Au^ff)^  =  <‘^CS2^^^  *  ^eff  =  <78b) 

After  the  substitutions: 

Kau  =  (l-S  »  ‘0-*)-‘/2KcSj  (79) 

(3)  4 

or  equivalent,  that  the  x  the  gold  particle  is  -  10  times  that  of  CS2»  which  is  about 
-8 

10  esu.  The  phase  conjugation  reflectivity  of  the  matrix  is  insignificant  at  low 
intensities,  but  becomes  comparable  to  that  of  the  particles  at  higher  intensities  because 
of  a  saturation  of  the  nonlinear  susceptibilities  of  the  particles.  Thus  the  suspension  may 
be  "tailored"  to  vary  the  total  optical  nonlinearity  in  the  sample. 

Using  the  parameters  of  CW  absorption  and  picosecond  optical  phase  conjuga¬ 
tion  measurements  the  linear  and  nonlinear  properties  of  gold  colloidal  suspensions  were 

13) 

studied.  The  theoretical  values  of  x  for  the  gold  spheres  can  be  obtained  from 
Equation  49  (page  46)  using  the  fitting  paramatars  for  CW  absorption  in  Saction  5.i»  and 
is  on  the  order  of  10~^^  esu.  This  is  much  too  low  compared  to  the  experimental  value. 
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This  indicates  that  the  free  electron  contribution  which  is  accounted  for  in  the  model  is 
not  the  main  nonlinear  mechanism  in  the  gold  microparticle  composites. 

5.3  COHERENT  TWO-WAVE  MIXING 

5.3.1  Transient  Coherent  Two- Wave  Mixing  in  Kerr  Medium 

Recently  there  has  been  considerable  interest  in  coherent  beam  combination 
via  two-wave  mixing  in  nonlinear  media  due  to  its  potential  applications  to  optical  signal 
processing  and  high  energy  laser  physics  (Ref.  34).  This  process  proceeds  as  follows:  two 
laser  beams,  a  pump  and  a  probe  wave,  enter  a  nonlinear  medium  and  induce  an  optical 
index  grating.  If  the  phase  of  this  grating  lags  that  of  the  interference  field  pattern 
formed  by  the  two  incident  laser  beams  I  and  II,  energy  will  be  coherently  diffracted 
from  beam  I  to  beam  II.  The  diffracted  beam  I  acquires  the  phase,  frequency  and 
propagation  vector  of  beam  II. 

To  date,  most  research  on  coherent  beam  combination  via  two-wave  mixing 

has  focused  on  photorefractive  crystals  (Ref.  35)  using  degenerate  CW  beams,  where  a 

phase  lag  arises  due  to  spatial  nonlocality.  Many  of  the  details  of  this  process  can  be 

influenced  by  varying  the  orientation  of  the  crystal  relative  to  the  laser  beams.  Another 

scheme  for  creating  a  phase  lag  between  the  laser  interference  pattern  and  the  light 

induced  index  grating  is  to  utilize  nondegenerate  beams  in  optical  Kerr  media  or  other 
(3) 

x'  '  materials  with  finite  respK)nse  time  (Ref.  36).  For  this  case,  energy  is  coherently 
transferred  from  the  high-frequency  beam  to  the  low  frequency  one  and  the  gain  is  a 
maximum  for  Ot  =  1.  Here  Q  is  the  difference  in  laser  frequencies  and  t  is  the  medium 
response  time.  Finally,  nonlinear  instantaneous  Kerr  film  coated  on  a  metallic  substrate 
such  as  gold  or  silver  can  transfer  energy  between  surface  polar  itons  and  laser  beams 
(Ref.  37). 

The  role  of  beam  polarization  on  transient,  degenerate  beam  combination  via 
two-wave  mixing  in  CS2  utilizing  Gaussian  laser  pulses  whose  pulsewidth,  Xp,  is  not  too 
much  greater  than  x  was  examined.  For  this  situation,  the  finite  response  time  of  the 
medium  and  the  width  of  the  laser  pulses  are  responsible  for  the  appearance  of  the  re¬ 
quired  phase  lag.  Since  this  is  a  transient  phenomenon  dependent  upon  temporal  non¬ 
locality,  the  energy  transferred  from  one  beam  to  the  other  will  intimately  reflect  the 
time-dependent  characteristics  of  the  system.  Prior  research  has  focused  on  beams  that 
are  linearly  polarized  in  the  same  direction  (Ref.  38). 


90 


In  the  two-wave  mixing,  beam  combination  experiments  a  pulsed  Nd:Yag  laser 
operating  at  4  Hz  was  used.  The  laser  system  consisted  of  a  passively  mode-locked  oscil¬ 
lator  that  provided  ~  20-ps  pulses  at  532  nm  with  ~  1  m3.  In  all  cases,  the  beam  diameter 
was  2  mm  and  the  transverse  profile  was  essentially  Gaussian.  The  pulse  is  split  into  two 
components,  corresponding  to  the  pump  and  probe  waves  with  the  final  intensity  ratio 
equal  to  10.  A  delay  stage  in  the  pump  beam  path  controlled  the  relative  arrival  times 
between  the  pump  and  probe  pulses  in  the  sample.  The  angle  between  the  pump  and  probe 
wavevectors  was  4.5®,  and  the  optical  path  length  was  5  mm.  The  polarization 
dependence  of  the  coupling  between  the  two  beams  was  examined  using  a  configuration 
consisting  of  a  right-circularly  polarized  pump  pulse  and  a  linearly  polarized  probe  pulse. 
The  pump  reference  and  polarization  discriminated  transmitted  probe  energies  were 
measured  by  pyroelectric  energy  meters  as  a  function  of  their  relative  arrival  times. 

The  Kerr  medium  was  CS2  which  has  a  characteristic  response  time  of  1  to  2  ps. 

Experimentally,  it  was  observed  that  the  coherent  signal  due  to  energy 
transfer  was  asymmetric  with  respect  to  probe  delay  time  for  the  case  where  the  output 
polarization  was  parallel  to  the  input  probe  polarization.  When  the  pump  pulse  arrived  in 
the  sample  slightly  after  the  probe  pulse,  i.e.,  delay  time  <  0,  there  was  a  net  energy 
transfer  from  the  probe  into  the  pump.  This  energy  transfer  was  on  the  order  of 
20  percent  of  the  initial  probe  energy.  When  the  probe  arrived  in  the  sample  slightly 
after  the  pump  pulse,  i.e.,  tq  >  0,  it  gained  a  small  fraction  of  its  initial  energy.  Energy 
transfer  to  the  probe  beam  in  the  orthogonal  polarization  direction  was  also  measured. 

In  this  case,  the  energy  transfer  curve  was  symmetrical  with  respect  to  the  pulse  delay 
time  and  achieved  a  maximum  value  of  about  25  percent. 

Next,  a  theory  is  presented  for  coherent  beam  combination  via  transient  two- 
wave  mixing  in  optical  Kerr  media  with  degenerate  laser  pulses.  The  transient  dynamics 
of  most  Kerr  media  are  reasonably  well  described  by  the  Debye  relaxation  model.  The 
laser  beams  are  governed  by  the  wave  equation 

?  n2  .2  .  .2 

^  A  lE(r,t)  =  4  H  Pu,  (r.t)  (80) 
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where  ,t)  is  the  nonlinear  polarization  and  n  is  the  linear  refractive  index.  Since 
the  dynamics  of  the  nonlinear  polarization  are  governed  by  Debye  relaxation,  the  i— 
component  of  the  nonlinear  polarization  is  given  by 

=  E.(r.t)  /  df  exp[  -  -^^^1  E^(r.f  )E.(r.f )  (81) 

where  is  the  ajki—  component  of  the  third-order  optical  susceptibility  tensor. 

First  the  parallel  polarization  case,  in  which  the  pump  is  circularly  polarized 
and  the  probe  wave  is  linearly  polarized  in  the  x-direction  is  considered.  They  are  nearly 
colinear  and  prop>agate  in  the  z-direction.  The  energy  of  the  transmitted  probe  polarized 
in  the  x-direction  is  measured.  If  the  pump  is  not  significantly  depleted,  polarization 
changes  in  both  beams  may  be  neglected  so  that  the  electrodynamics  of  beam 
combination  involves  independent  quantities:  the  pump  (probe)  intensity  I|  (I2)  and  the 
phase  difference  ♦  =  -  ^2*  phase  difference  between  the  two  pulses  is  driven  by 

their  nonlinear  interaction  with  one  another  through  their  coupling  to  the  Kerr  medium. 
Within  the  slowly-varying  envelope  approximation 


L  Ii(z,t)  =  - 

(82a) 

L  l2(2,t)  =  +  2^/2(^^^/nc)d  [Iil2]^^^FXz,t) 

(82b) 

L  ♦(z,t)  =  2^/2^iia./nc)d[l/Ii  -  I/I2]  (Iil2)^^^F"  <2»t) 

(82c) 

where  L  =  (a/3z  +  (n/c)a/3t)  and  d  =  is  real.  Equations  82a  through  82c  describe 
energy  transfer  processes  involving  the  x-component  of  the  probe  wave.  The  coherent- 
coupling  grating  term  F(z,t)  =  FXz,t)  +  iF'Xz,t)  is  defined  by 

f(2.t)  -  7  /  <lt'|Ii(2.t')l2(2.t')|*/2e«p(l(»(2.t')  -  ♦(2.t)) 

—  (83) 
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which  gives  rise  to  coherent  energy  transfer  between  the  beams.  The  quantity,  41(7, t')- 
^(z,t)}  +  i  (t-t*)/T  is  driving  the  phase  difference  between  the  laser  interference  pattern 
and  the  index  grating  created  within  the  nonlinear  medium  by  the  laser  beams 
themselves.  For  steady-state  situations,  ^  is  independent  of  time,  and  energy  transfer 
will  proceed  only  if  the  laser  beams  are  nondegenerate.  However,  for  transient 
situations  ^  is  time-dependent,  and  a  transient  phase  lag  will  develop  between  the  index 
grating  and  the  laser  interference  pattern  due  to  the  finite  response  time  t  of  the 
medium.  Note  that  the  depends  upon  (I^  -  I2),  so  that  locally  the  direction  of  energy 
transfer  will  depend  upon  the  relative  beam  intensities. 

The  laser  beams  are  modeled  as  Gaussian  pulses,  with  a  delay  time  tq  between 
them  and  identical  pulse  widths  tp  =  IOt.  Numerically  solving  Equations  82  and  83  yields 
the  transmitted  pulse  energy  of  the  pump  and  probe  waves.  Figure  38  compares  theory 
to  experiment  for  the  fractional  change  in  energy  of  the  probe  versus  pulse  delay  time. 
The  solid  line  is  theory,  the  diamonds  are  experimental  data  points  with  error  bars  in¬ 
cluded,  where  laser  fluctuations  were  the  dominant  contribution  to  experimental  error. 
An  examination  of  this  figure  reveals  that  theory  and  experiment  are  in  excellent 
agreement,  with  theory  predicting  the  positions  of  the  main  features  correctly;  e.g.,  the 
zero  crossing  and  the  maximum  fractional  energy  transfer  and  are  in  good  overall 
quantitative  agreement.  The  greatest  discrepancy  occurs  when  the  energy  transfer  is 
small  between  the  probe  and  the  pump. 

The  physical  content  of  these  results  can  be  stated  as  follows.  If  the  pump 
wave  enters  the  medium  before  the  probe  and  the  time  delay  between  them  is  long  com¬ 
pared  to  the  medium  response  time,  their  overlap  is  insignificant  and  the  induced  index 
grating  is  negligible.  As  the  pulse  delay  decreases,  the  beam  overlap  grows  and  the  index 
grating  they  generate  becomes  substantial.  If  the  pulse  width  is  not  too  long  compared 
to  the  medium  response  time,  the  phase  difference  does  not  approach  steady-state  and 
the  phase  lag  between  the  index  grating  and  the  laser  interference  pattern  is  signifi¬ 
cant.  This  manifests  itself  in  the  form  of  increased  energy  transfer.  This  trend  con¬ 
tinues,  until  the  two  beam  phases  are  the  same  over  most  of  the  overlap  region. 

Although  the  index  grating  becomes  deeper,  the  phase  difference  is  only  weakly  driven  by 
the  nonlinear  coupling  and  the  phase  lag  begins  to  decrease.  This  behavior  is  in 
accordance  with  Equations  82a,  b,  and  c  and  manifests  itself  as  a  reduction  in  energy 
transferred  from  one  beam  to  the  other.  Eventually,  the  pulse  delay  reaches  a  point  in 
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Figure  38.  Percentage  of  energy  transfer  versus  pulse  delay  in  CSj  for  a 

pump/probe  intensity  ratio  of  10,  an  optical  path  length  of  5  mm  and 
parallel  polarization.  The  initial  pump  intensity  was  1  GW/cm*, 
d  =  8  X  10“‘*  cm2/GW  and  t-  =  20  ps.  The  diamonds  are  experimental 
data  points  with  error  bars  ^d  the  curve  is  theory. 


which  no  energy  is  transferred.  Once  the  probe  beam  leads  the  pump,  the  situation 
reverses  itself,  the  plutse  lag  changes  sign  and  energy  flows  from  the  probe  to  the  pump 
wave.  This  trend  continues  for  increasing  time  delays  until  a  point  is  reached  in  which 
the  beam  overlap  is  once  again  too  weak  to  induce  a  meaningful  index  grating  and  energy 
is  no  longer  diffracted  from  the  probe  wave  to  the  pump. 

Next,  the  cross-polarization  case  with  a  right-circularly  polarized  pump  wave, 
a  linearly  polarized  probe  beam  initially  in  the  x-direction  and  energy  transfer  in  the 
y-direction  of  the  probe  is  measured.  Again  changes  in  pump  polarization  are 
neglected.  The  dominant  grating  arises  from  the  coupling  between  the  pump  beam  and 
the  x-component  of  the  probe  wave  which,  through  medium  nonlinear  anisotropy,  twists 
the  polarization  of  the  diffracted  beam  into  the  y-direction.  Denoting  the  y(x) 
component  of  the  probe  intensity  by  l2y  (l2x) 

L  Ii(2,t)  s  -2^^^*»/nc)  d/6[Iil2x]*^^F^2,t)  (84a) 
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(84b) 


L  l2y(2,t)  =  +23/2(ii<o/nc)  (d/6Xlil2x]^^^P'<z,t)} 

L  4.(2, t)  =  ^  F"(z,t)/2i2x] 

where  off-diagonal  components  of  d^jj^j^  are  related  to  the  diagonal  ones  and  F  is  given  by 
Equation  83.  Note  that  the  phase  difference  does  not  depend  on  the  relative  intensities 
of  the  pump  and  probe  waves  so  that  energy  is  always  transferred  to  the  y-component  of 
the  probe  beam.  The  numerical  results  are  shown  in  Figure  39.  For  this  case,  the  phase 
lag  will  not  reverse  sign  and  energy  transfer  is  symmetric  with  respect  to  tq  =  0.  As  the 
delay  time  between  the  probe  and  pump  increases,  the  index  grating  decreases  with  a 
corresponding  reduction  in  energy  transfer. 


-20  -10  0  10  20 
Time  Delay  /  t 


Figure  39.  Percentage  of  energy  transfer  versus  pulse  delay  in  CSj  for  a 

pump/probe  intensity  ratio  of  10,  an  optical  path  length  of  5  mm  and 
cross  polarization.  The  initial  pump  intensity  was  I  GW /cm  2, 
d  =  8  X  10"*®  cm  2/GW  and  Xp  =  20  ps.  The  diamonds  are  experimental 
data  points  with  error  bars  OTd  the  curve  is  theory. 


An  examination  of  these  figures  reveals  that  theory  and  experiment  are  in 
excellent  agreement  for  all  but  the  shortest  pulse  delay  times.  For  the  smallest  pulse 
delay  times,  the  index  grating  involving  the  y-component  of  the  probe  beam  is  significant 
and  since  Equations  84a,  b,  and  c  ignore  this  effect,  they  underestimates  the  energy 
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transferred.  Note  that  although  the  fractional  energy  transfer  approaches  25  percent, 
pump  depletion  is  still  negligible  and  therefore  the  pump  polarization  remains  essentially 
circular. 

The  temporal  pulse  shapes  at  output  are  Gaussian  for  Gaussian  input  pulses 
with  pulse  widths  being  10  t.  Frequency  chip  is  negligible  for  the  optical  path  lengths 
used. 

To  summarize:  polarization  resolved  coherent  two-wave  mixing  experiments 
have  been  performed  using  degenerate,  frequency-doubled  Nd:Yag  picosecond  laser 
pulses  in  CS2,  with  a  circularly  polarized  pump  and  a  linearly  polarized  probe  wave. 
Theoretical  predictions  are  in  excellent  agreement  with  the  measurements,  and  the 
important  role  of  beam  polarization  in  this  process  can  be  seen  clearly  by  comparing  the 
parallel-polarized  and  cross-polarized  cases. 

5.3.2  Transient  Coherent  Two-Wave  Mixing  in  Metallic  Colloids  (Gold  and  Silver) 

As  a  comparison  to  the  Kerr  medium  discussed  in  Section  5.3.1,  the  coherent 
part  of  the  third-order  polarization  for  probe  wave  at  W2  in  a  Kerr  medium  is  worked  out 
with  the  Debye  relaxation  time  characteristics  and  is  given  by 

1/2 

O(t.Tp)  P  (Tg/t  ) 

(85) 

where  the  overlap  function  CXt,!^)  and  the  dimensionless  polarization  prefactor  P(Tg/T) 
are  given  by  Equations  68b  and  68d.  In  Equation  85,  both  input  pulses  are  taken  to  be 
Gaussian  and  are  defined  as  in  Equation  64:  t  is  the  phenomenological  relaxation  time  and 
tq  is  the  probe  delay  time.  Comparing  Equations  85  and  67,  the  inhomogeneous  part  of 
the  time-dependent  third-order  polarization  for  microparticle  composites  is  the  same  as 
a  Kerr  medium  with  the  Debye  relaxation  characteristics.  The  homogeneous  part  has 
additional  time  asymmetric  factor  resulting  from  the  additional  fourth  field  enhance¬ 
ment  factor  involved. 

In  Figures  40  through  42,  the  dimensionless  time-dependent  polarization 
prefactors  are  plotted  as  functions  of  the  probe  delay  times.  One  can  see  the  effective 
decay  times  are  much  shorter  in  the  heavily  damped  spheres,  viz.  v  =  0.1  Up.  When  the 
decay  time  is  much  shorter  than  the  pulse  widths,  the  CW  limit  for  shorter  pulses  takes 


p(3C),.x  _  dyki  -^“2^  ^ 
T  ®  CT, 


/2n  1. 
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(a)  The  homogeneous  part  versus  probe  delay  time,  with  v  =  0.02  Up. 


^0  ^  ^0 

(b)  the  inhomogeneous  part  versus  probe  delay  time,  with  v  =  0.02  Up. 

Figure  40.  The  time*dependence  prefactors  versus  probe  delay  time  for  silver 

colloids  in  water  (eu  =  1.7689)  and  glass  (e^  =  2.2^)  with  radius  of  50  A  or 
100  A,  pulse  widths  10  ps  and  t  s  x/v. 
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(a)  The  homogeneous  part  versus  probe  delay  time,  with  v  =  0.02  Wp. 


(b)  the  inhomogeneous  part  versus  probe  delay  time,  with  v  =  0.C2  utp. 

Figure  41.  The  time-dependence  prefactors  versus  probe  delay  time  for  silver 

colloids  in  water  (cw  s  1.7689)  and  glass  (e|^  s  2.25)  with  radius  of  30  A  or 
100  A,  pulse  widths  1  ps  and  t  s  x/v. 
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(c)  the  homogeneous  part  versus  probe  delay  time,  with  v  =  0.1  Up. 


”^0  ^  "^0 

(d)  the  inhomogeneous  part  versus  probe  delay  time,  with  v  =  0.1  Up. 

Figure  41.  (Concluded) 
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(a)  The  homogeneous  part  versus  probe  delay  time,  with  v  =  0.02  Up. 


(b)  the  inhomogeneous  part  versus  probe  delay  time,  with  v  =  0.02  Up. 

Figure  42.  The  time-dependence  prefactors  versus  probe  delay  time  for  silver 

colloids  in  water  (cu  s  1.7689)  and  glass  (e^  =  2.23)  with  radius  of  30  A  or 
100  A,  pulse  widths  0.4  ps  and  t  s  x/v. 
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effect  and  the  time-dependent  behaviors  are  solely  determined  by  the  overlap  of  the 
incident  pulses.  Note  that  there  is  more  asymmetry  in  the  homogeneous  part  than  the 
inhomogeneous  part. 

The  transient  behavior  for  the  metallic  microparticle  composites  for  the 
inhomogeneous  (driven)  third-order  polarization  is  the  same  as  the  Kerr  medium 
considered  in  Section  5.3.1}  since  they  are  driven  by  the  nonlinear  polarization  with  the 
same  time-dependence,  c.f.  Equations  85  and  67.  For  the  homogeneous  contribution  the 
additional  time  asymmetry  will  enhance  the  gain  in  the  positive  probe  delay  time.  This 
was  indeed  observed  in  the  transient  two-wave  mixing  experiments  in  gold  colloids  (see 
Figure  30).  Also  all  the  homogeneous  contributions  are  larger  than  the  corresponding 
inhomogeneous  counterparts.  This  is  the  reminiscence  of  the  additional  enhancement 
field  factor  in  the  homogeneous  parts. 
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6.0  A  ONE-DIMENSIONAL  CONJUGATED  POLYMER 


6.1  ONE-DIMENSIONAL  VIBRON-LATTICE  MODEL 

The  conjugated  polymers  with  n-electrons  have  been  demonstrated  to  have  high 
electrical  conductivity  and  unusual  nonlinear  optical  properties.  The  potential  high  speed 
(subpicoseconds)  and  large  spectral  range  are  very  attractive  for  nonlinear  optical 
processing.  The  nonlinearity  associated  with  ^-electrons  is  the  nonresonant  type,  which 
provides  the  large  spectral  range.  The  ir-electron  response  time  is  limited  by  the  one¬ 
dimensional  diffusion  time,  which  is  much  shorter  than  the  three-dimensional  one  and  is 
on  the  order  of  subpicoseconds.  The  band  gaps  will  be  determined  by  the  shortest 
wavelength  in  the  applications. 

Here  we  use  a  simple  and  tractable  one-dimensional  vibron  model  (Ref.  39)  is 
used  within  the  framework  of  lattice  dynamics  to  describe  conjugated  polymers  and  were 
able  to  predict  the  values  of  off-resonant  with  reasonable  accuracy. 

Particular  attention  will  be  given  to  Polyacetylene  (PA)  which  contains  a  high 

density  of  v-electrons  and  promises  an  ultrafast  response  time.  It  has  been  measured 

(Ref.  40)  to  have  =  10~^^  esu  at  off-resonance  spanning  from  1.17  eV  to  1.5  eV 
I  (3)  -9 

with  a  0.1  ps  response  time  and  have  x|  '  =  10  esu  at  the  narrow  two-photon 
resonance  peak  centered  at  0.91  eV  with  a  response  time  of  a  few  ps. 

The  v-electrons  of  PA  are  modeled  by  vibrons  with  harmonic  on-site  potential 
and  nearest  neighbor  interaction  between  them.  The  o-backbone  of  PA,  which  is  formed 
by  CH  groups,  is  modeled  by  a  one-dimensional  (ID)  lattice  within  the  harmonic  lattice 
approximation.  The  interaction  between  the  v-electrons  and  o-backbone  is  also  restricted 
to  nearest  neighbors.  To  be  more  specific,  the  model  Hamiltonian  of  PA  is  chosen  as 
follows: 


H  =  H^  +  Hg  +  Hjj^^  + 

In  Equation  86,  the  Hamiltonian  H^  for  v-electrons  (or  vibrons)  is  given  by 


(86) 


?  "e  I  "S  *  1  -  I '•(n,.)  (|,<|„ 

n  n  m 


(87a) 
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where  and  are  the  meiss,  velocity  and  displacement  of  the  n-th  ir-electrons.  The 
on-site  potential  v(q^)  for  n-th  vibron  is  assumed  to  be  harmonic  type  and  is  given  by 

v(q„)  -  \  «7b) 

where  the  eigen  frequency  is  closely  related  to  the  bandwidth  W  of  the  n-electron  band 
structure.  The  interaction  L(n,m)  among  the  ii-electrons  is  reduced  to 


J  L(n.»)  =  L  I  (q  j  +  q  j)  q„  (g7c) 

nm  n 


for  the  nearest  neighbor  interaction  among  them.  The  nearest  neighbor  interaction 
constant  is  related  to  the  hopping  integral  to  in  SSH  Hamiltonian  (Ref.  viz. 
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(88) 


where  t^  *  2.5  eV  for  PA. 

The  Hamiltonian  H^  for  o-backbone  (or  1-D  lattice)  is  given  by 

"o  -  5  «  I  “n  *  I  "  I  <V  «»> 

n  n 

where  M,  and  are  the  mass,  velocity  and  displacement  of  the  n-th  CH  group  (or 
the  n+l-th  lattice  site)j  and  K  is  the  effective  spring  constant  of  the  1-D  lattice.  The  K 
is  equal  to  21  eV/A  (from  Reference  41). 

The  interaction  Hamiltonian  H^^^  between  the  vibrons  (n-electrons)  and  the 
lattice  (o-backbone)  is  given  by 


‘  1  «Vl  -  “n)  Vl  ♦  <“n  -  Vl)  Vl> 


'n 


(90) 


105 


where  the  nearest  neighbor  interaction  is  assumed  and  the  vibron-lattice  coupling 
constant  x  is  related  to  the  electron-phonon  coupling  constant  a  in  SSH  Hamiltonian 
(Ref.  41),  viz. 


X  =  a  (91) 

where  a  =  4.1  eV/A  for  PA. 

The  response  to  the  external  optical  field  is  assumed  to  be  electric-dipole 
interaction,  viz. 

Vext  =  «lnE  <«> 

where  E  is  the  electric  field  of  the  external  optical  pulse. 

From  the  model  Hamiltonian,  the  equation  of  motion  for  PA  is  given  by 


"e  "n  *  Vo  On  *  *  "n.!) 

-  '("n.!  -  “n)  0,^1  *  (“n  '  Vl)  Op-il  -  -eE 


(93a) 


J  -^K(2U  -  U  1-U  ,)-^2x(q  ,-q  ,)q  -  0 
n  '  n  ihI  n-1'  ^n-1'  ^n 


(93b) 


which  are  respectively  for  the  x-electron  and  CH  group  at  the  n-th  site.  Since  M  »  m^ 
and  hence  0^«  q^  •  adiabatic  approximation  can  be  used  to  reduce  Equation  93b  as 
follows. 

KWUn  -  Un.l  '  Un-1>  *  '  On-l)  %  =  «  W 

Equation  94  can  be  integrated  and  have 


U  •  U  1  *  xr“  (Q  +  Cl  1 )  Q  +  U 
n  n+l  K  '^n  ^n-1'  ^n  c 
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where  is  an  integration  constant,  and  is  equal  to  zero  or  0.04  A  for  undimerized 
(unconjugated)  or  dimerized  (conjugated)  PA,  respectively. 

This  model  Hamiltonian  can  be  employed  to  design  and  study  other  conjugated 

polymers. 


6.2  THORD-ORDER  NONLINEAR  SUSCEPTIBILITY 

The  equation  of  motion.  Equation  93  within  the  adiabatic  approximation  for 
the  third-order  susceptibility  of  PA  will  be  solved  in  this  section. 

Substituting  Equation  95  into  Equation  93a,  the  equation  of  motion  for  the  n-th 
11-electron  is  obtained. 


■2^n  ♦  Vo^n  -  ‘Vl  *  Vl' 

2 

-  T  Vl  Vl>  Vl*  * 


(96) 


4x2 

where  is  the  effective  nonlinear  coefficient  for  tte  v-electrons,  and  this  particular 
combination  is  the  consequence  of  the  electron-photon  interaction. 

Equation  96  can  be  used  to  solve  for  the  linear  and  nonlinear  response  of  the 
n-electrons  in  the  presence  of  an  optical  field  at  w,  viz. 


E  =  E^e"^^  ♦  C.C. 


(97) 


For  a  I-D  lattice,  the  linear  displacement  of  the  n-th  v-electron  can  be  written 


as 


«  _  «  oUkan-«»t)  .  -  _ 

%  =  <Iie  +  c*c. 


(98) 


where  a  is  the  lattice  qwcing  (or  the  placing  between  CH  group).  The  amplitude  qj 
obeys  the  following  linearized  equation! 
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(99a) 


Inig(w^-  J)  -  4  (Lq+  xUq-  cos  ka)J 


which  yields 


eE 


0  -1kan 
e 


"i  “  ‘SIR 

with  the  linear  dispersion  relation  A^((i>)  at  u  is  given  by 
A^Cu))  =  ntgCwo-  -  4  (Lq+  xUq)  coska 

The  linear  dip>ole  moment  at  u  of  a  PA  chain  with  N-sites  is  given  by 


"  ik.n 


_  /  \  _  _  V  _  .iKan  o 

PlW  lie  = 


The  corresponding  polarization  without  local  field  correction  is  given  by 


where  the  linear  susceptibility  (w)  along  the  chain  is  given  by 


with  A  as  the  cross-section  of  the  PA  chain,  and  A  =  3.125  x  10”^^  cm^. 

The  measured  dielectric  constants  (Ref.  42)  at  9.1  GHz  for  trans-PA  are 
e|  =  5.7  and  Cj^  =  4.0,  respectively.  From  Equation  102  and  the  measured  e| 

Aj(w  ■  0)  =  1.62  X  10^  erg/cm^ 

Equating  Equations  100  and  103 


s  11.91  eV 


(99b) 


(100) 


(101a) 


(101b) 


(102) 


(103) 


(104) 
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which  is  comparable  to  the  bandwidth  w  =  10  -  12  eV. 


The  static  dispersion  relation  A|(u  =  0)  is  also  related  to  the  energy  gap  of  the 
ii-electrons»  viz. 

Pn  2 

Ai(«)  »  nig  (^)  (105) 

which  yields  Eg  =  2.78  eV  with  Equation  103  and  is  comparable  to  the  observed  optical 
gap  2.0  to  2.2  eV. 

The  displacement  of  the  1-D  lattice  up  to  the  third-order  at  u  can  be  written  as 

%  =  <«Il  +  «13>  <106) 

where  and  are  the  amplitudes  of  linear  and  nonlinear  displacements.  The  nonlinear 
amplitude  q^  obeys  the  following  equation, 


2  2 


)  -  4  (L^+  xUq)  coskal  q3 


(2  +  3  coska  +  cos2ka)  |qj|^  q^ 


(107a) 


where  q^  is  given  by  Equation  99b.  Then 


0.2  |E-|^  E 

q,  »  -  (2+3  coska  +  cos2ka)  - = — 5-^ - 


(107b) 


Then  the  third-order  nonlinear  polarization  at  u  without  local  field  correction 


is  given  by 

P3(«)*x<^^(«)|EJ2Ej 


(108) 


109 


where  the  third-order  susceptibility  along  the  chain  is  given  by 


v(3)/  \  e^  8(2-t-3  coska  +  cos2ka) 
“  aA  K  \  u  /  m2 


In  Equations  107  and  lOS*  the  condition  of  wavelength  being  long  compared  to  the  lattice 
spacing,  i.e.,  ka  «  1,  can  be  used.  In  Equation  109  every  parameter  is  fixed,  and  yields 


(u  =  0)  =  3.04  X  10“^°  esu 

(110a) 

x^^^  (hw  =  1  ev)  =  5.3  x  10“^®  esu 

(110b) 

x^^^  (hu  =  1.5  ev)  =  1.2  X  10“^  esu 

(110c) 

which  is  in  reasonable  agreement  with  the  measured  values 

1.17  eV  to  1.5  eV. 

=  10”^®  esu  spanning 

This  model  can  be  extended  to  include  the  coulomb  interaction  and  to  study 
other  conjugated  polymers. 
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7.0  CONCLUSIONS 


This  report  reviews  theoretical  and  experimental  research  into  the  active 
nonlinear  optical  processes  primarily  in  the  microparticle  composites  and  also  in 
conjugated  polymers. 

The  experimental  studies  done  at  the  University  of  California  at  Irvine  focused 
on  the  transient  effects  in  the  phase  conjugation  and  coherent  beam  combination  (or 
coherent  two-wave  mixing)  using  picosecond  pulses.  Beside  the  time-resolved 
experiments,  the  polarization  resolved  experiments  were  performed  to  demonstrate  the 
coherent  coupling  nature  in  the  two-wave  mixing  processes. 

The  theoretical  studies  emphasized  the  physical  mechanisms  for  the  nonlinear 
properties  of  the  metallic  microparticle  composites  and  the  conjugated  polymers.  The 
modeling  of  the  nonlinear  optical  processes,  viz.  the  phase  conjugation  and  coherent 
beam  combination,  was  carried  out  for  the  Kerr  media  and  the  metallic  microparticle 
composites.  The  theory  for  the  transient  effects  of  the  Kerr  medium  and  the 
microparticle  composites  is  discussed  using  either  the  time-dependent  third-order 
polarization  appearing  in  the  Maxwell's  equations  or  by  solving  the  Maxwell's  equations 
and  the  dynamic  equation  of  the  medium  simultaneously.  The  additional  asymmetry 
observed  with  respect  to  the  probe  delay  time  in  the  two-wave  mixing  experiments  for 
microparticle  composites  is  predicted  in  the  calculated  time-dependent  third-order 
polarization. 

In  the  past  much  attention  has  been  paid  to  the  study  of  linear  optical 
properties  of  small  metallic  particles.  The  mechanisms  of  the  surface  plasma  resonance 
and  of  the  broadening  of  the  linear  absorption  band  in  small  spherical  particles  have  been 
widely  investigated.  The  calculation  for  the  linear  optical  properties  was  first  performed 
assuming  independent  electrons  in  the  local  approximation  and  improvements  were  then 
brought  by  introducing  the  nonlocal  character  of  the  response  of  electrons.  Recently, 
the  nonlinear  optical  phenomena  present  in  such  particles  have  been  studied  using  optical 
phase  conjugation  in  metallic  colloids.  The  surface  plasma  resonance  occurs  for  this 
nonlinear  process  in  the  same  way  as  for  linear  properties. 
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A  solid-state  nonlocal  continuum  theory  has  been  developed  for  the  nonlinear 
electrodynamics  of  the  microparticle  composites.  The  calculation  of  the  Kerr-type 
polarizability  for  the  free  electron  part  of  small  metallic  spheres  using  the  quasi¬ 
hydrodynamic  theory  has  been  done  in  this  report.  It  is  understood  that  some  approxima¬ 
tions  have  to  be  introduced  in  order  to  perform  tractable  analytical  calculation.  The 
quasi-hydrodynamic  theory  which  includes  the  possibility  of  single  particle-collective 
mode  interactions  is  a  simpler  alternative  to  a  classical  kinetic  theory  or  a  quantum  self- 
consistent  field  approach.  However,  the  interband  transition  contribution  of  the 
microspheres  by  a  frequency-independent  background  dielectric  function  e^j  to 
renormalize  the  free  electrons  has  been  taken  into  account.  The  free  path  effect 
approximation  is  used  for  the  surface-scattering-contributed  additional  damping  inside 
the  microparticles.  Silver  spheres  can  be  described  accurately  by  this  model.  That  is, 
the  calculated  third-order  susceptibility  from  dressed  free  electrons  and  the  measured 
third-order  susceptibility  are  comparable  in  silver  spheres,  which  are  both  10~  esu.  Gold 
spheres  can  be  explained  by  this  model  with  two  fitting  parameters,  i.e.,  the  background 
dielectric  constant  and  the  damping  factor,  for  their  free  electron  contribution. 

However,  the  calculated  third-order  susceptibility  from  free  electrons  in  gold  spheres  is 
much  too  small  compared  to  the  measured  one,  which  indicates  the  important  interband 
contribution  in  the  optical  nonlinearity  for  gold. 

A  semi-classical  lattice-dynamic  theory  was  derived  to  predict  the  large  non¬ 
resonant  optical  nonlinearity  below  the  band  gap  for  1-D  conjugated  polymers.  The 
characteristic  electron-photon  interaction  driven  electronic  nonlinearity  is  clearly 
identified  in  the  formalism  of  the  Kerr-type  polarizability  calculation  for  conjugated 
polymers. 

The  theoretical  studies  concentrated  on  the  electronic  nonlinearity  of  the 
microparticle  composites  which  is  caused  by  free  (or  conduction)  electrons.  This 
nonlinearity  was  also  the  focus  of  the  short  (picosecond)  pulse  experiments.  The  findings 
will  have  applications  in  high  bit-rate  transmission  and  communications. 

The  electronic  nonlinearity  involved  in  the  microparticle  composites  is  a 
resonant  type,  which  results  from  the  local  field  enhancement  of  the  excited  surface 
plasma  modes  and  always  is  accompanied  by  an  enhanced  linear  absorption  for  lossy 
materials.  Since  it  is  a  resonant  nonlinearity,  there  is  limited  spectral  range  for  the 


112 


nonlinearity,  which  is  0.05  ym  or  0.1  ym  for  silver  or  gold  particles.  Because  of  the 
enhancement  in  the  linear  absorption,  the  nonlinear  efficiency  must  be  optimized  with 
some  compromise  in  the  absorption  loss,  which  can  be  achieved  by  the  small  volume 
fractions  or  thin  samples.  The  optimal  thicknesses  of  the  samples  are  a  few  milli¬ 
meters.  Small  volume  fraction  is  also  needed  to  eliminate  the  unwanted  scattering  loss, 
which  has  to  be  less  than  10  .  Because  of  the  potential  excessive  Mie  scattering  from 
large  particles,  the  particle  size  below  100  A  must  also  be  restricted. 

In  the  experimental  studies,  the  microparticle  composites  were  used  for 
picosecond  phase  conjugation  as  well  as  weak  probe  signal  amplification.  The  time- 
resolved  and/or  polarization  resolved  degenerate  phase  conjugation  and  two-wave  mixing 
experiments  were  performed  in  the  microparticle  composites  with  a  volume  fraction 
between  4  x  10~^  and  10“^  with  particle  radii  <  120  A.  The  degradation  of  the  conjugated 
or  amplified  images  from  the  scattering  by  particles  is  negligible  with  aforementioned 
limits  of  the  sizes  and  volume  fractions.  However,  the  absorption  from  the  particles  is 
influential  under  all  circumstances,  which  limits  the  sample  thicknesses  to  1  to  2  mm. 

The  relaxation  time  of  the  free  (or  conduction)  electron  part  of  electronic 
nonlinetu'ity  in  the  metallic  microparticle  composites  is  on  the  order  of  tens  of  femto¬ 
seconds  according  to  the  quasi-hydrodynamic  theory.  However,  in  the  case  of  the  gold 
the  bound  (or  valence)  electron  contributed  electronic  nonlinearity  contributed  by  the 
interband  transition  will  have  much  longer  relaxation  time.  For  the  two-wave  mixing 
experiments  the  relaxation  time  was  on  the  order  of  tens  of  picoseconds. 

The  conjugated  polymers,  which  have  nonresonant  nonlinearity  below  the  band 
gap  with  potential  subpicosecond  relaxation  time  because  of  fast  one-dimensional 
diffusion  characteristics  were  also  investigated.  Because  the  nonlinearity  is  a 
nonresonant  type,  the  spectral  range  of  the  conjugated  polymer  is  large.  However,  there 
is  the  shortest  wavelength  limit  in  the  nonlinear  conjugated  polymer  applications,  which 
is  set  by  the  band  gap.  The  limits  on  the  band  gap  of  the  polymers  can  be  easily  tailored 
by  the  chemical  engineering,  e.g.,  by  doping  or  using  a  different  backbone  element.  The 
PAs,  considered  here  to  have  a  band  gap  of  2  eV,  are  suitable  for  IR  applications.  The 
polysilanes  and  polygermanes  have  band  gaps  in  the  UV  region  and  can  be  useful  for  all 
visible  applications. 
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APPENDIX  A 


To  solve  Equation  7a  for  a  sphere,  expand  the  electron  density  in  the  following 

series 

nj(r,«j)  =  i:A^„(qja)  (A.l) 

where  a  is  the  radius  of  the  sphere  and  ip(z)  is  the  n^^  order  modified  spherical  Bessel 
function  of  qj  =  q((tij).  In  spherical  coordinates  r  =  (r.e,^)  and  Yj^(e,^)  is  the 
spherical  harmonic  of  order  (l,m).  Since  the  i«)|  laser  beam  is  elliptically  polarized, 
nj^(r,iitj)  consists  of  the  1  =  1,  m  =  ±  1  term.  The  W2  beam  is  linearly  polarized  in  the 
02  -  (xcose2  zs1ne2)  direction  and  n|(r,ui2)  contains  the  1  =  I,  m  =  ±1  (see  Fig.  I). 
Applying  the  boundary  conditions  at  the  surface  yields  the  following  solutions  for  the 
first-order  electron  density  at  the  two  Uiser  frequencies 

n^(r,<m)  =  N(r,«j)Aj(nCose  +  sine  s1n*)/Dj^(w^)  (A.2a) 

®b 

ni(»*i«2)  *  4^  N(r,«2)A2(s1ne2COS0  +  cose2Slnecose)/D^(w2)  (A.2b) 

The  boundary  conditions  yield  the  following  expression  for  N(r,ii)j) 

where  i'j(z)  implying  the  derivative  with  respect  to  z,  and  the  denominator  Dj(oij)  is  given 
by 


DlCwj)  •  ^^*^**^ 

2 

with  fij  =  **j(*"j  ***  V  la  assumed  to  be  frequency  independent.  The  real  part 

of  this  important  denominator  vanishes  at  w^.  This  denominator  dominates  the 
electrodynamics  of  microparticle  composites  whenever  the  surface  (or  Frohlich)  dipole 
mode  is  excited  by  Uj  s 

The-first  order  electric  field  Ej^(r,uj)  is  given  by 
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Ei^(r,uii)  =  Aj^(nCose  +  sine  s1ne)/D^(u^)  (A.3a) 

=  ^y.{r,a2)  A2(s1n02  cose  +  cose2Slnecose)/D^(u2)  (A.3b) 

Ele(^.wi)  =  Aj(-ns1ne  +  cose  s1ne)/Dj(«j)  (A.3c) 

Eie(r,«2)  ®  ig(r,w2)A2(cose2Cosecose  -  s1ne2Slne)/Dj^(u2)  (A.3d) 

Ej^(r,«i)  -  rQ(r,a»j)Aj  cose/Dj(«j)  (A.3e) 

Ei^(r,«2)  *  -rg(r,w2)A2COse2Slne/0j(«2)  (A.3f) 

where 

z^(r,iaj)  »  3ej,(«j)  [l-lJ(qjr)/fljl|(qja)J  (A.3g) 

re(»‘»"j)  “  Il-1j(qjr)/0jqjrlJ(qja)I  (A.3h) 

The  first  order  electron  drift  velocity  Vj(r»wj)  is 

Vi^(r,«i)  -  J^(r,«j)VjQ(ncose  +  sine  s1ne)/Dj(caj)  (A.<^a) 
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Vir(r,w2)  =  z^(r,u2)V2Q(sin02Cos0  +  cos02Sln0cos^)/Dj(u)2)  (A.4b) 


VigCr.ui)  =  Zg(r,u)^)V^Q(-nSln0  +  cos0  s1n^)/Dj(uj) 


(A.4c) 


=  rg(r,w2)V2Q(cos02COS0COS#  -  sin02Sln0)/Dj(a>2)  (A.4d) 


=  ig(r,u2^)Vj^QCOS4/Dj^(wj^) 


(A.4e) 


Vi^(r,w2)  =  -rg(r,«2)V2QCOS02Sln*/Dj(w2) 


(A.4f) 


where  Vj^,  =  -ieAj/m(«j  +  iv)  and  the  variously  are  defined  as 


r^(r,«j)  »  3e^(wj)[l  -  1'i(qjr)/i'j(qja)l 


(A.4g) 


20(r.wj)  -  3e^(«j)Il  -  1j(qjr)/1'j(qja)qjr) 


(A.4h) 


Outside  the  microparticle,  the  first-order  electric  field  is  dipolar  in  nature.  In 
the  near  (static)  zone,  the  electrostatic  potential  of  a  sphere  is  given  by 


•^(r.w^)  -  r  + 


2 —  )  A,e, 


<A.5) 


where  the  angular  functions  e|  are  given  by 
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s  ncose  -t-  s1nes1n^ 


62  *  s1n62Cose  cose2Slnes1n^ 


In  Equation  A.3»  the  first-order  dipole  moment  at  U|  is  given  by 


..(1)/  \  -3  o(l),  . 

u'  '(01^)  *  3-  a  P  («^) 


where  a  is  the  radius  of  the  sphere  and  the  first-order  polarization  is  given  by 


-  * 


(A.6a) 

(A.6b) 

(A.7) 

(A.8) 
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APPENDIX  B 


First,  the  different  components  S*  (r)  of  the  driver  term  S“(r)  in 
Equations  2.11  depend  only  weakly  on  frequeilf^  and  do  not  reflect  the  surface  mode 
resonances  of  the  microparticles,  and  are  defined  below: 


S*  (r)  =  S*  (r)  -3lrQ(r.«j)r^(r,±w2)  +  rQ(r.±u2)E^(r.uj) iL'a/r 
00  20 


(B.la) 


(r)  *  N^(r,«j)E^(r.±«2 


)  +  Nj(r,±a»2)j:^(r.u)j) 


+  Nj(r,wj)E^(r,±u2)uj/(±w2'*‘1v)+Nj(r,±u)2)E^(r,u)j)u2/(wj^+iv) 


-  L■l2^(r.«^)Ey'(r.+u.2)-^Q(r,«^)E^(r.±w2)a/r)l 


-  L-lE^(r,±w2)E'(r.«j)-EQ(r,±u)2)E^(r.wj)a/r] 

+  2/3(6/«pa)^(Nj(r,«^)Nj'(r,±o)2)+Nj(r.±«2)Ni'(r.u,j)l) 

(B.lb) 


S*  (r)  »  S"  (r)  +  2L*  J  E  (r,±  «-)  (B.lc) 

®20  ’21  rate  i. 

Sj  (r)  -  I  S^(r,wj,±w2)  -  S^(r,±w2.»i)  L"  ^  Eg(r,±w2) 

(B.ld) 


(B.le) 
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(B.lf) 


^022^’^^  '  ^020^’^^  r  -  “2^ 

^♦21^*^^  =  S^Cr.w^.tu,^)  +  S^(r.±u2.wj) 

^♦22^^^^  *  \  ^^(’^•“1-  “2^  "  f  ^4,^’^*  “  “2*“l)  r  “  “2^ 

and 

^♦^’^•“l*““2  ^  *  Nj(r,uj^)rQ(r,±w2)  +  Nj(r,wj^)ZQ(r,±u)2)u^/(±u)2+  iv) 

-2^(«".«l)2'0(r.±«2)L" 

+  2/3(6/Mpa)^Nj(r,wj)N2(r,±M2)a/r  (B.li) 

where  the  dimensionless  linear  electron  density  Nj(r,«i)j)  is  given  in  Equation  A.2c  and  the 
function  is  defined  by 

L~  *  <dp[  (^2'*’  ^  (B.2a) 

and  the  primed  functions  are  defined  as 

F'(r,wj)  «  a  aF(r,uj)/8r  (B.2b) 


(B.lg) 

(B.lh) 
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Then,  to  solve  Equations  2.8  for  a  sphere,  expand  n2(r,ui^)  in  terms  of  sherical 
harmonics  as  n|(r,uj). 


n2(r,u.p  =  z 


+  Jdr' [65(q+lr.r  )vS^(r')l/6^ 


(B.3a) 


where  the  integral  is  an  inhomogeneous  term,  the  coefficients  are  determined  by  the 
boundary  conditions.  The  Green's  function  G^(q^|r,r')  is  given  by 


G5(q^|r,r')  =  q+i:1i(qi«^<)Ki(q+r^)Yi/(e',^')Y^Je,^) 


(B.3b) 


where  *  implies  a  complex  conjugate.  The  sum  and  difference-frequency  fields  obey  the 
Laplace  equation 


(r,«+) 


(B.4a) 


outside  the  microparticle,  and  the  Poisson  equation 


v^*2(r,w^)  -  4x602 (r,a)^)/ejj 


(B.4b) 


within  it.  Equations  B.2  will  be  solved  by  expanding  in  spherical  harmonics.  Thus, 
outside  the  sphere, 
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(B.5a) 


♦2(r.co^)  =  2C^^(q^a) 


The  electrostatic  potential  inside  the  microparticle  is  driven  by  the  second-order  charge 
density,  n2(r,up.  Thus,  within  this  region,  the  second-order  electrostatic  potential  is 
written  as 


(B.5b) 

where  GQ(r,r')  is  the  Green's  function  for  the  Laplace  equation.  The  coefficients  for  the 
various  multipole  fields  and  are  determined  by  the  electrostatic  boundary 
conditions  at  the  surface  of  the  microsphere.  The  first-order  fields  exhibit  a  dipolar 
character  outside  the  microparticle,  therefore,  the  electrostatic  fields  consist  of  1  =  1, 
m  =  ±1  at  wj,  and  i  =  1  m  =  ±1  at  W2«  Thus,  the  second-order  fields  should  have  monopole 
and  quadrupole  characters.  The  second-order  electrostatic  field  should  consist  of  1  =  0 
and  the  1  s  2,  m  =  (0,±l,  ±2)  terms.  However,  a  monopole  field  vanishes  by  charge 
conservation.  Hence,  only  the  coefficients  for  the  quandrupole  contribute,  i.e.,  D20 
(C20).  *^2±1  ^^2±1^  *^2±2^^2±2^  non-zero.  Finally,  the  second-order  drift 

velocity,  V2(r,Uj)  can  be  expressed  in  terms  of  n2<r,up  and  *2(r,(i)^)  via  Euler's 
Equation  and  given  by 


V2(r.«+)  -  v"(r,«^)  +  VgCr.u^) 


(B.6) 


where 


and 


V2(r.a)^)  =  -  ^  n^(r,(i.^)v^(r.±u.2)  -  —  n^(r.±a)2)vj(r.a)j) 
o  o 


(B.6b) 


Equations  B.l  through  B.3,  along  with  the  boundary  conditions  on  the  drift  velocity  and 
the  electrostatic  potential,  completely  determine  the  unknown  coefficients  B|j^,  and 
It  is  a  tedious  but  a  straightforward  task  to  obtain  the  following  expression  for  the 
density  n2(r,u)+); 


n2(r.«+)  =  r+{(4Tt/9)^/^n2^Q0(r.u.^)nsin82YQQ(e) 

+  (16iT/45)^/^n2^20''''*“+)''^’"V20(®) 

+  1  (2Ti/15)^/^n2^21+^’"»"±)^^"®2^^21^®**^ 

-1(2Tr/l5)^^^n2^22^’^*“±)  cos02lY22(e.*)  -  Y2_2(0.4»)l> 

(B.7a) 


where 


D(«i)D(w2) 


(B.7b) 
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2  1/2 

where  =  lup  -  u+(u+  +  1v)l  /e  .  The  quantity  r+,  which  has  units  of  density, 
exhibits  resonant  behavior  if  and/or  ai2  lie  near  the  surface  dipole  mode  frequency 
This  feature  of  the  second-order  density  response  follows  directly  from  the  drive  term, 
S^(r).  An  examination  of  Equation  B.7a  reveals  the  important  role  of  the  initial 
geometry  as  manifested  in  the  angle  62  and  the  initial  polarization  as  manifested  in  n. 

The  dimensionless  functions  in  monopole  and  quadrupole  radial  density 

n2,oo^''»“+)  and  n2^2m^''*“±^ 


where  the  denominator  D2(«>P>  which  appears  only  in  the  quadrupole  radial  density  terms 
in  Equation  B.7d,  is  given  by 


D2K)  =  -  1)  +  (q^  -  2i2(q+a)/q+ai^(q+a)l  (B.8) 

2 

where  =  “+(«+  +  •  Note  that  D2(up  is  a  function  only  of  w^,  not  U|  or  102 

individually.  This  feature  follows  because  D2<wP  reflects  the  possibility  of  exciting  a 
surface  mode  if  coincides  with  a  quadrupole  Frohlich  surface  resonant  frequency 

In  Equations  B.7,  the  functions  3£,j,(r,up  are  defined  as 


where  (l,m)  >  (0,0),  (2,0)  and  (2,1),  two  functions  Jqos  and  ^2111®*'®  by 
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(B.9b) 


"  (3eh(“+)  +  2ejj)  (Q+-  1)  J2in^**“±^ 
+3e^(u.+)J2ms(a.«+) 

with 

''2l«S<''*“t>  “  T  /  '^‘"■''"‘'■'2lll<'''-“±> 

a  0 

The  functions  j£m^(r»«4)  and  given  by 

Joo*('-.»+)  -  (qt/«^)/o'^<lr‘r'^1'g(cijr')S  (r’) 

Joo*(''.“±)  -  (r') 

00 


(B.9c) 

(B.9d) 

(B.9e) 

(B.10a) 

(B.lOb) 
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j2in^(r.««»P  =  (q+/a^)Jo'^dr'r'^i'2(q+r*)S;^(r') 

+  (l/a^)Jo''dr'r'i2(q+r')S;  (r')  (B.lOc) 

“  ^  -  ®2m 

and 

J2in^(r.u±)  =  (q+/a^)//dr'r'^k2(q+r')S;^^(r') 

+  (l/a^)//dr'r'k2(q+r')S*  (r')  (B.lOd) 

~  2in 

where  the  (r)  are  given  in  Equations  B.l. 

“«n 

The  next  focus  is  on  the  electric  component  of  the  second-order  electric  fields 
oscillating  at  the  sum  and  difference-frequencies  in  the  vicinity  of  the  microparticle. 
Inside  the  particle,  the  various  components  of  the  electric  field  are  given  by 

®  E2q{(4w/9)  ^ 

-  (2x/15)^^^cos02nlY2j  (e,4)-Y2_l 
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+  1(2i./15)^/2s1n02[Y2i{e.*)  +  Y^.^Cb.^)] 

-  1  (2ir/15)^/^cose2lY22(e,^)  -  Y22(e.«)] 

(B.lla) 

^20^'^*“±^  “  ^20  ^  $11162^  —  Y2q(0)  ^0^20^'^^ 

-  (2ii/15)  ^  €0502^  “  Y2_2(B»t)]  ig  21-^*^^ 

+  1  (2i/15)^/2  |_  iY^^(e,^)  +  Y2.i(e.«)l  2e.21+<^) 

-  1  (2v/15)  ^  COS02  [^22^®*^^  ”  ^2-2^®*^^^  ^0  22^*^^^ 

(B.llb) 

E2^(**»“±)  ®  E2q{-(2»/15)  ^  cos  02n  3^  IY2]^(0»^)  -  Y2_j^(0t^)]  ^0^21-^’’*^ 

+  1  (2,/15)1/2  s1ne2^|^  IY2i(0.O  ^  Y2.i(0.Ol  i;.21^(’^) 

-  1  (2,/15)1/2  cos02^|^  IY22(e.4)  -  Y2.2(e.01 

(B.llc) 
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where 


^20  ~ 


(B.12) 


Thus,  the  second-order  electric  fields  have  r,e  and  it-components:  This 
complex  polarization  structure  arises  from  the  fact  that  the  incident  fields  for  the 
noncollinear  configuration  break  the  cylindrical  symmetry  of  the  system,  a  point  which  is 
underscored  by  the  fact  that  the  ♦-component  of  the  electric  field  vanishes,  viz,  the 
cylindrical  symmetry  is  restored,  if  the  two  laser  beams  are  collinear  (i.e..  Go  =  ii/2)  and 
copolarized.  The  functions  z  .-(r)  are  dimensionless  field  strengths  which  are  defined 

Cl  g  Jblll 

by: 


^r,2in^^)  "  (B.13b) 

'^9.2m<''>  ■  +  ■:ai6(''.“±)i/3  (B.13c) 

where 

Q+  •  ^®b^ 

The  functions  and  ^2inG^’^*“+^  defined  as 

i2ii|(»*»“±)  *  3/5I(a/r)^J2„^(r,»^)-32„(a,«^)l3(q^r)/lJ(q^a)D2K)] 

(B.14a) 
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(B.14b) 


^2inG^’^*“±^  “  3/5lq+rJ2||^(r,w^) 

-  J2„(a.u+){1j(q+r)-ij(q^a)r/a}/i^(q^a)D2(u)+)l 


with  3qq5>  32mS  ^2mS  g^ven  in  Equations  B.7b,  B.7c  and  B.7e  and 


*^2in6 


(B.15) 


The  electric  field  inside  the  microparticle  has  two  possible  frequencies  where 
a  surface  mode  resonance  can  be  excited.  Specifically,  the  presence  of  the  factor 
implies  possible  resonant  enhancement  of  the  field  if  (i)|  and/or  ti>2  coincide  with  the 
surface  dipolar  mode.  Further,  the  D2(<o^)  resonant  denominator  in  Equations  B.14 
implies  that  the  quadrupole  surface  mode  at  u>^  can  also  amplify  the  field.  In  actual 
practice,  all  three  frequencies  (w|,  ^2  and  cannot  coincide  with  a  Frohlich  resonance 
for  the  second-order  processes  discussed  here. 

Outside  the  microparticle,  the  electric  fields  are  quadrupolar  in  nature.  In  the 
near  (static)  zone,  the  electrostatic  potential  of  a  sphere  is  given  by 


»2(r.«j.)  -  l/3{(16./45)l/^y20(e)Qjj 

-  (2./15)*/^  |V2l(e,*)  -  Y2.i(e,,))q^j)/r^  (B.16) 


The  components  of  the  quadrupole  moment,  in  terms  of  the  dimensionless  quantities 
U2j  are  given  by 
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(B.17a) 


^zz  “  -l6iier^s1ne2U2Qa®/3(3ej^(u^)  2e|j) 

=  -4iter^cos02U2ja^/(3e^(«^)  +  2£|j)  (B.17b) 

where  is  the  dimensionless  function  in  the  quadrupole  moments  and  given  by 

“an  ■  7  j"?  "2.ai(''-“t> 

cl 

=  lJ2in<**“±)^3(V^/‘^2^“±^  (B.lSb) 

As  expected,  the  quadrupole  moments  and  the  electrostatic  field  outside  the  micropar* 
tide  are  also  surface  enhanced  if  either  of  the  frequencies  or  U2  coincides  with  the 
surface  dipole  resonance  (from  or  with  the  surface  quadrupole  resonance  from 
D2(u^))« 

H 

The  second-order  electron  drift  velocity  V2  (r.w^)  in  Equation  B.6a  is  given 
by 

V2r('^.“±)  •  V20  [(4»/9)^/^ 

+  (l6./45)^/^1ne2nY2o(e)i;  20<’") 

-(2w/15)^/^cose2n(Y2j(e,^)  -  Y2_i(e,4))iJ  2i-(’^) 

+  1  (2»/15)^/^  s1ne2(Y2i(e,4)  +  Y2_i(e,#)) 

-  1  (2i,/15)^/2  cose2(Y22(e,0)  -  Y2.2(e.O)  2^,22^'*)^ 


(B.19a) 


vJeCr.u,^)  =  V2o((16„/45)^/2  i_  2e.2o('‘^ 

-  (2ir/15)  ^  C0S62n  (^21^®*^^  ~  ^2-1^®*^^^  ^6  21-^*^^ 

+  1  (2h/15)^/2  |_  (Y2i(e.4.)  +  Y2.i(e.^)) 

-  1  (2,r/15)^/2  00502  I3  ^e.22('*)> 

(B.19b) 

''2«^'^*“±^  “  '^20^"  (2ir/15)^^  ^°^®2'*  !♦  “  ^2-1^®**^^  ^♦,21-^’^^ 

+  1  (2»/15)  ^  ®^''®2  sTnefe  ^^21^®**^  ■*■  ^2-1^®****^^  ^♦,21+^’"^ 

(2«/15)  ^  ^®®®2  TTn0  (^22^®*^^  ”  ^2-2^®*^^^  ^♦,22^'^^^ 

(B.19c) 

where 

V20  “  -  16^  V"o(“+  ^'»)  <B.20) 

And  the  dimensionless  function  £*  ,-(»*)  is  given  by 

Q  f  jmi 

(B.21a) 

"J.2ii<»‘)  •  ^201  12<V)/^2<V)°2<"±)  - 

+  4,2m^^)  (B.21b) 
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(B.21c) 


q+a 

'i.21<'')  ‘  “21<''>  *  4t  5;,2i('')  <B-21<i) 

q+a 

where  Jqqj,  Jqq,  J2uj,  Dg,  ^0^2m  21  are  defined  previously  in  Equations  B.7, 

B.8  and  B.9,  E“  is  the  dimensionless  field  strengths  given  in  Equations 
B«13y  q'l'p  *  Up/6  t  and 

4<'')  •  "2.2ll<'’'“±>/V  «i.2l>(''>  ‘'tF''‘’±  <B-22) 
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APPENDIX  C 


To  solve  Equations  2.12  for  a  sphere,  expand  n^(r,iDj)  in  terms  of  spherical 
harmonics  as  n£(r,Uj)  and  n2(r,u^). 


n3(r..,j)  .  I  8„(qja)1,(qjr)V,^(e.*) 


+  /dr' [6j(qjr,r')7*T(r'>j)|/6^ 


(C.la) 


where  the  integral  is  an  inhomogeneous  term  and  the  coefficients  are  determined  by 
the  boundary  conditions.  The  Green's  function  G5(qj|r,r')  is  given  by 


Gs(qjlr.r')  =  (qjr^)Y^/(e' )Y^  Je,^) 


(C.lb) 


where  *  implies  a  complex  conjugate.  The  third-order  fields  at  wj  obey  the  Laplace 
equation 


vSjCr.uj)  =  0 


(C.2a) 


outside  the  microparticle,  and  Poisson  equation 


V  *  4nen3(r,uj)/ejj 


(C.2b) 
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within  it.  Equations  C.2  will  be  solved  by  expanding  in  spherical  harmonics.  Thus, 
outside  the  sphere,  we  have 


(C.3a) 


The  electrostatic  potential  inside  the  microparticle  is  driven  by  the  third-order  charge 
density,  Thus,  within  this  region,  the  third-order  electrostatic  potential  is 

written  as 


♦3(r.uj)  =  -  (4iie/ejj)Jdr'G^(r.r' )n3(r' 

(C.3b) 

where  G^Crjr')  is  the  Green's  function  for  the  Laplace  equation.  The  coefficients  for  the 
various  multipole  fields  and  Fjj^  are  determined  by  the  electrostatic  boundary 
conditions  at  the  surface  of  the  microsphere.  The  first-order  fields  exhibit  a  dipolar 
character  outside  the  microparticle,  viz,  the  electrostatic  fields  consist  of  1  =  1,  m  =  ±1 
at  Up  and  &  =  1  m  =  ±1  at  U2.  The  second-order  fields  have  monopole  and  quadrupole 
characters.  The  second-order  electrostatic  fields  consist  of  1  =  0  and  the  1  =  2, 
m  =  (0,±1,±2)  terms.  However,  a  monopole  field  vanishes  by  charge  conservation.  Hence, 
only  the  quadrupole  contributions  are  nonzero.  Then,  the  third-order  field  should  have 
dipole  and  octupole  characters,  since  the  third-order  electrostatic  fields  consist  of  t  =  1, 
m  =  0, 1  1  and  4  =  3,  m  =  0,  m  =  ±l,m  =  ±2,  m  =  ±3  terms.  Finally,  the  third-order 
drift  velocity,  V3(r,Uj)  can  be  expressed  in  terms  of  n3(r,Uj)  ancto^Cr.uj)  via  Eu'er's 
Equation  and  given  by 


V3(r,uj)  =  V3(r,uj)  +  vjcr.uj) 


(C.4) 


where 
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(C.5a) 


VaCr.uj)  =  -1[T(r;«j)  -  (0^6/111)7* (r.wj)  +  n3(r.a)j) ]/nQ(uj+  iv) 


V3(r.u»j)  = 


-  ^  nj^(r,u»2)v2(r,u_)  -  ^  n2(r,w_)Vj(r,u2) 

0  0 


~  n  ”  n  *^2^******+^  '^]^(*‘»““2^ 


(C.5b) 


Equations  C.  1  through  C.3,  along  with  the  boundary  conditions  on  the  drift  velocity  and 
the  electrostatic  potential,  completely  determine  the  unknown  coefficients  E|j^,  F|j^, 
and  G|^.  It  is  a  tedious,  but  a  straightforward  task  to  obtain  the  following  expression  for 
the  density  n3(r,Uj): 


n3(r,o)j)  = 


(C.6a) 


where 


^3J  =  e5Ajqja^/4»B^B^Dj(uj) 


(C.6b) 
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2  1/2 

where  Qj  =  [up  -  (Oj(uj  +  Im)]  '  /b  .  The  third-order  density  exhibits  resonant 
behavior  if  a)|  and/or  U2  lie  near  the  surface  dipole  mode  frequency  Ug.  This  feature  of  the 
third-order  density  response  follows  directly  from  the  drive  term,  T(r;u)j).  In  r 

Equation  C.6a,  only  the  dipole  terms  which  contribute  to  the  third-order  susceptibility 
are  written  out  more  explicitly.  , 

The  dimensionless  functions  in  dipole  radial  density  n^^  |Q(r,a)j)  and  n^j 
are  given  by 


(C.6c) 

where  the  denominator  D|(uj),  which  appears  in  the  dipole  radial  density  terms  in 
Equation  C.6c,  is  given  by  Equation  A.2d. 

In  Equations  C.6,  the  functions  L|^jjj(r,Wj>  are  defined  as 

where  (4,m)  =  (1,0)  and  (1,1),  function  is  defined  by 

*  ‘b>  (“j- 

+2.^(.j)Lj^(a,«j)  (C.7b) 

with 
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(C.7c) 


'■l«st'--“j)  =  ^  /  -“j)  <C-7<« 

&  o 

The  functions  aj^j^^(r,wj)  and  are  given  by 


+  {l/«^);r*<lr'r'kj(qjr')  (T^  ('•■:i«i)l 

^  T  m  ^  ’1  m  ^ 


(C.8b) 


CXitside  the  microparticle,  the  electric  fields  are  dipole  and  octupole  in 
nature.  In  the  near  (static)  zone,  the  dipole  part  of  electrostatic  pc  '^ential  of  a  sphere  is 
given  by 
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+  i  (2ir/3)^/2  [Yii(0.^)  -  Yj_j(6.^)lu{3)(a,  ) 


(C.9) 


(3) 

The  components  of  the  third-order  dipole  moment,  U|||  ' ,  in  terms  of  the  dimensionless 
quantities  are  given  by 


2  ^ 

4ner,.  y  P 

3Ji£l 


(3) 


Qj  (2eh(“j)  + 


(C.lOa) 


/ox 

where  Pj  is  the  dimensionless  function  in  the  third-order  dipole  moments  and 

given  by 


(C.lOb) 


(C.10c) 

As  expected,  the  third-order  dipole  moments  and  the  third-order  electrostatic  field 
outside  the  microparticle  are  also  surface  enhanced  if  either  of  the  frequencies  U|  or  tt>2 
coincides  with  the  surface  dipole  resonance. 

In  Equation  C.iOc,  the  first  and  second  terms  are  contributed  by  the 
homogeneous  cuid  inhomogeneous  solution  of  the  third-order  density  given  in  Equation 
C.6a.  The  homogeneous  part  has  an  additional  local  field  enhancement  factor  as 
compared  to  the  inhomogeneous  part  for  the  third-order  polarization  in  the  frequeixry 
domain.  This  difference  between  the  homogeneous  and  inhomogeneous  parts  will  appear 
in  the  time-dependent  third-order  polarization  to  be  discussed  in  Section  3.3.2. 
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The  three  factors  of  the  local  field  enhancement  are  due  to  the  third-order 
nonlinear  processed  in  the  small  spheres  and  the  additional  fourth  factor  appearing  in  the 
homogeneous  part  occurring  through  the  boundary  condition  of  the  sphere  for  the 
homogeneous  part  of  the  third-order  density. 


,(3) 

in 


Finally,  the  nonlinear  third-order  polarization  of  a  single  sphere  is  given  by 


*1 


'1 


(C.lla) 


where 


(3)  _  H 

^0  4Tr(2ejj(uj)  +  ejj) 


2  2 
e  a 

n.2,4 

m  6 


(C.llb) 
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